

DOE/NASA/0084-1 

NASA CR-1 68292 nasa-CR- 168292 

19840011117 


Optimization Methods Appiied to 
Hybrid Vehicle Design 


John F. Donoghue and James H. Burghart 
Department of Electrical Engineering 
Cleveland State University 


June 1983 


LIHRMY COPY 


IV;At< 8 1984 

UNGLEY RESEARCH CENTER 
IIBRARY, NASA 
HAy:^TOy, VIRGINIA 


Prepared for 

NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 
Lewis Research Center 
Under Grant NAG 3-84 


for 

U.S. DEPARTMENT OF ENERGY 
Conservation and Renewable Energy 
Office of Vehicle and Engine R&D 



DISCLAIMER 


This report was prepared as an account of work sponsored by an agency 
of the United States Government. Neither the United States Government 
nor any agency thereof, nor any of their employees, makes any warranty, 
express or implied, or assumes any legal liability or responsibility for the 
accuracy, completeness, or usefulness of any information, apparatus, 
product, or process disclosed, or represents that its use would not 
infringe privately owned rights. Reference herein to any specific 
commercial product, process, or service by trade name, trademark, 
manufacturer, or otherwise, does not necessarily constitute or Imply its 
endorsement, recommendation, or favoring by the United States 
Government or any agency thereof. The views and opinions of authors 
expressed herein do not necessarily state or reflect those of the United 
States Government or any agency thereof. 


Printed in the United States of America 


Available from 

National Technical Information Service 
U.S. Department of Commerce 
5285 Port Royal Road 
Springfield, VA 22161 

NTIS price codesi 

Printed copy; AO 9 
Microfiche copy: A01 


“iCodes are used for pricing all publications. The code is determined by 
the number of pages in the publication. Information pertaining to the 
pricing codes can be found in the current issues of the following 
publications, which are generally available In most libraries: Energy 
Research Abstracts (ERA); Government Reports Announcements and Index 
{GRA and I); Scientific and Technical Abstract Reports (STAR); and 
publication, NTIS-PR-360 available from NTIS at the above address. 



L KlM/lMHaH' 
Ul^ruHY IH/il-L 


iSSUt H HHUb CU IS, , oin^n 

/I WAS u26« 168292 MAG3-84 uL^hUji'" / 


loy 


1/ UUWUUnUi^? 


I L*-J llrUi I’J. 


lU 


d^^SlQn ILbr- rlndl Kepor 


, ? D/ DUr\^J5nr1i\ I 


r I I <d\ V'J- Lc^. ^ - 

I I T ! ?. !-r T f-^ !~"Ar>s UP 

IVHlU.iMl lO • 1 ! 


y. H. 

CSS: (Dept, of Electrics 


.t I'd Hidd! i I Id. / 


/'dESIgB'aWALVSIS^PROP^ time (COMPUTERS) 


\ ne use o T OF i. 
of hybrid vehi< 
t§chni. ^.ues mer ^ 
ii-ioiuht) hesi e'f 


opfiffitzatton [uethods as an effective desiyti^ toot 
vehicle propulsion systems is demonstrated^ 

! were used to select values for three desiun 

■at engine power rating and power split between Lhe rwo on-uoaro 


T|-,g principal conclusions are as ton.ou.fc. rirat^ ■ i, -'^'l -IP 
Htrateqv used to split the required power be luieen ine iwo dri-Doaru^erverdy 
sources can have a significant effect on life cycle coai ano^peiroieum 


POi !_ i ui is 


Second? the optimization program snouia oe coiiHi.ruci.eo 


cn 7U 



DOE/NASA/0084-1 
NASA CR-1 68292 


Optimization Methods Applied to 
Hybrid Vehicle Design 


John F. Donoghue and James H. Burghart 
Department of Electrical Engineering 
Cieveland State University 
Cieveland, Ohio 


June 1983 


Prepared for 

National Aeronautics and Space Administration 
Lewis Research Center 
Cieveland, Ohio 44135 
Under Grant NAG 3-84 


for 

U.S. DEPARTMENT OF ENERGY 
Conservation and Renewable Energy 
Office of Vehicle and Engine R&D 
Washington, D.C. 20545 

Under Interagency Agreement DE-AI01-77CS51044 






TABLE OF CONTENTS 


Page 

1.0 SUMMARY 1 

2.0 INTRODUCTION 4 

2.1 Research Objectives 4 

2.2 Previous Work 5 

2.3 Major Tasks 6 

2.4 Report Organization 7 

3.0 VEHICLE DESCRIPTION AND OPERATION 8 

3.1 Propulsion System 8 

3.2 Vehicle Operation and Ratio Selection 10 

3.3 Power Split Strategy 11 

4.0 VEHICLE DRIVING REQUIRETffiNTS 13 

5.0 ENERGY AND COST CALCUIJ^TIONS 16 

5.1 Yearly Electrical Eneigy Consumption 16 

5.2 Yearly Petroleum Consumption 17 

5.3 Life Cycle Cost Calculation 18 

6.0 VEHICLE SIMULATION 22 

6.1 General Coirputational Procedure 

and Trip Characteristics Design 22 

6.2 Simulation Time Increment Selection 25 

6.3 Power Overload 25 

6.4 Other Features of the Simulation 28 

7.0 OPTIMIZATION AS A DESIGN TECHNIQUE 30 

iii 



Page 


8.0 THE VEHICLE DESIGIV PROBLEM AS AN OPTIMIZATION PROBLEM 


9.0 SELECTION OF OP^flMIZATION fffiTHOD . 


37 


10.0 IMPLEMENTATION OF THE GRG PROGRAM 
FOR VEHICLE DESIGN PROBLEMS. . . 


il.O RESULTS AND DISCUSSIONS OF RESULTS 

11.1 Vehicle Power Source Sizing 

for Minimum Life Cycle Cost 

11*2 Vehicle Power Source Sizing 

for Changing Petroleum Costs 

11.3 Vehicle Power Source Sizing 

to Minimize Total Petroleum Consumption 

11.4 Power Split Optimization 

11.4.1 Parameterized AiResearch 

Power Split Function 

11.4.2 Polynomial Power Split Function , 

11.5 Summary of Results 


. 43 

. 43 

. 46 

. 48 

. 52 

. 54 

• 55 

. 57 


12.0 CONCLUSIONS 


59 


APPENDIX A - DICTIONARY OF VARIABLES 

APPENDIX B - CALCULATION OF VEHICLE MASS • . . . 

APPENDIX C - CALCULATION OF MANUFACTURING, 

MAINTENANCE AND REPAIR COSTS ’ . . . 

APPENDIX D - CALCULATION OF ACQUISITION COST 

APPENDIX E - VEHICLE SIMULATION 


XV 



Page 


E.l Conponent Models 91 

E.2 Calculation of Power Required at the Axle During 
the Acceleration and Cruise Phases of the STC 
and During Highway Driving 98 

E.3 Drive Train Power Calculations for Acceleration 

and Cruise Phases of STC and for Highway Driving 99 

E.4 Calculation for Coasting Phase of the STC 103 

E. v5 Calculation of Drive Train Power Regeneration 

During the Braking Phase of the STC 103 

APPENDIX F - BATTERY MODEL 106 

F. l Scaling the Battery Model 106 

F.2 Calculation of Battery Voltage and State of 
Charge During Discharge in Highway Driving 
and in the Acceleration and Cruise Phases 

of an STC 106 

F.3 Battery Charging at End of Trip 110 

F.4 Calculation of Battery Voltage and State of Charge 

During Regeneration in the Braking Phase of an STC. . . . 112 

APPENDIX G - CVT RATIO SELECTION PROCEDURE IIA 

APPENDIX H - AIRESEARCH POWER SPLIT STRATEGY 121 

APPENDIX I - COMPUTER PROGRAMS 122 

REFERENCES 


V 




1.0 SUf^.lARY 


The overall objective of this work was to demonstrate the use of optimi- 
zation methods as an effective design tool in the design of hybrid vehicle 
propulsion systems. 

Specifically, optimization techniques were used to select values for 
three design parameters (battery weight, heat engine power rating and a para- 
meter through which the required wheel power at each time instant was split 
between the two on-board energy sources) such that various measures of vehi- 
cle performance were optimized. The performance measures considered were ac- 
quisition cost, life cycle cost, petroleum consumption, or various weighted 
combinations of these. For specific values of the design parameters, vehicle 
performance was determined using a vehicle simulation. The optimization was 
carried out by coupling the simulation to a sophisticated, nonlinear program- 
ming code (GRG2) and running the two in combination. 

The optimization approach was successful in caning up with designs which 
were often significant improvements over hybrid designs already reported on 
in the literature. Furthermore, it was found that the strategy used to split 
the required wheel power between the two on-board energy sources can have a 
significant effect on life cycle cost and petroleum consumption. More impor- 
tantly, however, the work showed that it is quite feasible to couple a com- 
plex optimization program to a vehicle simulation program for the purpose of 
propulsion system design. It was also found that the optimization program 
should be constructed so that performance measures and/or design variables 
can be easily changed; in the conduct of the research it was often desirable 
to conpare, say, the minimum life cycle cost design to the minimum petroleum 



consumption design. 

Another area in which important results were obtained was in understand- 
ing the effect the vehicle simulation program design has on both the computer 
run time of the overall optimization program and on the ability of the opti- 
mization program to even arrive at a meaningful solution. First, it was 
found that the conputer run time could be significantly reduced by proper de- 
sign of the types of trips the vehicle takes in a one year period. Specifi- 
cally, if a trip consists of a combination of test cycles (accelerate- 
cruise-coast-brake cycles) and highway driving, then it is important that all 
the test cycles precede the highway driving, and that the same number of cy- 
cles are at the front end of each trip, no matter how long the trip. 

Second, it was found that care must be taken in designing the cost and 
constraint expressions which are being used in the optimization so that they 
are relatively smooth functions of the design variables. Otherwise, the 
function surfaces may be severly distorted and create many local minima which 
will prevent the program from even approaching the global minimion. This is 
particularly true with regard to how the battery replacement cost is handled. 

Finally, it was found that proper handling of constraints on battery 
weight and heat engine rating is particularly important for the success of an 
optimization study . These design variables cannot be allowed to be reduced 
to the point where the combined on-board power cannot meet the power demanded 
at the wheels. It was found that the way in which this "overload" constraint 
is handled strongly influences whether an optimization run gets "hung-up" at 
a non-optimal solution or finds the constrained, global optimum. 

The principal conclusion reached frcm this research is that optimization 
methods provide a practical tool for carrying out the design of the 
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propulsion system of a hybrid vehicle. Once the optimization program and the 
vehicle simulation are combined and operating, it is relatively easy to 
change objective functions and/or design variables so that designs based on 
different performance criteria (e.g. minimum life cycle cost design v.s. min- 
imum petroleum consumption design) or different driving conditions, component 
characteristics or propulsion systan configurations can be quickly carried 
out and compared. Use of optimization methods also force designers to make 
explicit their objective functions and constraints; thus designs, carried out 
by different researchers are more readily and meaningfully compared. 

This work has also lead to the conclusion that the way in which the de- 
manded power at any instant of time is split between the two on-board power 
sources has a significant impact on vehicle cost and petroleum consumption. 
Although some preliminary results in this area are presented in this report, 
it would appear that future hybrid vehicle work should look at this question 
of power split on a more systematic basis. 
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2.0 INTRODUCTION 


The energy crisis, dramatized by the oil embargo of 1973, has become one 
of the most serious long-term problems confronting the United States. In re- 
sponse to this crisis there has been a significant increase in the research 
and development of electric and hybrid vehicles in an attempt to reduce our 

dependence on foreign oil. The work discussed in this report is a part of 
that on-going effort. 

2.1 RESEARCH OBJECTIVES 

The primary objective of this work was to study the use of optimization 
methods for the design of propulsion systems for hybrid vehicles. The term 
hybrid vehicle as used in this report means a vehicle which contains a heat 
engine and gas tank as well as an electric motor and batteries. The configu- 
ration used here was a parallel hybrid with a continuously variable transmis- 
sion. Of particular interest in this work was determining the ease or diffi- 
culty with which a sophisticated, nonlinear programing algorithm could be 
coupled to a vehicle simulation program for purposes of optimally designing 
the vehicle propulsion systan. The design variables were heat engine power 
rating, the weight of the on-board batteries, and a parameter through which 
the power required at the wheels at each time instant was split between the 
two on-board power sources. The study considered optimizing various vehicle 
performance measures to determine how different performance measures influ- 
enced the final propulsion system design. The performance measures used were 
acquisition cost, life cycle cost, and petroleum consumption. 
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2.2 PREVIOUS WORK 

There have been no studies reported in the literature which attempt a 
vehicle propulsion system design using optimization methods. Two recent 
studies [1,2] have looked at different aspects of the design of various types 
of hybrid vehicles. In each case, however, a final design is arrived at us- 
ing a sound, but nevertheless ad hoc, selection process. Furthermore, the 
performance measure used to guide the selection is often not clearly stated, 
is modified or changed in the middle of a study, or is quite different fran 
study to study. These factors cc»nbine to make it difficult to compare the 
results of one study to those of another. In contrast to this, use of opti- 
mization methods force standardization by requiring a designer to explicitly 
state the performance measure to be used, and to maintain the sane measure in 
effect throughout the process of selecting final values for the design vari- 
ables. Furthermore, with an optimization approach the selection process it- 
self is standardized, and it also generally identifies the "best" values for 
the design variables so that the performance measure used is optimized (maxi- 
mized or minimized, whichever is appropriate). 

From the perspective of previous work on the problem of hybrid vehicle 
propulsion system design therefore, the work reported here represents the 
first time an attempt has been made to apply optimization methods to the de- 
sign problem in order to come up with what may be considered the best or "op- 
timal" design. 
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2.3 MJOR TASKS 


The conduct of the research fell naturally into a group of major tasks 
which are listed below: 

1. Developnent of a vehicle simulation to be used to predict 

petroleum and electrical energy consumption as the vehi- 
cle is driven over the specified driving pattern. These 
consumptions provide a basis for computing overall per- 
formance measures such as life cycle cost. 

2. Selection of a nonlinear programming algorithm (and as- 
sociated program) which will be used to perform the op- 
timization calculations. 

3. Coupling the vehicle simulation program to the optimiza- 
tion program to allow optimal propulsion system designs 
to be carried out. In carrying out an optimal design 
calculation, the optimization program must run the simu- 
lation repeatedly. Consequently, a good deal of the 
work involved in coupling the two programs is in tailor- 
ing the vehicle simulation so that it runs efficiently, 
and so that vehicle performance is a smooth function of 
the design variables. 

4. Translating various design goals into optimization prob- 
lems. For each design goal this would include specify- 
ing (i) the system performance measure to be optimized, 

(ii) the design variables, and (iii) the proper con- 
straints. 
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5. Generation and analysis of results using the combined 
optimization/vehicle simulation program. 

In the sections which follow, each of the major task areas outlined 
above will be considered in depth. 

2.4 REPORT ORGANIZATION 

The report is organized as follows. Sections 3 throv^h 6 present de- 
tailed information about the vehicle and its simulation. In particular, Sec- 
tions 3 and 4 describe the vehicle and its propulsion system as well as the 
vehicle's basic method of operation along with the driving requirements it 
must meet over a ten year period. Section 5 details the way in which the 
various vehicle performance measures (life cycle cost, etc.) are calculated, 
while Section 6 describes the vehicle simulation. 

The optimization design approach as it relates to hybrid vehicles is 
covered in Sections 7 through 9. Section 7 presents a general discussion of 
optimization methods which is useful background for subsequent sections. 
Section 8 deals with formulating hybrid vehicle design problems as optimiza- 
tion problems, while Section 9 describes some of the features of the actual 
optimization algorithm used in the study, as well as the reasons for its se- 
lection. 

Section 10 discusses how the vehicle simulation program and optimization 
program were combined, and Section 11 presents and discusses particular vehi- 
cle design problems and the results obtained using the optimization approach; 
conclusions are presented in Section 12. 
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3.0 VEHICLE DESCRIPTION AND OPERATION 

The vehicle selected to be studied in this research was essentially the 
same as the vehicle studied in [1]; it is a five-passenger sedan with a test 
mass of 2,049 kg. The details of the vehicle weight calculations are given 
in Appendix B. 

3.1 PROPULSION SYSTEM 

The propulsion system is a parallel hybrid configuration as shown in 
Figure 3-1. The traction motor is a brushless, DC, permanent— magnet motor op- 
erating from a battery voltage of 168 V. The motor is rated at 20 kw and has 
a top speed of 14,000 rpm. The ouput of the motor is connected to the trans- 
mission input through a gear which provides a 3.5:1 reduction of the motor 
speed. 

The heat engine is assumed to be a four cylinder engine. Since in the 
optimization work the kilowatt rating of the engine is a design parameter 
which is varied, no specific rating for the engine is given. However, in the 
optimization studies that were conducted , the heat engine would normally 
start out at a rating of 65 kw. The output of the heat engine is directly 
connected to the input of the transmission (no speed reduction). 

The transmission is a continuously variable transmission (CVT) which 
provides a range of ratios (input/output) from a 0.3:1 overdrive to a 3.33:1 
speed reduction. The output of the CVT feeds the differential through a 12:1 
speed reducer. The differential is assumed to have a ratio of 1:1. Thus, if 
the traction motor was running at its maximum speed of 14,000 rpm, and the 
CVT ratio was 0.3:1, the differential input would have a speed of 1,111 rpm. 
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The batteries are assumed to be lead-acid batteries with a system volt- 
age of 168 V. In the optimization runs in which the battery weight was var- 
ied it was assumed that the system voltage remained the same and that a vari- 
ation in battery weight meant the plate area of the battery cells was in- 
creased or decreased. A detailed description of the battery model used in 
the vehicle simulation is given in Appendix F. 

Finally, it was assumed that a clutch existed between the heat engine 
and the CVT input to allow the heat engine to be decoupled from the system 
when it was not supplying power. 

Section 6 and Appendix E (Vehicle Simulation) provide a detailed de- 
scription of the simulation and the models used for all of the components in 
the propulsion system. Reference to that section and that appendix can give 
the reader a more thorough understanding of the components, their interrela- 
tionships, and how typical power consumption calculations were made. 

3.2 VEHICLE OPERATION AND RATIO SELECTION 

In order to avoid engine stalling, it was assumed that when the vehicle 
was starting from rest all power was supplied by the traction motor until the 
vehicle speed reached 4.82 km/hr (3.0 mph). 

Whenever the vehicle was to brake to slow down, it was assumed that re- 
generative braking was used to the extent possible. If the maximum power 
rating of a ccxnponent limited the fraction of the avialable power which could 
be regenerated , then it was assumed that normal braking was used to dissipate 
the excess energy. 

In the early stages of the work, the CVT ratio was selected such that as 
the vehicle speed changed, the traction motor always ran at its most 
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efficient operating speed. The disadvantage of this approach in the hybrid 
propulsion system studied is that the heat engine speed is also determined by 
the CVT ratio, and a ratio that operates the traction motor most efficiently 
often operates the heat engine very inefficiently. In fact, it was found 
that the traction motor efficiency is always in the 90% or above range inde- 
pendent of the ratio, whereas the heat engine efficiency can often be sub- 
stantially improved by proper choice of the ratio. For this reason the CVT 
ratio was always selected to maximize the heat engine efficiency (except when 
the heat engine was not delivering any power). We note that a better ap- 
proach might have been to base the ratio selection on the amount of power 
each source is delivering, as well as on the relative efficiency improvanent 
possible from each source. A detailed description of the ratio selection 
procedure, along with the engine and traction motor data used, is given in 
Appendix G. 

3.3 POWER SPLIT STRATEGY 

In order to operate a hybrid propulsion system it is necessary to spec- 
ify a strategy for splitting the required wheel power between the two power 
sources. The power split strategy must operate continuously, deciding at 
each instant of time how the total power is to be split. A realistic strat- 
egy generally is a function of vehicle velocity and acceleration, and battery 
state of charge. 

An important part of this research was to use optimization methods to 
identify "optimal" strategies, or to select optimal values for parameters in 
a given strategy based on common sense. An example of a connmon-sense strat- 
egy which has been reported in the literature is the strategy proposed 
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by AiResearch [1], The AiResearch (AR) strategy is often used in the work 
discussed in the following sections; for this reason it is described in de- 
tail in Appendix H. 
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4.0 VEHICLE DRIVING REQUIREMENTS 


In order to ccmpute petroleum or electrical energy consumption (impor- 
tant factors in th^nselves , but also imporant factors in computing life cycle 
cost) it is necessary to specify how the vehicle is to be driven over its ten 
year life. Table 4-1 lists the lengths of eight different trip types for the 
vehicle, as well as the number of times the vehicle takes each trip type in 
one year. 

Some portion of every trip is composed of Special Test Cycle (STC) driv- 
ing. The STC is shown in Figure 4-1 and is a modified SAE J227A, D cycle; it 
contains a 14 second constant power acceleration period, followed by a 50 
second, 72 km/hr (45 mph) cruise period, a 10 second coast period, and a 9 
second braking period. Referring to Table 4-1, trips 1, 2 and 3 (trips less 
than 80 km in length) are made up of just a sequence of STC's, while trips 4 
through 8 are made up of a combination of STC's and highway driving. Each of 
these latter trips consists of 28 STC's followed by enough highway driving to 
reach the total trip length. Highway driving is constant speed driving at 90 
km/hr (56 mph) . 

The structuring of the longer trips to consist of a fixed number of 
STC's followed by varying lengths of highway driving had a big impact on re- 
ducing the computer run time of the overall vehicle simulation/optimization 
program. This is discussed in detail in Section 6 which deals with the vehi- 
cle simulation. 
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Table 4-1 



Yearly Distribution 

of Trips 


Trip Length 

Nunter of Trips 

Km 

Mi 

Per Year 

10 

6.2 

130 

30 

18.6 

85 

50 

31.1 

57 

80 

49.7 

54 

130 

80.8 

12 

160 

99.4 

7 

500 

311.0 

3 

800 

497.0 

1 
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5.0 ENERGY AND COST CALCULATIONS 


This section details the methods used to calculate the yearly electrical 
energy consumption, yearly petroleum consumption and life cycle cost. 

5.1 YEARLY ELECTRICAL ENERGY CONSUMPTION 

The yearly electrical energy consumption calculation is based on first 
determining the total ampere-hours out of the battery for each of the trip 
types described in the previous section. The battery charging characteristic 
is then used to determine the energy required to replace those amp-hours (in- 
cluding a charging efficiency). Each of these results is multiplied by the 
number of such trip? each year for each trip type, and the results then added 
to produce the total yearly electrical energy consumption. The details of 
the calculation are given below. 

As explained in Section 6 on the vehicle simulation, it is necessary 
when calculating the electrical energy and petroleum consumptions for a par- 
ticular trip to calculate the powers delivered by the battery and the heat 
engine at each time instant throughout the entire trip; it is not reliable, 
in general, to do it for just a few STC's and a segment of highway driving 
and then extrapolate the results for the entire trip. In the simulation, a 
small time step size is used and, for the electrical enej^gy calculation, at 
each time increment in the trip the battery model is used to determine the 
battery voltage, current and state of charge for the power that must be de- 
livered to the motor in that time increment (the calculation takes into ac- 
count the chopper/inverter efficiency). In this way a total net amp-hours 
out of the battery for each trip type is determined: the total net amp-hours 
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out includes the effect of regenerating capacity back into the battery during 
the braking phase of each STC. This total net capacity removed, together 
with an assumed charging current of 25 amps, is used to determine an average 
charging voltage from the battery charging characteristic (Appendix F). The 
average charging voltage is then multiplied by the total net capacity removed 
to determine the total energy that must be restored to the battery to bring 
it to a fully charged condition. Finally, to determine the total wall-pli:^ 
energy, the above total energy is divided by a charger efficiency factor. 
This calculation is done for each trip type. The total, yearly electrical 
energy consumed is finally calculated by multiplying the energy consumed for 
each trip type by the number of such trips per year, and then adding these 
results . 

If we let AHj^, i=l,..., 8, be the net amp-hours out of the battery for 
trip type i, be average charging voltage to fully charge the battery after 
delivering AHj^ amp-hours, EFF be the charger efficiency, be the number of 
trips of type i made in one year, then the total, yearly electrical energy 
consumed, Bj-Qp, in megajoules, is given as 

8 

Etot = (3600/10^) I (NiAHpVi/EFF) (5-1) 
i=l 


5.2 YEAPILY PETROLEUM OONSUMPTIONS 

As noted above, when calculating electrical energy and petroleiam con- 
sumptions it is necessary to calculate the powers delivered by the battery 
and heat engine at each time instant throughout the entire trip. In the sim- 
ulation, the power delivered by the heat engine, the engine speed and the 
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torque at each time instant are determined. The heat engine map (Appendix E) 
is then used to determine the fuel flow rate into the heat engine; multiply- 
ing this flow rate by the simulation time increment gives the fuel consumed 
over that time increment. 

The procedure used to calculate the total, yearly petroleum consumption 
parallels the one used for electrical energy consumption: the total petro- 
leum consumed for each trip type is determined by adding the fuel consumed 
over each time increment in the trip, these are multiplied by the number of 
trips of each trip type made per year, and the results then summed to produce 
the total petroleum consumption per year. 

5.3 LIFE CYCLE COST CALCUIJ^TION 

The life cycle cost is the sun of the vehicle acquisition cost, the pe- 
troleum, electrical energy, maintenance and repair, and battery replacement 
costs over a ten year period, minus the salvage value of the battery, the ve- 
hicle body and the power train. 

It is assumed that the costs of petroleum and electricity increase in a 
linear fashion from year one to year ten. The arrays containing this infor- 
mation are GASCST( I ) and ravCST(I ) ; the particular values used are given in 
Appendix A. 

Life cycle cost is computed in 1976 dollars, and it is assumed that 
there is no inflation. The cost for a future year is discounted to a present 
value using a discount factor (DISCNT in Appendix A). 

The salvage value of the vehicle body and power train is taken to be 10% 
of the marked-up manufacturing cost, i.e. 10% of the difference between the 
acquisition cost and the list price of the battery (0.1*(ACQCST-BTLIST), see 
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Appendix A for definitions of variables). 

The years in which a battery replacement is required depend upon the 
battery weight and the way in which the battery and heat engine are used. In 
calculating battery life it was assumed that there were 840 deep discharges 
in the lifetime of a battery [1], where a deep discharge is defined to be a 
discharge to 20% state of charge or below. The calculation procedure was to 
determine how far the battery was discharged when the vehicle made each of 
the trips defined in Section 4 (this was determined during the simulation of 
the vehicle to determine petroleum and electrical energy consumptions). If 
the battery state of charge dropped to only , say , 40% at the end of a parti- 
cular short trip, then that trip represented 75% of a deep discharge. In 
this way the total number of deep discharges of the battery per year was de- 
termined and used to calculate the number of years in the battery life. If 
the battery state of charge fell below 20% during the part of a trip in which 
STC driving was being performed, regeneration back into the battery and sub- 
sequent use of the battery until its state of charge again fell below 20% was 
allowed, but only one deep discharge was attributed to that trip. The bat- 
tery replacement cost had a mark-up factor which was twice that used in cal- 
culating the initial cost of the battery. The expression used to determine 
the replacement cost is 

BTCOST = BATCST*(1.0 + 2.0*BTfQOJP) (5-2) 

where 

BATCST = BATCER+WGBATT (5-3) 

and where (see Appendix A) WGBATT, BATCER, BATCST, BTMKUP and BTCOST are bat- 
tery weight, cost estimating ratio (CER), manufacturing cost, mark-up 
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factor and replacement cost, respectively. The way in which the battery re- 
placenent cost was used in calculating life cycle cost was influenced by the 
method used to handle battery salvage value. In particular, battery salvage 
value was taken to be the cost of the last battery replacenent multiplied by 
the ratio of the number of years left in the battery (at the end of the ten 
year vehicle lifetime) to the number of years in the battery life. Further- 
more, instead of adding the total battery replacement cost to the year in 
which the replacement took place, each replacement cost was evenly distri- 
buted over the ten years in the lifetime of the vehicle. Although this is 
not quite representative of what would actually take place, it was necessary 
in order to produce a life cycle cost expression for the optimization studies 
which was a relatively smooth function of battery replacement cost. Further- 
more, if it is assumed that the vehicle owner meets the battery replacement 
cost by borrowing money which is paid back over several years, then the above 
approach begins to more closely approximate reality. Based on the above dis- 
cussion, the battery replacement cost assigned to each year in the ten year 
lifetime of the vehicle is given by 


BATRYR = {[10.0/BATLFE)-1.0]*BTC0ST}/10.0 (5-4) 

where 

BATLFE = DISLFE/DEEP (5_5) 

where BTCDST is as calculated above, and where BATRYR, BATUE, DISLFE and 
DEEP are battery cost assigned to each year, battery lifetime, number of deep 
discharges in the battery lifetime, and the number of deep discharges per 
year, respectively. 
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Based on the discussion above, in Sections 5.1 and 5.2 and in Appendices 
C and D, the total yearly cost for operating the vehicle in year i is (see 
Appendix A for complete definitions of the variables) 

VHCOSTi = BpOT*KWCSTj, + PrroT*GASCSTi + BATRYR 

+ NIAINT + REPAIR (5-6) 

where E>pox and P-fOT are the yearly electrical energy and petroleum consumed, 
respectively, KWCST^ and GASCSTj^ are the cost factors, and MAINT and REPAIR 
are the maintenance and repair costs, respectively. Letting the salvage 
value of the vehicle body and power train be (see above) 

SALVGE = 0.1*(ACQCST - BTLIRT) (5-7) 

and the discount factor for year i be 

(i-1) 

Di = (1.0 + DISCNT) (5-8) 

where DISCNT is the discount rate, the life cycle cost is given as 

10 

LCCOST = ACQCST + I (VHCOSTi/D^^) - (SALVGE/Dii) (5-9) 

i=l 
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6.0 VEHICLE SIMULATION 


In the optimization studies the vehicle simulation was used primarily to 
calculate the electrical energy and petrolevim consumed by the vehicle in a 
one year period. Appendix E presents a detailed description of the modeling 
of the vehicle (the vehicle is described in general terms in Section 3) and 
of the way in which the petroleum and the electrical energy consumption was 
determined at each instant of time in the simulation interval. Section 5 
describes the way in which these individual consumptions were used to calcu- 
late the yearly consumptions. The objective of this section is to discuss a 
number of features of the vehicle simulation effort (other than vehicle mod- 
eling) which were considered to have a significant impact on the simulation 
run time, on the outcome of optimization runs, or which represented particu- 
lar features of vehicle operation which should be made explicit. 

6.1 GENERAL COMPUTATIONAL PROCEDURE AND TRIP CHARACTERISTICS DESIGN 

The purpose of this subsection is to describe the overall computational 
procedure used to calculate the fuel and energy consumptions for each trip in 
Table 4-1, and then to discuss how the structuring of these trips (STC and 
highway driving) can have a significant impact on simulation run time. 

As discussed in Section 4, each trip to be simulated is either a se- 
quence of just STC's or a canbination of 28 STC's followed by varying lengths 
of highway driving. Since in this study the power split was allowed to vary 
with battery state of charge and/or vehicle acceleration, the power supplied 
by either on-board energy source could vary significantly from one 
part of a trip to another, particularly on those trips which depleted the 
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battery. In the case of a depleted battery, the electrical energy supplied 
over the final STC would be zero, whereas that supplied over the first STC 
would be non-zero. Because of this, it was necessary to simulate the vehicle 
over the entire trip length for all trips: the fuel and energy consumption 
could not be calculated for one STC and a segment of highway driving, and 
then multiplied by appropriate factors to arrive at trip totals. 

The calcuation procedure used is as follows. The longest all STC trip 
is 35 consecutive STC's. The simulation was therefore run for 35 consecutive 
STC's, and the "state" of the vehicle at the end of each STC was saved. By 
"srate" is meant all the vehicle quantities necessary to resume the calcula- 
tion from that point forward in time (petrol eitn consigned, battery amp-hours 
out and regenerated, state of charge, voltage, instantaneous and average cur- 
rents, and whether the battery has discharged by the end of that STC). Know- 
ing these quantities at the end of each STC, the fuel and energy consumption 
for trips consisting of just STC's (trips 1, 2, 3 in Table 4-1) can be calcu- 
lated fron simulating just the longest all-STC trip (trip 3). 

To calculate the consumptions for trips involving highway driving, the 
simulation was initialized by the information saved at the end of the 28th 
STC, and highway driving was begun (in this way the 28 STC's preceeding high- 
way driving need not be resimulated). At each time step in highway driving a 
test was made to determine if the shortest STC and highway trip had been com- 
pleted; if so, the consumptions for that trip were computed and highway 
driving resumed until the next trip in Table 4-1 was completed. In this way, 
consiomptions for all trips consisting of STC's and highway driving can be 
calculated by simulating just the longest trip (trip 8, Table 4-1). 


23 



The above calculation procedure where fuel and energy consumptions for 
all highway trips can be calculated fron the longest trip depends on all 
highway trips having all of their STC's in the front of the trip, and all 
having the same_number of STC's up front. If either of these conditions is 
changed, it would become necessary to simulate each trip separately. This 
could represent a significant increase in simulation run time. For example, 
the time step size in simulating STC driving is 1 second, and that used in 
simulating highway driving is 60 seconds (see Section 6.2). Since each STC 
is 83 seconds long, and there are 35 STC's in the longest all-STC trip, the 
total number of drive train power calculations required for trips 1, 2 and 3 
in Table 4.1 is (1 x 83 x 35) = 2,905. The longest highway trip is 759 km of 
highway driving, driven at 90 km/hr. or a highway driving time duration of 
8.44 hrs. A 60 second simulation time step would therefore require 506 ad- 
ditional power train calculations for all of the highway trips (trips 4 
through 8), yielding a total of 3,411 calculations. If each highway trip had 
to be simulated separately (including the 28 STC's since they could be split 
differently between the beginning and end of each trip), the 28 STC's for the 
5 highway trips would require a total of (5 x 28 x 83) = 2,324 computations, 
and the highway parts of trips 4 through 8 would require 26, 60. 80, 306, and 
506 computations, respectively. The total number of computations would 
therefore be (26 + 60 + 80 + 306 + 506 + 2,324 + 2,905) = 6,207. Cbmpared to 
the above total of 3,411, it is seen that the number of drive train power 
computations has almost doubled. Since these computations are by far the 
major canputational burden of the simulation, and since the computation time 
is large to begin with, it is Important to design the driving pattern so that 
computational efficiencies can be achieved. 
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6.2 SIMULATION TIME INCREMENT SELECTION 

Another important factor in determining the length of time it takes the 
simulation to run is the time step size used in STC and highway driving. 
Since in a typical optimization run the yearly fuel and energy consumptions 
are computed 30 to 50 times, the choice of the time step size has a signifi- 
cant effect on overall computer run time. 

Separate time step sizes were used for STC driving and highway driving: 
since highway driving is at a constant velocity, it was felt that a lai^er 
step size could be used without missing significant changes in the vehicle 
state (e.g. battery state of charge). 

The time step sizes were selected by trying various values and observing 
the effects on petroleum consumption, electrical energy consumption and life 
cycle cost. In particular, the time step sizes were increased until a mean- 
ingful change in any of the above three quantities was observed; the values 
selected were set just below the values that produced a change. The STC time 
step size used was one second, and the time step size for highway driving was 
60 seconds. 

6.3 POWER OVERLOAD 

One of the most frequent optimization problans studied was to choose 
battery weight and heat engine rating to minimize life cycle cost. Without 
constraints on this problem the answer would be to set both design variables 
to zero. Of course, the vehicle would then not be able to meet the yearly 
driving requirements. In this research, the condition in which both power 
sources were reduced in size to the point where they could not supply the 
power required at the wheels was called power overload. This condition could 
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occur at any time in the yearly driving pattern, and it could occur for vari- 
ous combinations of heat engine rating and battery weight. The overload con- 
dition was handled in the optimization work by incorporating it as a con- 
straint: combinations of battery weight and heat engine rating which pro- 
duced overload were considered infeasible. 

One method (which did not work) of handling this constraint was to set a 
flag if an overload condition occurred anywhere in the yearly driving pat- 
terns. This flag would indicate to the optimization scheme that the current 
choices for the design variables were infeasible, and new values would be 
tried. The problem with this method was that the constraint was a discontin- 
uous (step) function of the design variables, and if it was violated the op- 
timization method (GRG2) would attempt to find exactly where, in the current 
search direction (see Section 7), the constraint was just satisfied. This 
attempt (using Newton's Method for finding the roots of an equation) would 
fail and the optimization scheme could not determine a useful, new search di- 
rection to locate the optiraimum. 

In an attempt to visualize the shape of the constraint in the battery 
weight/heat engine rating plane, a series of simulation runs were made in the 
region of the constraint boundary (varying battery weight and heat engine 
rating for each run, and using the AiResearch power split (see Section 3.3 
and Appendix H)). The result is shown in Figure 6-1 where it is seen that 
the constraint boundary is highly irregular. Even if the optimization method 
could find the boundary (for a particular search direction), once it began to 
move along the boundary it would get hung up at one of the many local minima. 
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Heat Engine Rating (Kw) 



The method finally adopted for handling the constraint was to approxi- 
mate it with two smooth curves as shown in Figure 6-1. For these smooth con- 
straint functions the optimization schone can easily find the location of the 
constraint boundary in a given search direction, and it can also move easily 
along the constraint boundary in search of the optimum. If we let — bat- 
tery weight (kg), and X 2 = heat engine rating (kw) , then the constraint equa- 
tions are given as 

PfjE > 37.0 - 0.02564 Wg (Wg < 273) (6-1) 

^HE ^ 77.23 - 0.2268 Wg + 1. 982x1 0”*^ (Wg > 273) (6-2) 

B 

The justification for using this approximate constraint boundary is that in 
any optimization study which chooses battery weight and heat engine rating to 
minimize life cycle cost, the objective is not necessarily to arrive at a de- 
sign which just borders on overload. Having the small "cushion" above over- 
load that use of the approximate constraint produces could be beneficial. 


6.4 OTHER FEATURES OF THE SIMULATION 

Several features of the simulation which should be mentioned are listed 
below. 

First, since a clutch existed between the heat engine and the CVT, the 
heat engine could be decoupled from the drive train when it was not being 
used. Hence it was assumed that there was zero fuel consumption by the heat 
engine during idling. 
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A second feature was that no matter how the power split strategy split 
the CVT input power, the heat engine was not allowed to deliver any power if 
the vehicle velocity was below 3 mph. This came into play as the vehicle 

started from rest at the beginning of a STC. 

A third feature (which is detailed in Appendix E) was that if either the 
battery or the heat engine could not deliver the power called for by the 
power split, the unmet power was assigned to the other power source and both 
the electrical and heat engine calculations were repeated. 
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7.0 OPTIMIZATION AS A DESIGN TECHNIQUE 


It is now appropriate to consider optimization as a technique which can 
be used with the vehicle simulation previously described to obtain desired 
design results. It will beccxne apparent in the discussion below that the 
role of optimization is that of systematically changing design parameters for 
use with the vehicle simulation until design objectives are met. Optimiza- 
tion techniques are based on a particular methodology and formalism which 
must be adapted to the design problon at hand. The purpose of this section 
is to define the terminology of optimization, briefly summarize the optimiza- 
tion problem mathematically, and discuss the general nature of solutions to 
such probloTis. 

Each optimization problon has an objective function (or criterion func- 
tion), f(x), which is an algebraic function of the n components of the vector 
X, and which is to be minimized (or maximized). The objective function in a 
design problem provides a measure of the quality of the design such that its 
minimum value occurs for a value of x which represents the best design. 

The components of x may be subject to upper and lower bounds. 


Xj. < Xj. 

min 


< Xv 


max 


(7-1) 


where, for a given variable Xj., one, neither or both of these bounds may be 
present. Variables often are given bounds for reasons of physical realiz- 
^•bility; e.g. , both battery weight and heat engine power rating must be posi- 
tive, and the percentage of power delivered from one source must be between 0 
and 100. 
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In addition to bounds which apply to the variables themselves, we must 
consider constraints , functions of these variables which also must be main- 
tained within limits. An optimization problem might contain m constraints 

g. < g.(x) < g. (j = l,2,..,m) (7-2) 

‘^min *^max 

where one or both of the limiting values on a particular constraint gj(x) may 
apply. The designer uses such constraint functions to calculate design char- 
acteristics which must be kept within upper or lower limits. A vehicle de- 
signer, for example, might want to place a lower limit on accelerating abil- 
ity, or an upper limit on acquisition cost. 

Since the objective functions and constraint functions applicable to hy- 
brid vehicle design are generally nonlinear functions of x, we have a nonlin- 
ear optimization problem. Converting the scalar variables in equations 7-1 
and 7-2 to vectors allows us to state the nonlinear optimization problem in a 
precise mathematical form: 

min f(x) 

X 

subject to 

%in < X < ^ax (7-3) 

and 

6min ^ S(x) < gj^jax (7-4) 

The notation min f(x) means that f(x) is to be minimized with respect to x. 

X 

Equations 7-3 and 7-4 are the vector versions of equations 7-1 and 7-2, 
respectively, and include all desired bounds and constraint relationships. 

Nonlinear optimization problems such as the one stated above have been 
well known for some time, and many techniques have been developed for finding 
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their solutions. Since f(x) and g(x) are assumed to be nonlinear, the 
techniques for solution of the problem are called nonlinear programming algo- 
rithms. Many of these have been developed over the past 25 years to provide 
a means of obtaining digital computer solutions to these problems [6, 7, 8, 9] 
Some of these methods are discussed in Section 9. 

Techniques for solution generally begin with an initial estimate, Xq, of 
the solution vector, which is feasible in that all bounds and constraints are 
satisfied, and proceed in an attempt to satisfy all of the necessary condi- 
tions for minimum f(x) while maintaining x feasible. The process terminates 
when these conditions are met, or when no further progress is possible. 

Each nonlinear programming method has its own set of assumptions regard- 
ing the mathematical properties of f(x) and g(x). Although there are small 
differences from method to method, it is usually necessary for both f(x) and 
the elanents of g(x) to be continuous, differentiable functions of x, with 
continuous first partial derivatives with respect to the components of x. 
These restrictions have important implications on the specifications of the 
objective and constraint functions for a design problem, as well as the 
development of simulation programs for optimization. 
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8.0 THE VEHICLE DESIGN PROBLEM AS AN OPTIMIZATION PROBLEM 


Having considered a general description of an optimization problem in 
Section 7, it is now possible to consider application of optimization methods 
to the design of the hybrid vehicle described earlier in this report. This 
involves formulating the vehicle design problem in terms of vehicle design 
parameters which will become the variables of the optimization problem, de- 
fining any bounds or constraint functions associated with these variables, 
and selecting an objective function which represents the designer's primary 
goal in a mathematical fashion. The key to successful formulation of a vehi- 
cle design problem as an optimization problem is in specifying all of the 
needed ingredients for the optimization problem (variables, objective func- 
tion, bounds and constraints) in a way which meaningfully includes all impor- 
tant aspects of the design problem. 

The primary design variables for the hybrid vehicle design problem dis- 
cussed here are the two power-sizing parameters, heat engine power rating (in 
kilowatts) and battery weight (in kilograms). Additional design parameters 
are introduced later (see Section 11) as part of a decision rule which is 
used on an instantaneous basis to decide how much power is to be used fran 
each of the two on-board power sources. 

Lower bounds of at least zero must be specified for both the heat engine 
power and the battery weight. This is because the optimization method will 
consider all possible values of the design parameters, even those which are 
physically meaningless, unless such bounds are imposed. From a designer's 
point of view, a hybrid vehicle which has a very low heat engine power rating 
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or negligible battery weight is a single— power— source vehicle and not a hy- 
brid vehicle; thus lower bounds were imposed to insure that the solution of 
the design problem by optimization produced a vehicle with both sources. The 
values of the lower bounds xxsed in this study were 15 kw for the heat engine 
power rating and 60 1^ for the battery weight. 

Two primary objectives were considered in this study: life cycle cost 
($) and life cycle petroleum consumption (liters). Both of these can be cal- 
culated using the calculation procedure discussed in Section 5, and the vehi- 
cle simulation discussed in Section 6. 

The principal constraint function used in this study is the power over- 
load constraint discussed in Section 6. It was necessary to impose this con- 
straint in order to ensure that the designs resulting from the use of optimi- 
zation represent vehicles which can produce sufficient power to successfully 
complete all the specified driving regimes. Other constraint functions can 
also be used for particular design problems, using as a constraint a function 
which could be an objective function for another design; for example, the de- 
sign might call for the hybrid vehicle which has the minimum life cycle cost 
subject to the constraint that its total petroleum consumption be less than 
some prespecified amount. 

There are many vehicle design problems which may be stated and solved as 
optimization problems. Listed below are several sample optimization problans 
for hybrid vehicle design. These same problems were solved as part of this 
study and are discussed later in this report where results are presented. In 
the statement of these problems, Clc is the life cycle cost (in $), Vp is the 
total petroleum volume (in liters) used over a life cycle, is the heat 
engine power rating (in kilowatts), Wb is the battery weight (in kilograms) 



and Pq is the power overload function (> 0 when no overload occurs). 


Choose PgE and Wg to minimize 


subject to 

^ Xln 


Po > 0 

Sample Design Optimization Problem No. 2 

Choose Pgg and Wg to minimize 


subject to 

BhE > Xln 


"B > \i„ 


Po > 0 

Sample Design Optimization Problem No. 3 


Choose Phe and W 3 to minimize VpCP^g, Wg) 
subject to PflE > 


in 


Wg > Wg . 

min 


Pq > 0 


ClC ^ ^LC 


max 


These three problems are typical of many which were considered in this 
study. All three are essentially vehicle power-sizing problems, where it is 
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desired to find the proper values of the heat engine power rating and the 
battery weight. They vary according to their objective functions and, in 

problem No. 3, the imposition of an additional constraint besides the power 
overload constraint. 

Other design problems considered later in this report involve constrain- 
ing the total petroleiom volume in the minimization of C^q, examining the ef- 
fect of petroleum pricing on design results, and parameterizing the power 

split decision rule. These problems are discussed with the results in 
Section 11. 
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9.0 SELECTION OF OPTIfllZATION METHOD 


Several methods including variable metric and conjugate gradient methods 
[6, 7, 8, 9] are available for direct solution of nonlinear programing prob- 
lems such as the type stated in Section 7 and applicable to the hybrid vehi- 
cle design problem. Indirect methods such as those based on the Newton- 
Raphson method are also available, but this class of methods was rejected be- 
cause of the likelihood of divergence with a poor initial estimate of the 
solution. Most direct methods such as gradient and modified gradient method 
were originally developed for unconstrained optimization problems (no bounds 
or constraints) and must be modified to include constraints. 

One proven method for successfully incorporating bounds and constraints 
into a direct search method was developed by Abadie and Carpentier [8]. The 
method, called the Generalized Reduced Gradient Method (GRG), uses linearized 
constraints, defines new variables which are normal to some of the con- 
straints, and transforms the gradient to a new basis. 

Lasdon and Waren, et. al. [10] have developed an extensively refined 
code (computer program) for general purpose implementation of the GRG optimi- 
zation method. The code which is now called GRG2, has been successfully 
tried on test problems in competition with other nonlinear programming codes, 
and has consistently ranked high in the primary characteristics desired of 
such a code (accuracy, efficiency, reliability, and ease of use) (11, 12]. 
GRG2 easily incorporates bounds on the variables and a variety of con- 
straints. It provides great flexibility which includes a choice of search 
method, a variety of stopping criteria, provisions for scaling of variables 
and constraint functions, and several options regarding the amount of detail 
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provided in the printed output. 

The designer who wishes to use GRG2 for optimization results must pro- 
vide [13]: 

1. a subroutine which canputes the objective function and all of the 
constraint functions, given the optimization variables. 

2. a data file which includes an initial point for the search, upper 
and lower bounds on variables and constraints (as applicable), and 
names for all variables and functions. 

The designer may provide: 

1. a subroutine for computing derivatives for gradients (otherwise the 
prc^rajn uses finite difference approximations) 

2. a subroutine for reporting needed output information which is not 
directly related to the optimization calculations 

3. a file containing scale factors for the variables and constraints, 
and control parameters which specify how the detailed operation of 
the optimization search is to be carried out (e.g., search direc- 
tion method, stopping criteria, constraint tolerances) 

After a vehicle design problem has been properly formulated as an optimiza- 
tion problem, it is a fairly straight-forward procedure to write the subrou- 
tines and determine the data files indicated above. The most difficult part 
of this task is in writing the subroutine which computes the objective and 
constraint functions from the design parameters, since this normally involves 

the entire vehicle simulation as well as the computation of various vehicle 
costs . 
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10.0 IMPLEMENTATION OF THE GRG PROGRAM FOR VEHICLE DESIGN PROBLEMS 


It was desired to combine the GRG2 optimization program with the hybrid 
vehicle simulation program in a way which would provide flexibility regarding 
specification of objective function, constraint functions, variables and 
bounds, so that a variety of hybrid vehicle design problems could be studied 
merely by changing a data file in the program. 

This goal was accomplished by defining five functions of P^e 
which would suffice for the calculation of all possible objective and con- 
straint functions which were anticipated at the time the computer program was 
written. These functions are: 

1. life cycle cost ($) 

2. acquisition cost ($) 

3. total petroleum volume in life cycle (liters) 

4. linear coribination of above three functions 

5. power overload function 

The fourth function was included to provide a weighted sum of three objective 
functions for the designer who wishes to compromise and not select a single 
objective fran among the first three functions listed. 

As described in Section 6.3, the power overload function was developed 
by simulating the hybrid vehicle for various pairs of values of the param- 
eters (PhE'^^b) using a particular power split strategy developed by AiRe- 
search Manufacturing Co. [1,5]. Regions where the hybrid vehicle was not 

able to meet the power demanded by the driving regime were identified and the 
power overload boundary was seen to be an irregular surface (Fig. 6-1) 
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which was not amenable to treatment by nonlinear programtiing techniques; 
these techniques normally require a smooth, differentiable function for each 
constraint. In order to meet this requirement in at least an approximate 
way, the overload constraint boundary function, approximated 

by two polyncsnial functions as follows: 

PjjE > 37.0 - 0.02564 Wg (Wg < 273) (6-1) 

Pgg > 77.23 - 0.2268 Wg + 1.982x10-4 w2 (Wg > 273) (6-2) 

B 

These functions were used with GRG2 for determining solutions to optimization 
problons which included the power overload constraint. 

For examination of the power flow or power split strategy problem, where 
we are trying to find the parameters of a decision rule which specifies how 

much power comes from each on-board power source at each instant of time, the 

GRG/simulation program was provided with three different ways to determine 
the power split between the heat engine and the electric propulsion system: 

1. the AiResearch algorithm [1] 

2. a polynomial function 

3. a parameterized version of the AiResearch algorithm 

More detail on the particular functions involved will be provided in Section 

11 . 

Further flexibility in the program was provided by defining a data file 
which could be used to specify which type of vehicle design problon was to be 
solved with a particular optimization run: (1) vehicle sizing, (2) power 

split optimization, or a conbination of the two. 
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The overall organization of the vehicle design computer program involves 
an interaction between the GRG2 optimization program and the hybrid vehicle 
simulation program discussed in Section 6. The interconnection is illustrat- 
ed in Figure 10-1. The reader should note that the GRG2 program selects val- 
ues of the parameter vector, x, for simulation in the Hybrid Vehicle Simula- 
tion Prograim, and receives the calculated values of the five functions (which 
may be the objective function or constraint functions) from it. GRG2 then 
systematically adjusts x in such a way as to minimize the objective function 
while satisfying bounds and constraints. The REPORT subroutine in GRG2 is 
used to calculate and print output information about the hybrid vehicle which 
is not essential to the GRG2 prc^ram. 

The next section discusses particular problems which were attempted 
using this program, and the results which were obtained. 
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11.0 RESULTS AND DISCUSSIONS OF RESULTS 


Four different design problems were solved using the optimization ap- 
proach discussed previously. Each of these problens is discussed in a sepa- 
rate subsection where the problan is defined and results are presented for 
that particular problan. 

The two major problem areas considered were those of vehicle power 
source sizing and power split i>arameter optimization. The first three design 
problems deal with the sizing problan, including various design objectives 
and constraints, and the final design problem is concerned with optimizing 
the power split decision rule. 

11.1 VEHICLE POWER SOURCE SIZING FOR MINIMUM LIFE CYCLE POST 

The design objective here was to determine the proper power source siz- 
ing parameters, heat engine power (kw) and battery weight (kg), so that the 
life cycle cost of the hybrid vehicle is minimized and the final vehicle is 
able to meet all driving requirements as specified by the annual trip re- 
quiranents and reflected in the vehicle simulation. 

The optimization problan for this design problem is very similar to Sam- 


ple Problem No. 1 in Section 8 of this report: 

Min Clc(Phe»%) (H-1) 

PhE'^B 

subject to > 15 (11-2) 

Wb > 60 (11-3) 

Po^%E»%) > 0 (11-4) 


and the AiResearch power split function which is specified in Appendix H. 
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Several initial estimates for the solution were tried, and all led to 
the same result: 

Min CLQ(Pfjg,Wg) = Cj^ (35.5, 60.0) = $17,100 (11-5) 

This point is at the intersection of the overload constraint boundary and the 
lower bound for Wg. See Figure 11-1, which shows the optimization as it pro- 
ceeded frcm the design parameters reconinended by AiResearch [1], Pgg = 65 kw 
and Wg = 386 kg. A1 thought this point was reached in two iterations, it 
should be noted that each iteration included many evaluations of the objec- 
tive function for the various pairs of (Phe» %) values which were being con- 
sidered. During the first iteration, the search proceeded in the negative 
gradient direction until the overload constraint boundary was reached, using 
longer and longer steps until the boundary was encountered. During the 
second iteration, the search moved along the overload constraint boundary un- 
til the lower bound on Wg was reached; several intermediate points were used 
along the constraint before the iteration was completed. The overall result 
is a reduction in life cycle cost from $20,796 to $17,100. 

It is not surprising that the optimal design in this case occin:*red on 
the overload constraint boundary, since any point above this boundary would 
correspond to higher value of %e, Wg or both, leading to a higher life cycle 
cost, and any point below this boundary, although yielding a lower life cycle 
cost, would correspond to a vehicle not able to meet the specified driving 
requirements. It is also not surprising that the optimal solution occurs at 
the minimum value of Wg, since the nominal petroleum pricing used for this 
problem (See Section 5.3) makes it desirable to use as little electric power 
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as possible if minimum life cycle cost is the objective. 

A modified version of this problem was solved after including an addi- 
tional constraint, 

Vp < r Vp* 

to the optimization problem specified in equations 11-1 through 11-4. V^* is 

the petroleum volume, 9,814 liters, associated with the solution to the 
original optimization problem, and r (0 < r < 1) specifies a certain fraction 
of that amount. The purpose of this design problem was to investigate how 
the optimal power sizing parameters would change, and how the life cycle cost 
would increase in the case where less and less total petroleum volume is 
available for the life of the vehicle. This is a situation which could be 
quite realistic as world petroleum resources are depleted or when foreign 
supplies beccxne unavailable. 

The results for this modified problan are shown in Figure 11-2 for vari- 
ous values of r, from r = 1.0 to r = 0.60. The results show a steady in- 
crease of Clq as r is decreased, rising from $17,100 at r = 1.0 to $18,679 at 
r = 0.60. These results demonstrate the utility of a hybrid vehicle when 
petroleum supplies are strictly limited, and more reliance must be placed on 
the electrical propulsion system. 

11-2 VEHICLE POWER SOURCE SIZING FOR CHANGING PETROIjEUM OaS TS 

The design problem considered here is similar to that of the last sec- 
tion, where life cycle cost is to be minimized, but where the price of petro- 
leum is varied and its effect on the results is observed. The price of pe- 
troleum was varied by introducing a gas factor, y, which scales all of the 
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(Petroleum Usage Constrained) 



petroleum costs, y > 1.0 (y multiplies the array GASCST(I)). 

The optimization problem for this design is: 

Min Clq(Phe»%»t) (11-7) 

Phe*% 

subject to equations 11—2, 11—3 and 11-4. The AiResearch power split func- 
tion (see Appendix H) is again used and the value of y is held constant dur- 
ing the optimization, so that this problem represents a family of optimiza- 
tion problems, each corresponding to a fixed value of y* 

The several optimization problems all used the same starting point, 

= 65 kw and % = 368 kg. The results are shown in Figtire 11-3 for y = {1.00, 
1.25, 1.50, 1.75, 2.00, 2.50, 3.00, 5.00}. All results are on the boundary 

of the overload constraint function, with higher values of y yielding designs 
which rely more and more on the electrical propulsion system. The life cycle 
costs associated with these results are listed in Table 11-1. It is clear 
fixxn these results that the ncxninal petroleum pricing (y — 1.0) which was 
used for much of this study is too low to justify a hybrid vehicle on the 
basis of life cycle cost; the optimization would yield Wb = 0 if this were 
permitted. Higher gas factors, however, show designs which have a much 
greater reliance on the electrical propulsion system. 

11*3 VEHICLE POWER SOURCE SIZING TO MINIMIZE TOTAL PETROLElff! (X)NStJ<PTION 

The design problem considered here differs from ,the problems discussed 
in the last two sections in that the design objective was to minimize the to- 
tal volume of petroleum used over the expected life of the vehicle. The ra- 
tionale for this objective was that very low availability of petroleum could 
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Table 11-1 

Optimum Life Cycle Cost For Various Gas Factors 


Optimum 


Factor , y 

Clc ($) 


WB(kg) 

1.00 

17,100 

35.46 

60.00 

1.25 

18,112 

35.46 

60.00 

1.50 

19,123 

35.46 

60.00 

1.75 

19,956 

33.98 

117.86 

2.00 

21,046 

25.98 

309.96 

2.50 

22,334 

25.97 

310.10 

3.00 

23,495 

23.34 

336.74 

5.00 

27,926 

16.46 

428.40 
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make a hybrid vehicle whose petroleum consumption was minimal very attrac- 
tive, even at a scxnewhat higher price. Thus this problem, at least in its 
initial formulation, did not take costs into account explicitly. 

The design optimization problem here is: 

Min Vp(PHE,W 0 ) (11-8) 

Phe»'*^b 

subject to equations 11-2, 11-3 and 11-4, with the Airesearch power split 
function and y = 1 . 0 . 

The result for this problem involving no costs was a hybrid vehicle with 
Pjjg = 15.0 kw, and Wg = 457 kg. This design point is on the overload con- 
straint boundary and sets the heat engine power at its lower limit, yielding 
total petroleum usage of 5,209 liters (compared to 9,814 when life cycle cost 
is minimized) and a life cycle cost of $19,337 (compared to $17,100 when life 
cycle cost is minimized). The results illustrate the trade off that occurs 
when the objective is low petroleum usage rather than minimum cost: the pet- 

roleum consumption decreases and the cost increases. 

These results motivated another design optimization problem which was 
defined, starting with the problem above, by adding an additional constraint, 

Clc(PhE»%) < 18,000 (11-9) 

This in effect asks for petroleum consumption which is superior to (lower 
than) that of the minimum-life-cycle-cost petroleum consumption, but which 
limits the cost increase. The result, as expected, was intermediate to the 
previous results for minimum life cycle cost and minimum petroleum consump- 
tion, with petroleum consumption of 7,2061 liters and a life cycle cost of 
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$17,963. The design variables for this result were = 31.2 kw and W 3 = 
226 kg, a point on the overload constraint boundary. 

The addition of this constraint illustrates the great flexibility of the 
optimization approach to hybrid vehicle design. It is easy for a designer to 
add such a constraint and repeat the optimization procedure in order to gain 
the desired result. 

A summary of the results for the two problems discussed in this section 
as well as the minimum life cycle cost results are listed in Table 11-2. 

Table 11-2 

Summary of Optimization Results 

Function 


Minimized 

Constraint 

PfjECkw) 

WB(kg) 


Cjjc;($) 

Clc 

none 

35.5 

60 

9,814 

17,100 

Vp 

Cjjc < 18 , 000 

31.2 

226 

7,206 

17,963 

Vp 

none 

15.0 

457 

5,209 

19,337 


11.4 PO\yER SPLIT OPTIMIZATION 

This section discusses some particular problems which are in a class 
which may be described as power flow strategy problems. The general problan 
could include any number of on-board power sources, including storage de- 
vices, and is directed toward specifying the direction and amount of power 
flow for each power source at each instant of time. The result of solving 
the problem is called the power flow strategy because it is to be expressed 
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as a set of strategy algorithms, or decision rules, dependent on some of the 
vehicle dynamic variables and internal variables, which permit calculation of 
each power flow as these variables change with time. 

For the particular case of a hybrid vehicle considered here, with the 
two on-board power sources, the power flow strategy problem reduces to that 
of determining, at each instant of time, how much of the wheel power is sup- 
plied by the electric propulsion system, and how much by the heat engine. 
Battery charging during regeneration is handled as a separate consideration 
(see Appendix F). The power flow strategy then amounts to determining the 
power split between heat engine power and battery power, and so the problem 
was defined using a power split function, s(t), 0 < s(t) < 1, where the ex- 
trames correspond to all-petroleiim (s = 0) and all-electric (s = 1) power. 

It is important to recognize the infeasibility of finding s(t) as the 
solution to an optimization problem. This would be possible only if we were 
to define a dynamic optimization problem with a prescribed vehicle driving 
regime (speed vs. time) for the entire life of the vehicle. This is unreal- 
istic froti two points of view: (1) the time interval involved would be so 
long that the computations could not be performed efficiently, and (2) it is 
not possible to predict the exact, second by second driving velocity of a ve- 
hicle realistically for the entire life of the vehicle, and for a substan- 
tial number of vehicles. 

Thus the approach taken here was based on the parameterization of a 
strategy function or decision rule, with the power split strategy expressed 
in terms of variables likely to influence the ideal power split function, and 
using parameters whose values could be obtained using the same methodology 
and computer programs used to determine optimal vehicle sizing parameters. 
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This part of the project included two different approaches concerning 
the parameterization of the power split function, s(t). The first involved 
the parameterization of the AiResearch power split function, and the second 
involved the use of polynominal functions to represent s(t). 

11.4.1 PARAtiETERIZED AIRESEARCH POWER SPLIT FUNCTION 

The power split function, previously discussed and used in the study by 
AiResearch Corporation [1] (also see Appendix H) , contains several constants 
which could be treated as parameters. One of these, which was thoi^ht to be 
one of the most critical, is the value of the power split function during ac- 
celeration, Sa» which was given the value 0.3 in the AiResearch study. An 
optimization problem was formed with s^ as a parameter: 


Min C[^(Sa) (11-10) 

subject to = 65 kw (11-11) 

Wb = 386 kg (11-12) 

0 < Sa < 1.0 (11-13) 


and the remainder of the AiResearch power split algorithm. Equations 11-11 
and 11-12 specify vehicle power sizing parameters which correspond to those 
of the AiResearch hybrid vehicle design; since this vehicle is easily able to 
meet the power demanded by the driving regime, it was not necessary to in- 
clude the power overload function as a constraint for this problem. 

This optimization problot was solved repeatedly for several values of 
the gas factor, y, which was discussed in Section 11.2. Results, which are 
listed in Table 11-3, show that no electric power should be used during ac- 
celeration for the lowest values of the gas factor, and that the power split 
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rises to 0.264 with y = 5.0. The optimal life cycle costs rises correspond- 
ingly with the gas factor, as expected. These results show that the optimal 
distribution of power between the two sources is strongly dependent upon the 
petroleum price, and that a power split value of 0.3 could be justified only 
for a very large value of the gas factor. 


Table 11-3 

Optimal Bower Split During Acceleration 


Gas Factor , y 

Power Split, Sg^ 

Life Cycle Cost, Clq($) 

1.0 

0.000 

20,000 

1.5 

0.000 

21,886 

2.0 

0.114 

23,303 

2.5 

0.194 

24,669 

3.0 

0.231 

26,021 

4.0 

0.248 

28,702 

5.0 

0.264 

31,364 


11.4.2 POLYNOIIAL P0\YER SPLIT FUNCTION 

Another approach to expressing the power split function in parameterized 
form is to use a simple polyncxnial function of vehicle variables with para- 
meters to be determined through optimization. This leads to a problem of 
constraining s(t) to its specified range, 0 < s(t) < 1, since any polynomial 
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function of vehicle variables is likely to exceed this range during at least 
part of the driving regime. This problem cannot be treated directly through 
constraints on s(t) because that involves a dynamic optimization problem 
which is unrealistic in this application as discussed earlier. 

The power split function was constrained through the use of limiting. 
Let Sp(t) be the value of the power split polynomial function. Then 

( 0. Sp(t) < 0 

Sp(t), 0 < Sp(t) < 1 (11-14) 

1. Sp(t) > 1 

The polyncmial function chosen for investigation involved two key vehi- 
cle variables, battery state of charge, c(t) and vehicle acceleration, a(t). 
Several polynomial functions were tried starting with the linear polynomial 

Sp(t) = OQ + ai c(t) + a.2 a(t) (11-15) 

and later including higher-d^ree polynomial functions. 

The search procedures of GRG2 did not produce worthwhile results with 
this approach. Typically, the search started with an initial estimate of the 
a parameters and failed to improve the objective function in several itera- 
tions, finally stopping because of lack of progress at a point very close to 
the initial estimate. Our conclusion about this is that it was not due to 
any deficiency in GRG2, but was due to a fundamental problen with the ap- 
proach, in particular the limiting process defined by equation (11-14). The 
limiting apparently led to erroneous values of the gradient so that the 
search directions determined by GRG2 were not effective. It seems as if the 
changes used for the parameters had little or no effect because the value of 
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Sp(t) calculated fran the polynomial such as that in equation (11-15) was 
promptly and frequently overridden by the limiting action as the vehicle sim- 
ulation was carried out. Approaches like the one described in Section 
11.4.1, where a proven and acceptable power split function was parameterized, 
seem much more promising. The main reason for this is that such an approach 
allows incorporation of bounds on the parameters thenselves in a direct man- 
ner, rather than the indirect approach used in the limiting of the polynomial 
function. 

11.5 SUMMARY OF RESULTS 

This section has treated the application of an optimization technique 
based on nonlinear prograrmiing to the design of hybrid vehicles. Four dif- 
ferent design problems were solved using the digital computer program which 
combined the optimization algorithm with the hybrid vehicle simulation. 

The method was successfully applied to vehicle power source sizing prob- 
lans for a variety of design goals and specifications. It was found that all 
solutions were on the power overload constraint boundary such that, in each 
case, any reduction in the battery weight or heat engine power rating would 
yield a vehicle which would not be able to meet the power demands of the 
specified driving regime. 

Two methods were considered for parameterizing the power split function 
which specifies, at each instant of time, what fraction of the total power 
demand is supplied by the electric propulsion syston, and what fraction by 
the heat engine. The parameterization of the AiResearch power split function 
was found to be much more successful than using a polynomial power split 
function with limiting. 
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The use of optimization as a systematic approach to a variety of hybrid 
vehicle design problans was shown by these results to be very effective. 
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12.0 CONCLUSIONS 


This study has led to several conclusions regarding the effective use of 
an optimization technique for hybrid vehicle design. 

Efficient and effective use of an optimization technique for hybrid ve- 
hicle design depends heavily on a vehicle simulation program which is de- 
signed for this purpose. One of the main necessary features of such a simu- 
lation is that it must permit continuous variation of the power plant sizing 
parameters, battery weight and heat engine power rating. 

Another area denanding attention in such a simulation is the way the 
specified driving pattern, including standard test cycles and highway driv- 
ing, are incorporated into the simulation; this has impact both on ccxnputer 
run time, which must be minimized, and on the ability to accurately simulate 
long trips involving a declining state of charge for the batteries, since 
this state of charge affects the power split. One cannot get accurate re- 
sults by merely simulating one standard test cycle and extrapolating the re- 
sults . 

Attention must be paid to the simulation time step sizes used in the 
simulation. They must be small enough to provide sufficient accuracy in the 
calculation of vehicle variables and costs, while remaining as large as pos- 
sible for the purpose of minimizing conputer run time; this is particularly 
important when the simulation is used in optimization, since the optimization 
program must run the simulation program dozens of times during the completion 
of one design problem. 

Since the main purpose of the simulation program is to calculate func- 
tions which are used as the objective or constraint functions for the 
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optimization program, the simulation program must be designed to produce 
these functions as smooth functions of the design parameters. Otherwise the 
optimization program is likely to calculate search directions which will not 
yield a successful result because of functional discontinuities, or the opti- 
mization is likely to stall at a false local minimum. 

This research demonstrated that the handling of the power overload func- 
tion is critical to obtaining successful results for power-sizing design 
problens. Design results for this problen are very likely to lie on the power 
overload constraint boundary , since the region above it represents increased 
life cycle cost and petroleum consumption, while the region below it repre- 
sents a vehicle which will not be able to meet the power denand associated 
with the specified driving regime. The exact location of the power overload 
constraint boundary need not be determined; it is more important that this 
constraint be a smooth function so that the optimization program can work ef- 
ficiently. 

It is important to choose the optimization program and to carefully de- 
sign the interface between the simulation program and the optimization pro- 
gram in order to provide considerable flexibility to the designer in the 
choice of the design objective, constraints and variables for a variety of 
design problems. This project demonstrated that this can be done success- 
fully. 

Standard procedures for improving the efficiency of the optimization 
program, including scaling of the design variables and functions, were used 
and were judged to be effective. Gradient calculations by finite differences 
required the use of double-precision arithmetic. These are points which the 
designer must keep in mind when developing optimization program applications. 
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The two methods of parameterizing the power split algorithm, which spe- 
cifies how much electric power and heat engine power is used at any time, de- 
monstrated the feasibility of this parameterization. The polynomial power 
split function which was limited in value did not produce useful results be- 
cause the optimization program did not have a free hand in varying the power 
split function due to the limiting. The parameterization of the AiResearch 
algorithm was much more effective, and this general area of parameterizing an 
established algorithm is a promising area for future research. Optimization 
of the power split function, in particular, is a research area which should 
be pursued. 

Optimization has been shown to be a useful tool for systenatic design of 
hybrid vehicles. This research showed that it is possible to couple a so- 
phisticated, general purpose optimization program to a complex hybrid vehicle 
simulation and produce a successful hybrid vehicle design using a reasonable 
amount of computer time. 


61 



APPENDIX A - DICTIONARY OF VARIABLES 


i iN4,L 


PKLiGRfiH HVf'OI’ 


%U 


cmmttmmuuuummmmmummnmmmuumtmmum 
cut 

cut PROGRAM HWOP = HYBRID VEfUCLt: f'ERl-ORMANCE EVAl.UAHOH PROGRAM 

c 
c 

cumnmminmumumnumntmmmummtmummuum 
ctuntMU uttutiu 

cmuum DICTIONARY OF OARIABLEC UtUmU 

iiuttutut uuutut 

c%'muutuuunu%uuunuuuuumunntnunuuuuuutuunu 
c 
c 


C VARIABLE 
C 


DESCRIPTION 


C AACTL 

i: 

c 

C AnCy( 

C AACfl 
U 

C ACCEi. 

C 


ACTUAL. CAPACITY WITHDRAUN 
I'RiJM BATfERY 


A 2S AMP CHARGE RATE) 

AVERAGE ACCELERATION OVER ONE 
DRIVING CYCLE TIME [NCREMEMT 


C 

C ACQC3T ACQOISITION COST OF VEHICLE 

n 

L- 

C ACR ARRAY OF CAPACITY REMOVED (USED 
C IN CALCULAriNG IIP AVERAGt: 

C CHARGING VOLTAGE) 

C 

C Al.LLTA(I) ARRAY (IF ILECTRiC ENERGY, 

C COMPUILT) BY SOMHINC ITIL CUKRI 

C OUT Of THE .liATT, IN TAF VARTd 

C 
C 
C 

i: 

C 

c 
c 
c 
c 
c 


IRIPS (FOR PRINTING) 

ALLLTPd) ARRAY OF ELECTRIC ENERGY, 

COMPUfEi) liY SUMMING THE POOER 
LUn Of the; BATT, in the VARIOU: 
fRIPS (FOR PRINTING) 

AESPED(l) ARRAY OF ENGINE SPEEDS USED IN 
HE MAP INTERPULAIION 


UNITS 

SOURCE 

VALUE 

uuuunnuutuuuuuutu 

AMP -HOURS 

CALC 

utuu 

E VOLTS 

IR 

CSU 

HAVE 

METERS PER 
SECOND SRUARED 

CALC 

tutut 

DOLLARS 

CALC 

unm 

AMP -HOURS 

NASA 

HAVE 

PER KILOGRAM 



MEGAIOOLS 

r 

CALC 

tuutt 

MEGAJOLILS 

CALC 

utuu 

RAD/SE.C 

NASA 

HAVE 


62 



C AHCHG 

AHP-HRS INTO BATTERY FROM HE 

AMP -HOURS 

CALC 

C 

CHARGING 



If 

C AHEKW(I) 
C 

ARRAY CONTAIMING HE POWER 
RATINGS (USED IN FUEL RATE 

KILOWATTS 

NASA 

C 

C 

CALCULATIONS) 



u 

c 

C AHOHT 

BATTERY AMP -HOURS OUT 

AMP-HOURS 

CALC 

L 

C ALLSTC 
C 
C 
C 

C ALPHAl 

LENGTH OF LONGEST TRIP WHICH 
USES JUST SrC'S (HO HIGHWAY 
DRIVING) 

METERS 

NASA 

WEIGHTING PARAMETER FOR LIFE 
CYCLE COST IN OBJECTIVE FUNCTION 

- 

CSU 

C 

c 

C ALPHA2 

WEIGHTING PARAMETER FOR PETROLEUM 
COST IN OBJECTIVE FUHCTIUN 

- 

CSU 

C 

c 

C ALPHA3 

WEIGHTING PARAMETER FOR ACQUISITION 

- 

CSU 

C 

COST IN OBJECTIVE r-UNCTION 



C 

C APETT(I) 

ARRAY OF PETROLEUM USED IN 

LITERS 

CALC 

C 

THE VARIOUS TRIPS (FOR PRINTING) 



C 

C ARECEN 

o 

AMP-HOURS REGENERATED 

AMP -HOURS 

CALC 

t 

C ASHCER 

CONSTANT USED TO CALC. 

DOLLARS PER 

NASA 

C 

THE ASSEMBLY COST 

KILOGRAM 


C ASHCST 
C 

C ATORQ(I) 

VEHICLE ASSEMBLY COST 

DOLLARS 

CALC 

ARRAY OF TORQUES CORRESPONDING TO 

NM 

NASA 


i; . me utandako n kh eKGifc used ru 

C COHI’LITE THE SCALED ARRAY OF TORQUES 

C USED FOR A DIFFEREHT SIZED ENGINt 

C 


C BATCER 
C 

c 

C BATCST 
C 

CONSTANT USED TO CALCULATE 
COST OF BATTERIES 

DOLLARS PER 
KILOGRAM 

NASA 

MANUFACTURING COST OF 
OATTEIUE3 

DOLLARS 

CALC 

if 

C BATLFE 
C 

C BAT R PI. 

C 

C 

NO. OF YEARS IN BATTERY LIFt. 

YEARS 

CALC 

NO. Of .BATTERY REPLACEMENTS IN 
VEHICLE LIFE 


CALC 


%%%%%%% 

T AHEKW(I) 
i 37. 

Z 56. 

3 75. 

-1 U2. 

nuttt 

50,001.0 

VARIES 

VARIES 

VARIES 

nutn 

mt%n 

0.1984 

tittut 

HAVE 

1.87 

%unu 


63 



o o 


C 8ATSAL 
c 

SALVAGE VALUE OF BATTERY 

DOLLARS 

CALC 

C BE.NRGV 
i: 

c 

BATTERY ENERGY USED TOR THIS 
INCREHeNT in THE DRIVING CYCLE 

JOULES 

CALC 

C 

i: 

c 

GRADE ANGLE ARRAY FOR COMPLETE 
DRIVING CYCLE 

DEGREES 

CSU/NASA 

C BLTA 

C 

c 

GRADE ANGLE AT THIS POINT 
IN THE DRIVING i:YCLE 

DEGREES 

CALC 

C BBFC 

CONSTANT USED TO CALC MOST 

RAD/SEC 

CSU/NASA 

i; 

EFFICIENI ENGINE SPEED FOR A 



c 

c 

GIVEN PERCENT OUTPUT OF HE 

X RATED POWER 


c inxosT 

c 

c 

C BTLIST 
C 

COST OF BATTERY REPLACEMENT 
DURING VEHICLE LIFE 

DOLLARS 

CALC 

LIST COST OF BATTERIES 

DOLLARS 

CALC 

C BTMAIM 
C 

c 

BATTERY MAINTENANCE COST 

DOLLARS PER 
KILOMETER 

CALC 

C BIHKUP 

C 

c 

BATTERY MARKUP FACTOR 
(A CONS IAN I BETWEEN 0 AMD 1) 

- 

NASA 

C CYiPCTY 
c 

CELL CAPACITY 

AHP-HOURS 

CALC 

C CAPKG 

ij 

C 

C CDA 

CAPACITY OF A BATTERY 
CELL PER KILOGRAM 

AMP-HOURS/KG 

CALC 

DRAG COEFFICIENT TIMES 

METERS 

NASA 

\] 

r. 

VEHICLE FRONTAL AREA 

SQUARED 


C CHARGE(I) ARRAY WHICH SAVES AT THE END 
C 01- EACH GTC THE TOTAL AMP-HRS 

AMP -HOURS 

CALC 

C 

c 

INTO THE BATTERY FROM HE CHARGING 


C CKGCER 

CONSTANT USED TO CALCULATE 

DOLLARS PER 

NASA 

G 

c 

COS! OF CHARGER 

KILOGRAM 

C CHGCST 
c 

COST OF BATTERY CHARGER 

DOLLARS 

CALC 

C CHCEFF 
c 

BATTERY CHARGER EFFICIENCY 

- 

NASA 

C CHGfiCl 

CONSTANT USED IN CALCULATING 

KILOGRAMS 

NASA 

C 

C 

C CHGRKW 
C 

MASS OF CHARF.ER 

PER WATT 

POWER RAVING OF CHARGER 

KILOWATTS 

NASA 


uuu* 

tuttn 

ALL 0. 

tuiut 

4.2 

uunt 

tnuu 

tuuu 

0.3 

tnut* 

**u*u 

0.6 

S.89 

nuut 

0.85 

0.0 

14.5 


64 



C ClllCi 
C 

c cHicrff 

c 

c 

C CIUC3T 
r 

CONSTANT OSEIt IN CALCULATING 
liAS!) ur CilOl’PER (NVEHTER 

KILOGRAMS 
PEP KW 

NASA 

0 9i 

CONSTANT USED TO CALCULATE 
CHOPPER INVERTER COST 

DOLLARS PER 
KILOGRAM 

NASA 

14.00 

COST Of CHOPPER INVERTER 

D01.LARS 

CALC 

uutu 

C C.HIKW 

1^. 

c 

POWER RATING OF CHOPPER 
INVER i'ER 

KILOWATTS 

NASA 

40.0 

u 

C CmHAH 

c 

c 

CHOPPER INVERTER MAINTENANCE 
COST 

DOLLARS PER 
K HOME i'ER 

CALC 

tutut 

U 

t CHCiOSE; 

c 

c 

c 

C CUCf-R 

'0 

r» 

IF EQUAL TO i, ELECTRICITY COST 
COMPUTED FROM ELECYKCI). If EQUAL 
TO fl, ELEC (I ) replaces ELECYRII) 

- 

CSU 

0.0 

CONSTANT USED TO CALCULATE 
THE COST Uf (TIE CLUTCH 

DOLLARS PER 
K IT.OGRAH 

NASA 

0.0 

U 

C Cl.TCST 
r 

COST Of CLUTCH FOR VEHICLE 

DOLLARS 

CALC 

uut%* 

C CtiUNT 

C 

C 

f' 

VARTALBLE USED IN DETERMINING 
WlllCIi YEARS BAIT. REPLACEMENT 
COSTS ARE ADDED 

- 

CALC 

nnt%* 

c Ci-; 
c 

CONSTANT IN ROLLING 
RESISrANCE FORCE CALCULATION 

- 

CALC 


C CVTC;l 

c 

c 

C CVTCER 
\: 

CONSTANT USED TO CALCULATE 
CVT MASS 

KILOGRAMS 
I'ER KILOWATT 

NASA 

1.07 

CONSTANT USED TO 
CAl.CULiTTE CVT COST 

DOLLARS PER 
KILOGRAM 

NASA 

1.4P 

U 

c cvrcsT 
c 

C CVTKW 

r 

COST OF CVT 

DOLLARS 

CALC 

tutut 

POWER RATING OF CVT 

KILOWATTS 

NASA 

90.0 

L 

C CVTHAR 

c 

C 

c cxa) 

r* 

CVT MAINTENANCE COST 

DOLLARS PER 
KILOHETEI! 

CALC 

nnttt 

ARRAY CONTAINING POWER SPLIT 
PARAMETERS 

MIXED 

CALC 

ttttttt 

C DLEP 
C 

NO. Of DEEP DISCHARGES OF 
DAffERY DURING VEHICLE LIFE 

- 

CALC 

tutttt 

1* 

C DELR 

INCREMENT FOR DECREASING CVT 

- 

CSU 

O.i 
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c 

fiAflO TU PREUEHl' HOfOR FROd 




c 

STALLING 




\ f 

C Dtl.SVC 

FRACTION OF AN STC USED 


CALC 

ntnu 

i; 

c 

IN CURRENf TRIP 




c mu 
c 

SIMULATION TIME STEF SIZE 

SECONDS 

CALC 

untn 

C DLL'(H« 

SIMULATION TIME INCREMENT FOR 

SECONDS 

CSU 

60.0 

0 

r 

HIGHWAY DRIVING 




Lr 

C DLLV 

VEHICU. VELOCITY DECREASE PER 

KILOMETERS 

CALC 

nnut 

C 

fIME INCREMENT DURING 3HEAKIHG 

PER SECOND 



C 

C 

PHASE Of STC 




C: DIGRKW 

POWER RATING OF GEAR BETWEEN 

KILOWATTS 

NASA 

90.0 

i: 

OliTERENTIAl. AND CVT (GEE ALSO 




c 

PGRAIE) 




c; DHFKW 

WHEN RLIN2i=I!, THE AMOUNT HEKW 

KILOWATTS 

ARB 


c 

IS INCfiEMENTED FOR NEXT CCOMP CAlL 




c 

c 

WHEN fiUN21=i THE TOLERANCE ON HEKW 




C !)!.FCER 

CONSTANT USED TO CALCULATE 

DOLLARS PER 

NASA 

2AZ 

ij 

[: 

GUST or A DirrGRENflAL 

KILOGRAM 



C D1FC3T 
c 

COST OF DIFFERENTIAL 

DOLLARS 

CALC 

uutu 

C ill! FCi 

CONSTANT USED IN CALCULATING 

KILOGRAMS 

NASA 

0.66 

i: 

f 

MASS Of DIFFERENflAL 

PER KW 



C DIFFKU 

POWER RATING OF DIFFERENTIAL 

KILOWATTS 

NASA 

90. 0 

c 

r 

(ALSO 3Q: PPRATE) 




L? 

C DlSCHd 

) ARRAY CONTAINING THE NUMBER 

.. 

CALC 

unut 


Oi' BAflERY DISCHARGES PER YEAR 




C 

c 

FOR THIS PARTICULAR TRIP TYPE 




C D)SCMG( 

BATTERY DISCHARGE CURVE 

VOLTS/CELL 

NASA 

HAVE 

C i>UCj lLr 

c 

F) 




C DIGCNT 
r 

DOLLAR DISCOUNT RATE 

- 

NASA 

0.02 

C OMIEFF 

FACTOR USED DURING DISCHARGE 


CALC 

iuut 

C 

TO CONVERT BATTERY TERMINAL 




c 

c 

ENERGY TO INTERNAL ENERGY 




C HI3FAC 

THE AMOUNT OF BATTERY DISCHARGE 


CALC 

tutut 

C 

DURING ONE fRIP OF THIS 




c 

PARTICULAR TRIP TYPE 





66 



840. 


C 

C DISLI'E KO. or DEEP DISHCHARGE8 IN •• NASA 

C BATFERY LIFE 

C 


C DLTSTC 
C 

SIMULATION TIME STEP SIZE 
FUR 3FC DRIVING 

SECONDS 

CSU/NASA 

i.O 

L 

C DWGBAT 
C 
C 
C 

C LBTOT 
0 
C 

C EfJTOTL(I) 

C 

C 

i: 

r* 

WHEN RllN2i=i, OR 2 THE AMOUNT 
WUBAFF Vo INCREHENFED FUR NEXT 
GCUHP CALL 

KILOGRAMS 

ARB 


TOTAL BATTERY ENERGY USED SINCE 
START OF TRIP 

MEGAJOULES 

CALC 

ttittu 

ARRAY WHICH SAVES AT THE END 
OF EACH STC THE TOTAL BATTERY 
ENERGY USED SINCE THE START 
OF ALL SrC'S 

MEGAJOULES 

CALC 

ttttttt 

L 

C EFF(KM, 
i: HA FIG) 

r* 

CVT EFFICIENCY HAP 
(A NUMBER BETWEEN 0 AN!) 1) 

- 

NASA 

HAVE 

Ir 

C EfFBC 

c 

r* 

COMBINED BATTERY CHARGER EFT. 
AND BATTERY TURNAROUND EFF 

- 

NASA 

a. 75 

If 

C EFF NTH (PERCENT MOTOR EFFICIENCY HAP 

C RATED PWR, 

C SPEED) 

r 


NASA 

HAVE 

If 

C ELCC(l) 
c 
C 
C 

c 

c 

c 

C ElECSTd) 
i: 

ARRAY CONTAINING WALL PLUG 
ELECTRICAL ENERGY (COMPUTED 
ON AHP-HR BASIS) REUUIRED TO 
MEET ELECTRICAL ENERGY 
CONSUMPTION FOR THIS PARTICULAR 
TRIP IVPE FOR ONE YEAR 

MEGAJOULES 

CALC 

tutu 

ARRAY CONTAINING ELECT. COST 
FOR YEAR I 

DOLLARS 

CA!.C 

tttttii 

\r 

C ELCCYR(l) 
t: 

c 

c 

c 

c 

ARRAY CONTAINING WALL PLUG 
ELECTRICAL ENERGY (COMPUIED 
ON POWER BASIS) REQUIRED TO 
MEEf LlECTiUCAL ENERGY 
CONSUMPTION FOR THIS PARTICULAR 
TRIP TYPE FOR ONE YEAR 

MEGAJOULES 

CALC 

uunt 

c 

C ENfLtC 

fj 

C 

C FNKGY 

ELECTRICAL ENERGY CONSUMED 
IN CURRENT TRIP (BASED ON POWER) 

MtGAJOULES 

CALC 

ttttttt 

ELECTRICAL ENERGY CONSUMED 

MEGAJOULES 

CALC 

tutu 
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C IN CUHRUNT TRIP (BA'olID t)M 

C AMP-HRS) 

i; 


C FA 
r 

VEHICLE ACCELERATION FORCE 

NEWTONS 

CALC 

tuun 

V 

c ri) 

r 

VEHICLE DRAG FORCE 

NEWTONS 

CALC 

mtm 

U 

C ri HLKWd) ARRAY CONTAINING PENALTY/ 

- 

NASA 

1 EFHEKH(I) 

C 

chproveki:;nt factors i-or 



i L.097i8 

c: 

VARIOUS Hf POWER RATINGS 



2 i.OS 

C 

(USED IN FUEL RATE 



3 i. 

C 

r 

CALCULATIONS) 



4 9S4S4S 

u 

c rcR 

f 

GRADE FORCE 

NEWTONS 

CALC 

tut%u 

C IR 
r 

ROLLING RESISTANCE FORCE 

NEWTONS 

CALC 

ttuttt 

1 . 

C rRSCLL 

FACTOR DEPENDENT ON HEKW ONLY 


CALC 

ttitttt 

l! 

USED IN MAKING DIFFERENT SIZED 




c 

r 

HE MAPS FROM A VS KW HE MAP 




L 

C f UMRT 

HEAT ENGINE FUEL RAIL 

LITERS PER 

CALC 

utuu 

C 


SECUND 



t: G(I) 

ARRAY CONTAINING WEIGHTED OBJECTIVE t- 

CALC 

tuun 

i: 

r* 

AND CONSTRAINT FUNCTIONS 




u 

C GAGCSKI) ARRAY CONTAINING PRICE OF 

DOLLARS 

NASA 

I GASCST(I) 

t) 

PETROLEUM FOR EACH YEAR 

PER LITER 


I 0.32 

C 




2 0.35 

C 




3 fl.,58 

c 




4 0.41 

c 




5 0.44 

c 




6 0.47 

c 




7 0.50 

t; 




8 C.S3 

c 




9 0 . 56 

c 




10 0.59 

C GRAv 

ACCELERATION OF GRAVITY 

METERS PER 

CSU 

9.807 

r 


SECOND 



c 

f 


SOOARED 



L 

C GRCi 

CONSVAN'f USED IN CALCULATING 

laLOGRAHS PER 

NASA 


c 

MASS iJr UFImR 

KILOWATT 



c 





grcfr 

CONSTANT USED TO CALCULATE 

DOLLARS PER 

NASA 

2.43 

i- 

COSf III A Cr.AP 

K (LOGRAM 



C GRIlCOI 

COST OF GEAR BETWEEN 

DOLL ARS 

CALC 

rMttM 
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i''ri-Ta’£Hn:AL amd 


V- 


C Gl'!'.r'F(i' 

i RCt Nl GEARS AND Dll ‘ rRF.N rifil. 

- 

NASA 

HAVE 

i.: , 

PlKC> N r FIT ‘ C ICHCY H.V!’ '■ (, 




c; Gra.is) 

(iN I I'L 1 Wi .i . M ZEi\ 0 




c 

Am) 1 0) 




c 





C GKMCSr 

nt)ST Of- GEAR lOWEEN fiLLCIRIC 

DOEE.ARS 

L'ALiJ 


i; 

c 

C HbCj 

MCmiR AND CVT 




CONSTANT USED IN CALCULATING 

KU.OGRAMS PER 

NhoA 

O.fi 

r 

HASS Di' s-.NGlNf 

XW SOUARICO 



C 

CONSlAtf! USED IN CALCUl.ATING 

K j LOGKAHS 

NASA 

2.?A 

i: 

r 

MAS;i OF LMSi'NE 

PEP. K'W 



w 

C fCCF.R 

CONClANi USED TO CALCUlATE 

[101. EARS PER 

NASA 

1.43 

i: 

C 

C HtCOGT 
i: 

c \{\-m 

UFA! EMU f ME CUGC 

K lLOIjRAH 



COST OF HEAT ENGINE 

DOLLARS 

CALC 

ttttut 

POWER RATING Of HEAT ENGINE 

laLOWATTS 

NASA 

65.0 

c 

STiiRi lNG OAi UE OF HEKW WHEN 




c 

1* 

RUN;Pi--F: 




C HLMAIK 

HEAT ENGINE HAINTENAfCE COST 

DOLLARS PER 

CALC 

tttttu 

c 

i; 

C HEHAP( 


'KILUHEVER 



Hf AT ENGINE HAP 

LITERS PER 

NASA 

HAVE 

i; TORiJUE, 

■•iPEED) 

SECOND 



c Hour a ) 
c 

ARRAY WHICH SAVES AT THE END OF 

EACH sn: HIE amp hours out of the 

AMP --HOURS 

CALC 

tutut 

c 

i) 

DATTERY SINCE THE START OF THE STC 

'S 



t; Kwryn. 

HIGHWAY VELOCITY FOR VEHICLE 

MFTERS PER 

NASA 

P.S.O 

n 

U 

c 


SECOND 



C IGOnOD 

BATTERY DISCHARGE FLAG 

.. 

CALC 


c 

(EOUAES E MEANS LtATFERY IS AT OP 




c 

r 

BELOW 807. DEPTH OF DISCHARGE) 




u 

C )AMi'. 

C 

r 

AVERAGE BATTERY CURRENT SINCE 
START UF STC DRIVING 

AMPS 

CALC 

tutttt 

f: IBATt 
z 

BATTERY OUTPUT CURRENT 

AMPS 

CALC 

unn* 

C IBIWX 
r- 

BATTERY MAX. CHARGING CURRENT 

AMPS 

NASA 

200.0 

L 

C })OOLD 

OLD VALUE OF BATTERY CURRENT 

AMPS 

CALC 

nmu 
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CALC 


iuntt 


t: 

C IDtlDd) ARRAY WHICH CAMES STAIF. OF 
C OECrii OF OISOIIARGE FLAO AT 

C EHD or EACH SVC 

C 

C lEfF BATTERY CELL CURRENT (ACCOUNTS AMPS CALC Uttttt 

C FOR CHANGE IN CO.L RE3IOTANCE 

C DUE TO CHANGE IN CELL HASS) 

C 

C iriTMX HAY VALUE OF BATTERY EFFECTIVE AMPS CSU 700. 

C CURRENT 

C 

C .ILL EC ELECTRICAL PATH CALC. H.AC (EQUALS CALC **:»it*!|:* 

i; i (F LLECTRiCAL RATH CALC WAS MADE 

C THIS ITERATION) 

C 

C H ILL BATTERY CURRENT NEEDED TO CHARGE AMPS CAl.C ttUtU 

C BAT fERY TO sac OF 1.0 

C 

C ITLLBK COUNTS NUMBER OF TIMES ACTUAL CALC UUtU 

C i'OUER SPLIT IS CHANCED FROM SPLIT 

C DETERMINED BY SPLIT FUNCTION 

C 

C .TiTECHG BATTERY CHARGING FLAG (EQUAL TO 1 - CALC tUUU 

i. MLANS BATThRY IS BEING CHARGED BY HE) 

C 

C liiWY HIGHWAY DRIVING FLAG (EQUALS 1 IF - CALC UnU* 

C DRIVING ON HIGHWAY RATHER THAN STC) 

C 

C .DCHG COUNTS ITERATIONS IN ATTEMPTiNG TO CALC nUtU 

C MtUCII CALCHLATEO BATTERY VOLTAGE TO 

C BATTERY VOLTAGE ASSUMED IN MAKING 

C A HE CHARGING OF THE BAT (FRY 

C CALCULATTON 

I.' 

C IITER COUNTS ITERATIONS IN ATTEMPTING - CALC UUtU 

L' TT.I MATCH CALCULATE 0 BA (TER Y 

C VOLTAGE TO BATTERY VOLTAGE ASSUMED 

C IN MAKING A BATTERY DISCHARGE 

C CALCULATION 

C 

C IIREGN COUNTS ITERATIONS IN ATTEMPTING TO •• CALC UUtU 

C MAfCH CALCULATED BAKERY VOLTAGE 

C TO BATTERY VOLTAGE ASSUMED IN MAKING 

C A POWER REGENERATION CALCuLAIION ■ 

C 

C .IMAX NUMBER OF POWER SPLIT PARAMETERS - CSU tUUn 

C 

C INOCHG HE CHARGING INHIBIT (LAG - CSU 0 

C (-1 FOR MO HE CHANGING) 
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c 


C IHOELE 

FLAG WHICH INHIBITS FURTHER USE 

- 

CALC 

ntuu 

C 

OF BATTERIFS IF SOC FALLS BELOW 




c 

0.2 DURING HIGHWAY DRIVING 




c 

p 

(-1 FOR INHIBIT) 




t 

C 10VRL.D 

POWER OVERLOAD FLAG (COUNTS NO. 

- 

CALC 

ttuitt 

C 

OF riHES BATTER CEG AND HEAT ENGINE 




c 

p 

CANNOT MEET POWER REtlUIRED BY VEHICLE) 



f. 

POWER SWITCH FLAG (EQUALS i WHEN WE 


CALC 

ttnut 

C 

HAVE ALREADY ATTEMPTED TO SWITCH 




C 

EXCESS HEAT ENGINE POWER TO 




C 

r 

ELECTRIC HOrOR) 




C nERMX 

MAX NUMBER OF ITERATIONS ALLOWED 


CSU/NASA 

iO 

c 

IN ATTEMPTING TO MATCH CALCULATED 




t; 

flAITERY VOLTAGE 10 BATTERY VOLTAGE 




c 

r 

ASSUMED IN MAKING THE CALCULATION 




w 

C IVESTCa) ARRAY WHICH SAVES AT THE END OF 

AMPS 

CALC 

tutu* 

C 

EACH SrC THE AVERAGE CURRENT OUT 




C 

OF THE BATTERY SINCE THE START 




C 

r 

OF rut: src'G 




V 

c ivc;ar 

■INDEX INDICATING PRESENT YEAR IN 

... 

CALC 

uutu 

c 

r 

LIFE CYCLE COST CALCULATIONS 




C JTRIP 

INDEX INDICATING PRESENT TRIP IN 


CALC 

uunt 

C 

p 

Lae CYCLE COST CALCULATIONS 




C KMf'RYR 
r 

KILOMETERS TRAVELED PER YEAR KILOMnERS 

NASA 

16,000. 

u 

C KHCST(I) 

ARRAY CONTAINING PRICE OF DOLLARS PER 

NASA I 

KUCST(I) 

1 

ELEC TRIG TTY FOR EACH YEAR 

HEGAJOOLE 

1 

.om 

C 



•I 

.0125 

c 



.3 

.01,59 

c 



4 

.0153 

c 



S 

.0167 

c 



F. 

oiai 

c 



7 

.0194 

c 



0 

.0200 

c 



9 

.0222 

c 

c 



(0 

.0236 

C LCCflST 

c 

LIFE CYCLE COST 

DOLLARS 

CALC 

uuut 

C LI NOTH 
r 

LENGTH OF THE CURRENT TRIP 

METERS 

CALC 

uutu 

C UNCOH 

LINEAR COMBINATION OF i.II E 

DOLLARS 

CALC 

tuuu 
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C CYCLE COSIj PEfKULEUtl COST 

C AND ACQUISITIOH COST 


L 

C tWlNTC 
V, 

c 

C MARKUP 
c 

VEHICLE. MAINTENANCE COST 
FOP ONE YEAR 

DOLLARS 

CALC 

tttuu 

DEALER MARKUP FACTOR 

- 

CALC 

%uuu 

C MC 
C 

C HCiiTN!') 

C 

C 

p 

MASS Of A BATTERY CELL 

KILOGRAMS 

CALC 

tuuu 

MASS OF THE BATTERY CELL USED 
.IN DEVELOPING THF BAT i FRY 
DISCHARGE CURVE 

KILOGRAMS 

NASA 

9.83 

W 

C MRFCST 
c 

HANUFACTUfilNG COST FOR VEHICLE 

DOLLARS 

CALC 

ittwn 

c in 
c 

p 

TOTAL MASS OF VEHICLE AND 
PROPULSION SYSTEM 

KILOGRAMS 

CALC 

ttttttt 

C MTGRKU 
i: 

C 

c 

POWER RATING OF GEAR BETWEEN 
KLECTRIC HO TOR AND CVT (ALSO 
SEE PHGRIE) 

KILOWATTS 

NASA 

40.0 

C MTKUPK 

i; 

c 

C Ml MAIN 

C 

c 

ELECTRIC MOTOR PEAK KILOWATT 
«Ar.(NG (ALSO SEE PHTMAX.) 

KILOWATTS 

NASA 

• 10.0 

ELECTRIC MOTOR MAINTENANCE COST 

DOLLARS PER 
KILOMETER 

CALC 

nnu* 

C MlRCi 

C 

c 

CONSTANT USED IN CALCULATING 
MASS (IF ELECfRIC MUiOR 

KILOGRAMS 
PER KW 

NASA 

0.6 

C M1RCLR 

C 

c 

CONSTANT USED TO CALCUl. ATE 
TOE COST OF THE ELECTRIC MOTOR 

DOLLARS PER 
KILOGRAM 

NASA 

IS. 95 

c i-nRcsi 
c 

C liTRKU 
c 

f 

COST Of ELECTRIC MOTOR 

DOLLARS 

CALC 

unm 

POWER RATING OF ELECTRIC MOTOR 
(ALGl) SD; PHTPID'’ 

KILOWATTS 

NASA 

20.0 

C HIRSPD AN ARRAY SPECIFYING MAXIMUM 

t: (PMVPCT) EFFICIENCY MOTOR SPEED FOR A 

c GiyKN tiurpur power 

p 

RADIANS PER 
SECOND 

NASA 

HAVE 

ij 

C HI 

i; 

c 

1 

NUMBER OF SERIES CONICCTED 
BATfERY Cai.S NECESSARY TO 
REACH SYSTEM VOLTAGE 

- 

NASA 

84 

C NCORt 
C 

DIMENSION OF GRG ARRAY Z 


CSU 

6900 
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C NCrCLE, 

c 

HUHFIEP or SIC CYCLi-S AiREADY 
SlhUl.ftrED 


CALC 


Lr 

C N(.:Vl..MX 
1. 

HAX NUHK1.E OF STC'S IN ANY 
t’ARfiCULAl'' TRIP 


NASA 

'SS 

c 

C NiirKli 

i; 

(OHEN RlJN2‘i •= E) THE NO. OF 
tU! rEUFNf HEKW VALUES USEB 


CSU 

VARIES 

i. 

C HOSTC 

TRUNCATED VALUE OF SICNO 

- 

CALC 

I:;****** 

L 

c m 

c 

c 

C Hal Ci K 

i 

C 

NUMBER 01 VARIABLE PARAMETERS 
niP URL 

- 

CSU 

VARIES 

NUMBER Of STC'S IN THf FRONT 
OF EACli HUP RtGUIlilNG 
HIGHWAY DRIVING 


NASA 

P8 

u 

c htihl: 
c 
c 
c 

C HTkII’ 

i: 

INTEGER INDICATING NUMBER OF 
TIME TNOIEHENTS 01 DRIVING 
CYCLE ALREADY SIMULATED 


CALC 

tnun 

INDEX INDICATING UHICH TRIP IS 
BEING USED 

- 

CALC 

num 

1 ., 

C NiRlPX 
C 

MAXIMUM NUMBER Of TRIP TYPES 

- 

NASA 

8 

C 

C t'AXt.t. 

POWER REOUIRED AT AXLE 

WATTS 

CALC 

%u%tu 

1 

C 

1' 

BATTERY MAX CHARGING POWER 

WATTS 

NASA 

33,600. 

C fSiOUl 

c 

c 

i; 

I * 

BATTERY OUTPUT POWER 

WATTS 

CALC 

uutn 

HEAT ENGINE POWER USED TO 
CHARGE BATTERIES 

WATTS 

CALC 

ttttttt 

K 

C f'(;yT 

POWER REOHIRED AT CVT INPUT 

WATTS 

CALC 

nttnt 

i; i'CHiHx 
f: 

MAXIMUM POWER RATING OF CVT 

WaTT S 

NASA 

90,000. 

C I'lJpCT 

1‘ 

C 

j • 

DIFFERf N'ilAL OUTPUT POWER AS A 
PERCLNi' OF DIFFTlRENl iAL RATED 
POWER 

PERCENT- 

CAl.C 

nmm 

C I'DRAUl 

c 

DIFFERENTIAL RATED POWER (ALSO 
SF.i: D.I.FFKW) 

WATTS 

NASA 

90,000. 

C I'DRIVF 
C 

POWER KEOUIRED AT DRIVE SHAFT 

WATTS 

CALC 

num 
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C PE.I.EC 
i; 

ELECTRIC POWER RE£|IIIRED AT 
IRPUr ri! COT 

WATTS 

CALC 


u- 

C PLTf:BT(l) ARRAY COM'IAINIWG PE-TRimih 

C CO’Vf'5 rsJR YfciAR I 

p 

DOLLARS 

CALC 

nuut 

U 

C PETROL 

1 ! 
r- 

PETROLt.UH VOLUME USED 
IN CURRENT TRIP 

LITERS 

CALC 


C PLTSTCai ARRAY WHICH SAVEH AT THE END 
C OF EACH ore THE' fOrAL PETRllLEUir 

C VHLllHr USE)) SINCE THE START OF 

C ALL STC'S 

C 

LITERS 

CALC 

utuu 

C ^'ETYR(I) 

C 

C 

r 

ARRAY WHICH CONTAINS THE YEARLY 
PEFfiOLEUM CONSUMPTION FUR THIS 

particular trip 

LITERS 

CALC 

tunn 

t 

c pum 

i: 

c 

C PGfCT 
\: 
c 

c 

C [‘GRATE 
C 

c 

POWER REQUIRED AT INPUT TO GEAR 
BErWEEN evr AND DlrTERENTIAL 

WATTS 

CAT.C 

tutut 

OUTPUT POWER OF GEAR BETWEEN 
COf AND DIFFERENflAL AS A 
PERCENT (IE GEAR RATED POWER 

PERCENT 

CALC 

nuut 

RATED POWER FOR GEAR BETWEEN 
evr ANO DlFEERENriAL (SEE ALSO 
DFGRKW) 

WATTS 

NASA 

90,000. 

C Pfll; 
c 

r> 

HEAT ENGINE POWER REQUIRED AT 
INPUT TO evr 

WATTS 

CALC 

uuut 

u 

C PHLHAX 
C 

C i'HLS'! 

V- 

r 

MAX OUTPUl POWER OT HEAT ENGINE 

WATTS 

CALC 

utuu 

ESTIMATED OUTPUT POWER OF HE 
USED IN OEIERHINING CVT RATIO 

WATTS 

CALC 

uuut 

C PILOSS 

INVERTER POWER LOSS 

WATTS 

CALC 

utuu 

L 

C flNP 

POWER REQUIRED AT OUTPUT OF 
INVERTER 

WATTS 

CALC 

utuu 

L- 

C PlNVrtX 

c 

C PfSGPCT 

C 

C 

n 

MAX OUTPUT POWER OF INVERTER 

WAT TS 

NASA 

40,000. 

POWF-R AT OUTPUT OF GF.AR BO WEEN 
HU FOR AND CVT AS A PERCENT OF 
GEAR RATED OUTPUT POWER 

PERCENT 

CALC 

uuut 

*-f 

C I'HCRTE 
C 

RATED OUTPUT POWER FOR GEAR 
BEfWEEM HOrOfi AND CVT (ALSO 

WATTS 

NASA 

40,000. 
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C SEt: HTGRKW) 

c 

C PHOTOR electric: MOTOR OUTPUT POWER 

C 

C I'HT ESTIMATED OUTPUT POWER OF MOTOR 

C USED IN DETERMINING THE CUT RATIO 

C 

C PMTMAX ELECTRIC MOTOR MAX OUTPUT POWER 

C (ALSO SEE MIMAX AND MTKWPK) 

C 

C PMTPCT MOTOR OUTPUT POWER AS A PERCENT 

C OF RATED POWER (USED FOR RATIO 

C CALC) 

C 

C PHTRTD RATED POWER OF MOTOR (SEE ALSO 

C HTRKU) 

C 

C POVRLD ACCUMHULA1 ION OF SQUARED VALUES 

C OF P SHORT 

C 

C PREGEN POWER AVAILABLE FOR REGENERATION 

C ARER ROAD LOSSES 

C 

C PROAD VEHICLE POWER LOSS DUE TO ROAD 

C LOSSES DURING BREAKING PHASE OF 

C STC 

C 

G PSHORI DIFFERENCE BETWEEN DEMANDED POWER 

C AND POWER PROPULSION SYSTEM CAN 

C SUPPLY 

C 

C PIIKRT(CURRENT/KG) PEUKERT CURVE 

VARIABLE USED TO PRINT VBATT 
AND VBi;HG 


TOTAL. POWER AVAILABLE FOR 
REGENERATION DURING THIS TIME 
INCREMENT IN THE BREAKING 
PHASE OF THE STC 

CVT RAT 10 (OUTPUT SPEED 
DIVIDED BY INPUT SPEED) 

MINIMUM ALLOWABLE CVT RATIO 

MAX ALLOWABLE CVT RATIO 

GEAR RATIO OF DIFFERENTIAL (AXLE 
SPEED 0.1VIDED BY SHAFT SPEED) 


C 

C PVBATT 

i; 

C 

C PWRSTP 
C 
C 
C 

c 

C RCVT 

c 

c 

C RCVTMN 
C 

C l-'CVTMX 
C 

C K!)1FF 

i: 


WATTS 

CALC 

nmtt 

WATTS 

CALC 

tttuit 

WATTS 

NASA 

40,000. 

PERCENT 

CALC 

tuutt 

WATTS 

NASA 

20,000. 

WATTS 

SQUARED 

CALC 

ttttut 

WATTS 

CALC 

tnn%* 

WATTS 

CALC 

ttttm 

WATTS 

CALC 

tttutt 

AHP-HOURS/KG 

NASA 

HAVE 

VOLTS 

CALC 

ttutn 

WATTS 

CALC 

ttiitu 


CALC 

nuu* 

- 

NASA 

0.0 

- 

NASA 

3.4 


NASA 

i.O 
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n o 


c 


C REGFFF 

r 

REGENERATIVE EFFICIENCY 

- 

CALC 

w 

C RlGtN(I) 
C 
C 
C 
r 

ARRAY WHICH SAVES AT THE END 
OF EACH SrC THE AHP HOURS 
REGENERATEO INTO THE BATTERY 
SINCE THE START UF ALL STC'S 

AMP -HOURS 

CALC 

u 

C REGHIN 

p 

HINIHUH POWER WORTH REGENERATING 

WATTS 

NASA 

L 

C RfPAlR 
1* 

VEHICLE REPAIR COST FOR ONE YEAR 

DOLLARS 

CALC 

U 

C REP CHI 

C 

C 

C RtPCVT 
C 

p 

CHOPPER INVERTER REPAIR COST 

DOLLARS PER 
KILOMETER 

CALC 

CVT REPAIR COST 

DOLLARS PER 
KILOMETER 

CALC 

U 

C REPMIR 
C 

P 

ELECTRIC MOTOR REPAIR COST 

DOLLARS PER 
KILOMETER 

CALC 

C REPRHE 
C 

P 

HEAT ENGINE REPAIR COST 

DOLLARS PER 
KILOMETER 

CALC 

C RGEAR 
r; 

RATIO OF GEAR BETWEEN CVT 
ANO DIFFERENTIAL (DRIVE SHAFT 

- 

NASA 

C 

p 

SPEED DIVIDED BY CVT OUTPUT SPEED) 


L 

C RGRH 
C 
C 
C 

C RtIO 
C 

r 

RATIO OF GEAR BETWEEN ELECTRIC 
MOTOR AND CVT (CVT INPUT SPEED 
DIVIDED BY ELEC. MOTOR SPEED) 


NASA 

AIR UEI(5m DENSITY 

KILOGRAMS PER 
CUBIC METER 

CSU 

If 

C RUN21 

A VARIABLE READ IN: 


ARB 

I.; 

0 FOR 1 GCOHP CALL (WITH ONLY A 

SUMMARY SHEET PRINTED) 

c 

1 FOR i GCOMP CALL (WITH B9 L A8 

L SUMMARY SHEET PRINTED) 

c 

c 

2 FOR MAKING A LCCOST GRID WITH HEKW 
AND WGBATT AS THE AXIS 


G 

C 

r 

3 FOR MAKING POWER SPLIT VS CONSUMPTION 
DATA 


1/ 

C RTIRE 

P 

TIRE RADIUS 

METERS 

NASA 

L 

C GAI.VGE 
C 

C GATORG(I) 

i: 

SALVAGE VALUE OF VEHICLE 

DOLLARS 

CALC 

ARRAY Of SCALED TORQUES 
USED IN FUEL KATE CALC 

NM 

CALC 


c 


tnut 

tuttu 

i,OflO. 

iutut 

uutti 

tittttt 

mmt 

tmm 

0.0833 

0.2857 

1.225 


0.50 

untn 
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G SFGO:) 

AKRAY Of WEIGHTING FACTORS 

- 

CSU 

VARIES 

G 

c 

c Gi-xa) 

FUl? ARUAiY G 




ARRAY Of WuIGHTING FAGTORS 
FUR i'ARAHETERS OTTIHIZED BY GRG 

- 

CSU 

VARIES 

G 

c 

c soc 

f' 

STATE OF CHARGE 


CALC 

nuttt 

L« 

G GOGl. 

THE HINIHUM STATE OF CHARGE 

- 

CALC 

tuiut 

c 

ATIAINEO HY (HE BATTERY DURING 




c 

ACCELERATION 




t 

C GOG MAX 

MAXIHUN STATE OF CHARGE TO 

- 

NASA 

0.3 

G 

WillCli HE CHARGES BACTERV 




G 

C S0GSTC<I) 

^ ARRAY liIHICH SAVES THE STATE OF 

- 

CALC 


1' 

GHARGE OF VUE BA TIER Y (VC THE 




C 

r 

END OF EACH STC 




C GIALLtWHT) A KIKCTIOM yOICH GIVES 

% RATED POWER 

NASA 

HAVE 

G 

C 

C 

G SIGDIS 
t 

C S i GNU 

n 

THE MAX POWER l-ROM (iOTOR 
EOR A GIVEN SPEED 




DISTANCE OF ONE STC 

METERS 

NASA 

i,4S0.0 

NUMBER 0( STC'S IN CURRENT TRIP 

- 

CALC 

%%n.uu 

b 

G (AGGEl. 

TIME AT WHICH ACCELERATION 

SECONDS 

NASA 

14.0 

G 

PHASE OF SrC ENDS 




C 

C lAXlE 

r 

AXLE TORQUE 

NEWTON METERS 

CALC 

unu% 

L 

C (EFAKE 

G 

C 

C (GO AST 

TIME At WHICH BRAKING PHASE 
OF STG IE! REACHED 

SECONDS 

NASA 

74.0 

TIME AT WHICH COAST PHASE 
OF STC IS REACHED 

SECONDS 

NASA 

64.0 

G 

C 

C ICVT 

c 

C TGYCL.E 

TORQUE AT INPUT TO CVl 

NEWTON METERS 

CALC 


TIME REQUIRED TO DRIVE ONE 

SECONDS 

NASA 

108.0 

• * 

U 

rt 

DRIVING CYCLE (INCLUDING IDLING) 




c 

C TDiUVE 
C 

C IGl AK 

DRIVE SHAFT TORQUE 

NEWTON METERS 

CALC 

nttm 

TORQUE AT INPUT TO GEAR 

NEWTON METERS 

CALC 

nttm 

G 

BETWEEN CVr AND DirFERENflAL 




c 

C THE 

HEAT ENGINE OUTPUT TORQUE 

KEI4TON METERS 

CALC 

mtm 


c 
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c line. 

SlhULATlON TIME 

SECONDS 

CALC 


%%%%%%% 

i 

C lHCjTO[> 

OUTPUT TORQUE OF ELECTRIC NEWTON METERS 

CALC 


nmn 

t: 

f 

MOTOR 





u 

C TOTBA’f 

FIRST YEAR EQUIVALENT TOTAL 

DOLLARS 

CALC 


ttttut 

C 

BAflEfiy REPLACEhENf COST 





c 

(SEE TOTI.LE) 





c mw 

FIRST YEAR EQUIVALENT BATTERY 

DOLLARS 

CALC 



c 

SALVAGE VALUE (THE AMOUNf OF 





c 

MONEY INVESTED IN THE ITRST 





c 

YLAR NECESSARY fU OBfAIN THE 





c 

AMOUNT Tlii. BATTERIES WOULD BE 





c 

I* 

SiiLO FOR IN THE TENTH YEAR) 





1- 

C T01F.LE 

FIRST YEAR EQUIVALENT TOTAL 

DOLLARS 

CALC 


tuun 


ELEC ffi [CAL ENERI.5Y COSTS, OOTAIMED 





C 

BY COMPUTING THE SUM OF THE FIRST 





t: 

YEAR EQU (VALENT YEARLY COSTS, 





c 

FOR ALL TEN YEARS, WFFJCfi 





c 

ARE AMOUNTS, WHICH IF INVESTED IN 





c 

THE FIRST YEAR AT A FIXED INTEREST 





c 

RATE, YIELD IHE COST OF ELECTRICAL 





c 

ENERGY FOR THE GIVEN YEAR 





C TOIPU 

FIRST YEAR EQUIVALENT TOTAL 

DOLLARS 

CALC 


tttuu 

C 

PETROLEUM COST USD) DURING 





C 

r» 

LIFE CYCLE (GEE TOTELE) 





L‘ 

c. TtirSAL 

FIRST YEAR EQUIVALENT VEHICLE 

DOLLARS 

CALC 



c 

SALVAGE VALUE (CEE T0TB3) 





C IRIPlMai 

ARRAY WHICH CONTAINS THE 

METERS 

NftSA 

I 

TRIPLN(I) 

i: 

LENGTFF OF ALL THE REQUIRED 



i 

10,600.0 

c 

TRIf’S 



2 

30,000. 

c 




3 

SO, 600. 

c 





86,000. 

c 




S 

130,600. 

c 




6 

160,000. 

c 




7 

500,000. 

c 

r 




ii 

(• 

360,000. 

L 

f; ikipwod) 

ARTTAY CONTAINING NUMBER OF EACH 

- 

NASA 

i 

TfilPNO(I) 

s; 

UF fiiE INOrVIOUAL TRIPS PEFT YEAR 



1 

130 

C 




'1 

t- 

35 

c 




3 

S7 

c 





54 

c 




S 

12 

c 




6 

7 
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ooo (-^oo c-joo f-^ooj-joo coo o o ro^o^o 


7 3 
0 i 


C TRIPTH 
C 

r» 

TIHE IT TAKES TO MAKE A 
PARTICULAR TRIP 

SECONDS 

CALC 


U 

C TRIPYR 
C 

r* 

NUMBER OF TRIPS OF A PARTICULAR 
TYPE PER YEAR 

- 

CALC 

tmm 

c 

C TSTOP 
C 

TIME AT WHICH VEHICLE MUST 
STOP AT ENO OF STC 

SECONDS 

NASA 

83.0 

c 

C V(I) 

c 

p 

DRIVING CYCLE VELOCITY ARRAY 
(ENTRIES CALCULAIEO IN PROGRAM) 

METERS PER 
SECOND 

CALC 

ttnttt 

U 

C VI 

c 

p 

VELOCITY AT END OF SIMULATION 
TIME STEP 

METERS PER 
SECOND 

CALC 

nnm 

t 

C V2 
C 

p 

VELOCITY AT BEGINNING OF 
SIMULATION TIME STEP 

METERS PER 
SECOND 

CALC 

tmm 

o 

C VAVE 
C 

p 

AVERAGE VELOCITY OVER ONE 
DRIVING CYCLE TIME INCREMENT 

METERS PER 
SECOND 

CALC 

tmm 

L 

C VBATT 
r 

BATTERY VOLTAGE DURING DISCHARGE 

VOLTS 

CALC 

tmm 

w 

C VBCHG 

BATTERY VOLTAGE DURING CHARGING 

VOLTS 

CALC 

tmm 

L# 

C VliTOL D 
C 

VALUE OF BATTERY VOLTAGE 
FROM PREVIOUS ITERATION 

VOLTS 

CALC 

mtm 

U 

C VCELL 

!' 

P 

BATTERY CELL VOLTAGE 
(WHEN CELL IS FULLY CHARGED) 

VOLTS 

NASA 

2.05 

o 

C VCHG( llATIERV CELL CHARGE CURVE 

C AHP--HRS/KG, 

C CURRENT) 

1' 

CELL VOLTS 

NASA 

HAVE 

'j 

c vdelt 
c 
c 
c 

p 

MAX ALLOWABLE DIFFERENCE BETWEEN VOLTS 

BATTERY VOLTAGE FROM LAST ITERATION 
AND PRESENT CALCULATED BATTERY 
VOLTAGE 

CSU 

1.0 

U 

C VLLO 
C 

r» 

VELOCITY AT CURRENT TIME DURING 
COAST PHASE OF SIC 

METERS PER 
SECOND 

CALC 

mutt 

L 

c VHcoara) array containing yearly 

C COGl OF OPERATING VEHICLE 

r« 

DOLLARS 

CALC 

tmm 

C VHWBTT 

VALUE or BATTERY VOLTAGE AT END 

VOLTS 

CALC 

mtm 
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c 

c 

f 

OF STC'S PRECEEOIKG ALL HIGHWAY 
TRIPS 




C VINT 
c 

BATTERY INTERNAL VOLTAGE. 

VOLTS 

CALC 

tnttt 

C VHINHE 

C 

C 

VEHICLE VELOCITY BELOW WHICH 
NO HEAT EN(5IN£ POWER CAN BE 
DEVELOPED 

MLTERS PER 
SECOND 

NASA 

1.34 

V/ 

C VNEXT 
C 
C 
c 

VELOCITY VEHICLE HAS COASTED TO 
AT END OF CURRENT TINE INCREHENT 
DURING COAST PHASE OF STC 

METERS PER 
SECOND 

CALC 

mmt 

W 

C VOLIOT 
0 
r 

TOTAL HEAT ENGINE FUEL VOLUME 
USED SINCE START Of TRIP 

LITERS 

CALC 

tmm 

w 

C WAXLE 

C 

c 

ANGULAR SPEED OF AXLE 

RADIANS PER 
SECOND 

CALC 

tmtu 

C WCONST 
C 

p 

CONSTANT TERN IN HASS FORMULA 
FOR VEHICLE 

KILOGRAMS 

CALC 

mmt 

w 

C WCURB 
C 

C HCVT 

C 

c 

CURB HASS OF VEHICLE 

KILOGRAMS 

CALC 

tmm 

SPEED AT INPUT TO CVT 

RADIANS PER 
SECOND 

CALC 

mmt 

C UDORIE 

C 

c 

DIFFERENTIAL RATED OUTPUT 
SPUED 

RADIANS PER 
SECOND 

NASA 

116.35 

c uDPcr 
c 
c 
c 

C WDRIVE 

C 

c 

DIFFERENTIAL OUTPUT SPEED AS 
A PERCENT OF DIFFERENTIAL 
RATED OUTPUT SPEED 

PERCENT 

CALC 

ttttttt 

ANGULAR SPEED OF DRIVE SHAFT 

RADIANS PER 
SECOND 

CALC 

tmm 

c «r 
c 

VEHICLE FIXED MASS 

KILOGRAMS 

NASA 

510.0 

C WfCER 
C 
C 
c 

CONSTANT USED TO CALCULATE 
THE COST OF THE FIXED WEIGHT 
OF THE VEHICLE 

DOLLARS PER 
KILOGRAM 

NASA 

3.13 

C WrCQST 

c 

p 

COST OF FIXED HEIGHT OF 
VEiliCLE 

DOLLARS 

CALC 

mmt 

C WGBATT 

C 

C 

c 

MASS OF BATTERIES 

WHEN RUNEl=i OR R = STARTING 

VALUE OF WGBATT 

KILOGRAMS 

NASA 

386.0 
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o O o 


C WGCHGR 

r 

HASS OF THE BATTERY CHARGER 

KILOGRAMS 

CALC 

. mtm 

U 

C UGCHi 

MASS Of IHE CffOPPEfi INVERTER 

KILOGRAMS 

CALC 

mmt 

t 

C WGCLT 
C 

C HGCMT 

MASS or CLUTCH 

KILOGRAMS 

NASA 

0.0 

MASS OF THE CUT. 

KILOGRAMS 

CALC 

ntttt* 

tr 

C UGDIFf 
C 

C UGEAR 

C 

C 

C MGENG 

p 

MASS OF THE DIFFERENTIAL 

KILOGRAMS 

CALC 

mmt 

ANGULAR SPEED AT INPUT TO GEAR 
CEfUEEN CVf AND DIFFERENTIAL 

RADIANS PER 
SECOND 

CALC 

mtm 

MASS OF THE ENGINE 

KILOGRAMS 

CALC 

mmt 

t 

C WGfXTR 
C 

MASS Of EXTRA COMPONENTS IN 
VEHICLE PROPULSION SYSTEM 

KILOGRAMS 

NASA 

0.0 

L 

C WGGRDF 
C 

MASS OF THE GEAR BETWEEN THE 
DIFFERENIIAL AND THE CVT 

KILOGRAMS 

CALC 

mtm 

L 

C UGGRh 

t: 

c 

C WGHOIR 
C 

C MGORTf. 
C 

MASS OF THE GEAR BETWEEN 
THE EI.fXTRIC motor AND THE CVT 

KILOGRAMS 

CALC 

mtm 

HASS OF THE ELECTRIC MOTOR 

KILOGRAMS 

CALC 

mtm 

RATED OUTPUT SPEED FOR GEAR 
BEIHEEN CVT AND DIFFERENffAL 

RADIANS PER 
SECOND 

NASA 

116.35 

C WtPCT 

c 

r. 

r 

OUTPUT SPEED OF GEAR BETWEEN 
CVT AND OIFFERENriAL AS A 
PERCENT OF RATED OUTPUT SPEED 

PERCENT 

CALC 

mtm 

t 

C WIIE 
C 

p 

HEAT ENGINE ANGULAR SPEED 

RADIANS PER 
SECOND 

CALC 

mtm 

U 

C WHEFF 
C 

P 

MOST EFFICIENT ENGINE SPEED 
FOR A SPECIFIED POWER OUTPUT 

RADIANS PER 
SECOND 

CALC 

tttmt 

L 

C WHEHN 
C 

C UKORTE 
c 

p 

MIN ENGINE SPEED 

RAD/SEC 

CSU/NASA 

104.72 

RATED OUTPUT SPEED FOR GEAR 
BETWEEN MOTOR AND CVf 

RADIANS PER 
SECOND 

NASA 

419.0 

C HflOTOR 
C 

r 

ELECTRIC MOTOR ANGULAR SPEED 

RADIANS PER 
SECOND 

CALC 

mtm 

1/ 

C MHT 

C 

C 

MOTOR SPEED COMPUTED IN 
RATIO SUBROUTINE 

RADIANS PER 
SECOND 

CALC 

ttmm 


81 



C WHTRHX 
C 

C Wl‘ 

Q 

MAX SPEED OF MOTOR 

RAD PER SEC 

csu 

I,67S.S 

POWER TRAIN WEIGHT 

KILOGRAMS 

CALC 

ttutt* 

C Wf'L 

C 

C 

C WPROP 

c 

r 

MAXIMUM DESIGN PAYLOAD 
FOR VEHICLE 

KILOGRAMS 

NASA 

41S.0 

HASS Of VEHICLE PROPULSION 
SYSTEM 

KILOGRAMS 

CALC 

uum 

U 

C «3 
C 
0 

MASS OF STRUCTURE AND 
CHASSIS FOR VEHICLE 

KILOGRAMS 

CALC 

%nnn 

C WSCER 
C 

C 

r 

CONSTANT USED TO CALCULATE 
THE COST OF THE STRUCTURE AND 
CHASSIS FOR 1HE VEHICLE 

DOLLARS PER 
KILOGRAM 

NASA 

3. 00 

w 

C WSCOST 

c 

r 

COST OF STRUCTURE AND 
CHASSIS FOR VEHICLE 

DOLLARS 

CALC 

tmm 

u 

C «1 

c 

p 

TOTAL MASS OF VEHICLE AND 
PROPULSION SYSTEM 

KILOGRAMS 

CALC 

tutm 

C HTt. 

p 

TEST PAYLOAD FOR VEHICLE 

KILOGRAMS 

NASA 

207.0 

C X'X 

c 

r 

ARRAY Of PARAMETERS 
OPTIMIZED BY 6RG 

MIXED 

CALC 

tttttu 

If 

C 2 

OVERALL GRG LARGE MAIN ARRAY 

- 

CALC 



c 

STOP 

END 

END 
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APPENDIX B - CALCULATION OF VEHICLE MASS 


This appendix details the method used to calculate the vehicle mass. 
The variable names used below correspond to the names used in the dictionary 
in Appendix A and also in the canputer program. To obtain actual values used 
for any constants, refer to Appendix A. The overall vehicle mass formula is 
the same as that used in [1, 2]. The component mass formulas and the values 
for constants used in those formulas substantially agree with those found in 
[1, 5]. In the following, all masses are in kg. 

Total Mass Calculation 

= constant term in the total mass formula 
= fixed mass of the vehicle (510 kg) 

= design payload (415 kg) 

= mass of propulsion systan 
= total mass of vehicle 
= test payload (207 kg) 

= Bm. 4. 0.23 (WPL) + WF (B-1) 

0.77 

= WCCWST + 

0.77 

The mass of the propulsion system, WPROP, i.s calculated below. 

Propulsion Systen Mass Calculation 

We let (where the masses are in kg, and the conponent power ratings are 
in kw): 


If we let 


then 


wcmsT 

WF 

WPL 

WPROP 

WT 

WTL 

wcmsT 


and 


WT 
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WGBATT 


battery mass 

WGENG 

= 

engine mass 

WGCHGR 

= 

battery charger mass 

WjGRJI 

= 

mass of gear between 
motor and CVT 

WGCLT 

WGCVT 

“ 

clutch mass 

CVTKW 

_ 

power rating of CVT 


CVT Mass 



WGMOTR 



HEKW 

= 

power rating of heat 

— 

traction motor mass 



engine 

WGGRDF 

= 

mass of gear between 

DIFFKW 

ss 

power rating of 



differential and CYT 



differential 

MTRKW 

= 

power rating of traction 

NfTGRKW 

= 

power rating of gear 



motor 



between motor and CVT 

CHIKW 


power rating of chopper/ 

DFGRKW 

=r 

power rating of gear 

WGCHI 


inverter 



between CVT and 
differential 


chopper/ inverter mass 

CHGRKW 

— 

power rating of 

WGDIFF 

— 

differential mass 



battery charter 


The quantities CVTCl, HEC2, MTRCl, CHICl, DIFFCl, GRCl, CHGRCl are all 
constants used below in the mass formulas; their values are given in 
Appendix A. 

Using the quantities defined above, the component masses and the pro- 
pulsion systen mass are given by 


WGCVT 

=S 

CVTCl+CVTKW 

(B-3) 

WGENG 

= 

HEC2+HEKW 

(B-4) 

WaiOTR 

= 

MTRCl+MTRKW 

(B-5) 

WGCHI 

= 

CHICl+CHIKW 

(B-6) 

WGDIFF 

= 

DIFFC1*DIFFKW 

(B-7) 

WGGRM 


GRCl+MTGRKW 

(B-8) 

WGGRDF 

= 

GRCl+DFGRKW 

(B-9) 

WGCHGR 

= 

CHGRC1*CHGRKW 

(B-10) 
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and 


WPROP = WGCVT + WGENG + WCaiOTR + WGCHI + WGDIFF + WGGRM 

+ WGGRDF + WGCHGR + WGBATT + WGCLT (B-11) 

In addition, the curb mass, WCURB, and the power train mass, WP, are 
given by 


WCURB = WT - WTL 


(B-12) 


and 


WP 


WPROP - WGBATT 


(B-13) 



APPENDIX C - CALCULATION OF MANUFACTURING, 
MAINTENANCE, AND REPAIR COSTS 


This appendix details the methods used to calculate the vehicle manufac- 
turing, maintenance and repair costs. The variable names used below corre- 
spond to the names used in the dictionary in Appendix A and also in the com- 
puter program. To obtain actual values used for any of the constants, refer 
to Appendix A. The formulas and values for the constants used in these cal- 
culations are taken frcxn [1, 4]. The component masses used below are calcu- 
lated in Appendix B. 


We 

WF 

WGCVT 

WGCHI 

WGDIFF 


WS 

WCURB 


Manufacturing Costs 

let (where the masses are in kg) 

= vehicle fixed mass WGENG = engine mass 

= CVT mass WGCHGR = battery charger mass 

= chopper/ inverter mass WGGRM = mass of gear between motor & 

CVT 


= mass of differential WGGRDF = mass of gear between 

differential and CVT 

= structure and chassis WGMOTR = mass of motor 
mass 

= curb mass WGCLT = mass of clutch 


In addition, HECER, CVTCER, CHGCER, CHICER, GRCER, DIFCER, WPCER, WSCER, 
MTRCER, CLTCER and AST.ICER are the various component cost estimating ratios 
(CER's) in $/kg. The component costs are then calculated as 


HBCOST = HECER*WGENG (C-1) 
CVTCST = CVTCER*WGCVT (C-2) 
CHGCST = CHGCER*WGCHGR (C-3) 
CHICST = CHICER*WGCHI (C-4) 
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GlttiCST = GRCER+WCXJRM (C-5) 

GRDCST = GRCER*WGGRDF (C-6) 

DIFCST = DIPCER*WGDIFF (C-7) 

WFOOST = WFCER*WF (C-8) 

WS = 0.3*(WPL + WPROP + WF) (C-9) 

WSCX)ST = WSCER*WS (C-10) 

MTRCST = MTRCER*WGMOTR (C-11) 

CLTCST = CLTCER*WGCLT (C-12) 

ASMCST = AaiCER*WCURB (C-13) 


The manufacturing cost is then the sum of these twelve component costs. Note 
that the battery cost is handled separately and is not included in the manu- 
facturing cost (it is included in the calculation of acquisition cost in Sec- 
tion 5). Note also that the above manufacturing cost includes (in addition 
to the component costs) an assanbly cost (ASfiCST), a structure/chassis cost 
(WSOOST) , and a fixed v/eight cost (WPCXIST), 

Maintenance Costs 

It is assumed that maintenance is performed on the heat engine, motor, 
battery and chopper/ inverter. REMAIN, MTMAIN, BTMAIN, CVTMAN and CHIMAN are 
conponent maintenance costs ($/km), HEKW and MTKWPK are component power rat- 
ings (kw), BTLIST, BATCST and BTMKUP are battery list price, battery cost and 
battery markup factor, respectively, WGBATT and WGCHI are the weights of the 
battery and the chopper/ inverter (kg), respectively, and BATCER and KMPRYR 
are the battery CER and the kilometers driven per year, respectively. With 
these definitions, the total yearly maintenance cost, MAINT, in dollars, is 
given as 


87 



MAINT = KMPRYR*(HEfiAIN + MTMAIN + BTMAIN + CVTMAN + CHIMAN) (C-14) 


where 


REMAIN = [0.18 + 0.005*(HEKW/0. 746) 1/160.93 (C-15) 
MTOAIN = [0.06 + 0. 002* (MTKWPK/0. 746) 1/160.93 (C-16) 
BATCST = BATCER*WGBATT (C-17) 
BTLIST = BATCST* (1.0 + BTMKUP) (C-18) 
BTMAIN = (0.004*BTLIST)/160.93 (C-19) 
CVTMAN = 0.063/160.93 (C-20) 
CHIMAN = [ (36.83/2. 0)*SCSlT(WGCHI) 1/(10. 0*KMPRYR) (C-21) 


Most of the above expressions are taken frcm [41. However, the chopper/ in- 
verter maintenance cost is taken from [11, but with the modification that 
half of the cost is assigned to maintenance and half to repair. 

Repair Costs 

It is assumed that repair is performed on the heat engine, motor, chop- 
per/inverter and CVT. The yearly repair cost, REPAIR, in dollars, is given 
as 


REPAIR = KMPRYR*(REPRHE + REIMTR + REPCVT + REPCHI) (C-22) 

where REHIHE, REFMTR, REPCVT and REPCHI are the component repair costs ($/km) 
calculated frcxn (HEKW, MTKWPK, WGCHI and KMPRYR are all defined above) 

REPRHE = [0.28 + 0. 008* (HEKW/0. 746) 1/160.93 (C-23) 

REFMTR = [0.09 + 0.002* (MTKWPK/0. 746) 1/160.93 (C-24) 
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REPCVT = [0.05 + 0.0013*(MTKWPK/0.746) 1/160.93 
REPCHI = [ (36. 83/2. 0)*SQRT(WGCHI) 1/(10. 0*KMPRYR) 


(C-25) 

(C-26) 
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APPENDIX D - CALCULATION OF ACQUISITION COST 


The acquisition cost is a canbination of the manufacturing cost of the 
vehicle drive train and body and the cost of the batteries. The vehicle ac- 
quisition cost, ACQCST, in dollars, is given as 

ACQCST = MARKUP*MNFCST + BATCST*(1.0 + BTMKUP) (D-1) 

where MARKUP is the manufacturing cost markup factor, MNFCST is the manufac- 
turing cost (calculated in Appendix C), BATCST is the battery cost, and 
BTMKUP is the battery markup factor. BATCST and MARKUP are calculated fran 
(the ex-pression for MARKUP is taken fran [4]) 


BATCST = BATCER+WGBATT (D-2) 

MARKUP = (29. 176*10-5 )*MNFCST + 1.40833 (D-3) 

and where WGBATT is the battery weight, and BATCER is the battery cost esti- 
mating ratio (CER). 
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APPENDIX E - VEHICLE SIMULATION 


For each choice of the design variables (battery weight, heat engine 
rating and power split strategy) it is necessary to simulate the vehicle as 
it is driven over the specified yearly driving requirements (Section 4) in 
order to calculate yearly electrical energy and petroleum consumptions, the 
number of deep discharges of the battery per year, and whether at any point 
in any of the trips the on-board energy is insufficient to meet the power re- 
quired at the wheels. This appendix gives the details of how that simulation 
was carried out. It includes a description of how each component in the pow- 
er train was modeled as well as a description of the overall computational 
procedure used to determine the power required frcm each of the on-board en- 
ei^ sources at each instant of time in driving an STC or in highway driving. 
Reference should be made to Section 3.1 and Figure 3-1 for a description and 
diagram of the overall propulsion system and how the various components de- 
scribed below are interconnected. 

E.l COMPONENT MODELS 

In this part of the appendix we describe the methods used to model each 
of the components in the drive train. It should be noted that this is a sim- 
ulation of the steady state rather than the dynamic behavior of the vehicle: 
the component models are therefore steady state rather than dynamic models. 

Battery Model : Of all the canponent models, the battery model is the 
most complex. Because of this complexity it is described in depth in its own 
appendix, Appendix F. 

Heat Engine Model : The heat engine model used was a heat engine map 
[1] which plotted fuel flow rate (liters/sec) against engine speed (rad/sec) 
for a given engine torque (Nm). The map consisted of a set of such curves 
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for different torques as shown in Figure E-1. The map was used to determine 
fuel flow rate corresponding to a known engine output torque and speed. The 
maps were stored as arrays, and double variable linear interpolation was 
used. The actual maps used in this study were the same as those used in [1], 
but with an added degree of complexity. In this study the heat engine rating 
was continually changed in a smooth fashion as the optimization algorithm 
searched for the optimum. Consequently, maps for any rating between a maxi- 
mum and minimum rating had to be easily generated to carry out the optimiza- 
tion. From interpolation studies using the three maps given in [1] for rat- 
ings of 56 kw, 75 kw and 112 kw, it was found that the other maps could be 
generated from the 75 kw map by leaving the engine speed axis alone, and mul- 
tiplying the fuel flow rate ordinate and the torque corresponding to each 
curve by a factor equal to the rating of the engine whose map is sought di- 
vided by 75. Thus, at the beginning of each simulation run, the map array 
for the given heat engine used for that simulation run would be generated and 
stored. Then at each simulation time increment, that array would be used 
with a double variable linear interpolation scheme to determine the fuel flow 
rate during that time increment. 

Chopper/Inverter : The characteristic of the chopper/ inverter that re- 
quired modeling for the instantaneous power calculations was the power loss 
in the device. The model used was the same as that used in [1,2] and is giv- 
en as 

PILOSS = 3.0*(PINV/VBATT) + 0.035*PINV (E-1) 

where PINV and VBATT are the chopper/inverter output power and battery volt- 
age, respectively. The ratio (PINV/VBATT) is an estimate of the current 


92 
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Heat Engine Speed (Rad /Sec) 


Figure E-1 Heat Engine Map 



being drawn from the battery. 

Electric Motor : The traction motor was a brushless, DC, permanent- 
magnet motor operating from a battery voltage of 168 V. It was rated at 20 
kw and had a top speed of 14,000 rp»n. The model for the motor was taken from 
[1] and is an efficiency map with motor efficiency plotted against percent 
rated output power (see BITRTD and BITPCT in dictionary) for a given motor 
speed. The map is shown in Figure E-2. However, in using the map it must be 
recognized that there is a stall-out speed for each percent rated output pow- 
er: for the given output power the motor cannot be run below this speed 
without stalling. The stall-out map is used in selecting the CVT ratio; the 
stall-out map actually used for this motor is given in Appendix G where the 
ratio selection is discussed. 

Continuously Variable Transmission (CVT) : The CVT provides a range of 
ratios (input speed/output speed) which vary continuously from a 0.3:1 over- 
drive to a 3.33:1 speed reduction. The CVT model [1] is an efficiency map as 

shown in Figure E-3 where efficiency is plotted against output power for a 
given ratio. 

Gears and Differential : The differential (with a speed ratio of 1:1) 
and the gears between the electric motor and the CVT input (providing a speed 
reduction from motor to CVT of 3.5:1) and between the CVT output and the dif- 
ferential input (providing a speed reduction of 12:1 fron CVT to differen- 
tial) were all modeled by the same efficiency map shown in Figure E-4 where 
gear efficiency is plotted against percent rated output power for a given 
percent rated input speed. The rated powers and speeds for the different 
gears are PDRATE, WDORTE, PMGRTE, WMORTE, PGRATE and WGORTE, and are found in 
the dictionary in Appendix A. 
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Figure E-2 . Electric Motor Efficiency Map 
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Efficiency 



Figure E-4 Gear Efficiency Map 


E.2 CALCULATION OF POWER REQUIRED AT THE AXLE DURING THE ACCELERATION AND 
CRlilSEi PHASES' t)'F~THg'STC'A^~i^^ HIGHWAV DRIVINCj 

The power required at the axle is determined by summing the forces act- 
ing on the vehicle and multiplying this by the average velocity over the pre- 
sent simulation time incranent. In particular, if the simulation time incre- 
ment is A and current time is t, then the average velocity over the current 


simulation time interval is 

Va = (1/2) [v(t) + v(t+A)] (E-2) 
where v(t) is the vehicle velocity at time t. The required acceleration over 
this time interval is 

a(t) = [v(t+A) - v(t)]/A (E-3) 
If Fq, Fj^, F^ and Fq are the dr^ force, rolling resistance force, accelera- 
tion force and grade angle force, respectively, then 

Fd = C^Ap (Va2/2) (E-4) 
Er = CRflfg (E-5) 
Cr = [10.0 + O.OlVa + (8 X lO-^) x 10-3 (E-6) 
Fa = (1.1) Mp Va (E-7) 
Fq = g sin (3) (E-8) 


where all the forces are in newtons, g is the acceleration of gravity, Mf is 
the vehicle mass (kg), $ is the grade angle at the current simulation time, 
C^A is the aerodynamic drag coefficient multiplied by the vehicle frontal 
area (CdA= 0.6 n^) and p is the air weight density (p = 1.225 kg/m^). The 
factor of 1.1 in the expression for F^ takes into account accelerating the 
rotating inertias in the vehicle drive train. The power required at the axle 
at the current time is therefore given as 
Pa “ (^D ^A ^g) ^a 
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(E-10) 

(E-11) 


The axle torque and speed, and are obtained from 
*A " ''a/% 

Ta = V^A 

where %> is the tire radius. 

E.3 DRIVE TRAIN POWER CALCULATIONS FOR ACCELERATION AND CRUISE PHASES OF STC 

AND FOR HIGHWAY DRIVING 

At each time step in the simulation a calculation is made to determine 
the petroleum and electrical energy consumed in that time step. This is done 
by moving the power required at the axle (see Section E.2) back through each 
component in the drive train. In the following discussion reference should 
be made to the component models given in Section E.l. It should be noted 
that the maps used in the calculations consist of data arrays, and a double 
variable or single variable linear interpolation routine (whichever is appro- 
priate) is used to locate points not in the arrays. 

Knowing the power, speed and torque at the axle, the differential effi- 
ciency map and ratio is used to determine the power, speed and torque at the 
input to the differential. The same exact calculation is then made to deter- 
mine the power, speed and torque at the input to the gear between the CVT and 
the differential. At this point in the power train calculation the CVT ratio 
is determined so that the heat engine is run at its most efficient speed for 
the heat engine's current power level (this method of ratio selection is dis- 
cussed in detail in Appendix G). Once the CVT ratio is known, along with the 
CVT output power, speed and torque, the CVT efficiency map and ratio can be 
used to determine the power, speed and torque at the input to the CVT. This 
CVT input power is then split between the heat engine and the electric motor 
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and batteries using the power split strategy currently in force in the simu- 
lation; no matter what power split strategy is in use, the outcome is an al- 
location of a fraction of the total CVT input power to the heat engine and 
the remainder to the motor and batteries. There is a maximum allowable heat 
engine power, and the simulation off loads any excess to the electrical leg 
if the power split calls for more than the heat engine can deliver. 

Once the heat engine output power and speed are known (the heat engine 
is directly connected to the CVT input so that its speed is the same as the 
CVT input speed), the heat engine output torque can be calculated (power /- 
speed) and used in conjunction with the speed to determine the fuel flow rate 
frcm the heat engine map. This fuel flow rate is then multiplied by the sim- 
ulation time step size to determine the fuel consumption for the present sim- 
ulation time increment. The ccmputation of the total yearly fuel consumption 
using these fuel increments is discussed in Section 5. 

The electrical leg canputation proceeds by first calculating the power, 
torque and speed at the input to the gear between the motor and CVT by using 
the gear output power (from the power split) and speed in conjunction with 
the gear map and the known gear ratio. The motor input power can then be 
calculated using this output power and speed and the motor efficiency map. 
This establishes the power which must be supplied by the chopper /inverter. 

The loss in the chopper/ inverter is modeled by an equation requiring the 
chopper/ inverter output power, which is known, and the battery current, which 
is not known at this point in the calculation. However, it is reasonable to 
assume that (i) if the simulation time step is small, the battery voltage 
does not change significantly from one simulation time step to the next, and 
(ii) the losses in the chopper/ inverter are small. Under these assumptions. 
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the battery current for the present time increment (for use in the chopper/- 
inverter loss equation) can be estimated by dividing the battery voltage at 
the last simulation increment into the present chopper /inverter output power 
and using this result in the chopper/ inverter loss equation. This procedure 
allows calculation of the required battery output power. 

To calculate the battery voltage, current and state of charge requires 
use of the battery model and an iterative procedure. An iteration procedure 
is required because the present battery current is not known, yet this cur- 
rent is required in order to use the battery model to calculate battery volt- 
age and state of charge. At this point in the calculation, the only battery 
quantities which are known are the battery present output power, and the bat- 
tery current, voltage and state of charge which existed at the last simula- 
tion time increment. The iterative procedure used to arrive at present 
values for these last three quantities is to guess at the present battery 
voltage (using the value from the last simulation time step; at start-up the 
battery is fully charged and the first value is the system voltage) and di- 
vide this into the known output power to obtain an estimate of the battery 
current. Knowing the ouput current, the battery model (Appendix F) is used 
to calculate a battery state of charge and voltage. This new battery voltage 
is then conpared to the assumed value: if it differs by less than a speci- 
fied limit (VDELT in Appendix A) the calculation is completed; if the differ- 
ence is greater, the new voltage is assumed to be the battery voltage, a new 
current is computed and tfie process is repeated. An upper limit on the num- 
ber of such iterations was used, but the process always converged in three or 
less Iterations. 

At this point the drive train power calculation is complete and the sim- 
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ulation moves to the next time interval. However, there are a number of 
special features in the drive train power calculation related to limiting 
which should be mentioned. 

First, there are limits on the heat engine power and the electrical en- 
ergy available at any instant of time in the simulation. If the power re- 
quired at the CVT input is split in such a way that either the electrical leg 
or the heat engine leg cannot meet the danand, a calculation is made to de- 
termine the maximum power available from the limiting leg, this is moved 
through any components between the source and the CVT input, and the other 
leg's demand is increased by the appropriate amount (with the other leg's 
power calculation redone if necessary) . If a situation arises in which both 
power sources together cannot meet the total demand, the vehicle is consid- 
ered to be overloaded. The topic of overload is discussed in detail in 
Section 6. 

In the battery calculation it is important to make sure the battery cur- 
rent does not exceed a maximum value (this limit must be scaled each time the 
battery mass is changed). As described above, the battery current is com- 
puted by dividing the battery voltage into the battery output power. If this 
current turns out to be excessive (at low battery voltages) when the battery 
calculation is finished, the available power from the battery must be esti- 
mated (by multiplying the maximum allowable current by the battery voltage 
just calculated), the battery calculation redone, this new battery output 
power moved through the electrical leg components to the CVT input, and the 
difference between the total power and the available electrical power made up 
by increasing the heat engine required power and redoing the heat engine leg 
canputation. 
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E.4 CALCULATIONS FOR COASTING PHASE OF THE STC 

In the coasting phase of the STC the vehicle simply coasts for ten sec- 
onds (see Section 4). The coasting phase always begins at the same cruise 
velocity of 72 km/hr, but the velocity the vehicle has at the end of the 
coast period (and hence the amount of energy available for regeneration) de- 
pends upon the battery weight and heat engine rating (the rating determines 
the engine weight); consequently the coast phase final velocity changes dur- 
ing a simulation run and the coast period must be simulated for each STC. 

The velocity at the end of the coast phase was computed as follows. For 
each time step in the coast phase the current values for the drag, rolling 
resistance and grade forces (using the vehicle velocity at the beginning of 
the time step) were computed (see Section E.2). Adding these, dividing by 
the vehicle mass and multiplying by the time step size (BELT) yielded the de- 
crease in vehicle velocity over the current simulation step. This was re- 
peated for each time step in the coast phase to determine the vehicle veloc- 
ity at the end of the coast phase. 

E.5 CALCULATION OF DRIVE TRAIN POWER REGENERATION DURING THE 

BRAiaNG PHASE OF THE STC 

In regenerating power back into the battery it was assumed that if the 
power available at any time fell below scxne minimum (REGMIN) , regeneration 
would stop and the vehicle would come to rest through normal braking. Also, 
if in the drive train calculation the regenerated power at the input to any 
component exceeded that component's maximum rating, the power passed on was 
assumed to be equal to that maximum, and the difference was assumed to be 
dissipated by normal braking. 
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The first step in the calculation in each simulation time incr«nent is 
to determine the regenerative power available at the axle for that time in- 
crenent. Knowing the velocity at the start of the braking phase, and the 
length of time allowed for braking, the velocity decrease and the average 
velocity for each simulation time step can be determined. The velocity de- 
crease over the time step is used to calculate the decrease in kinetic energy 
the vehicle must undergo in that time step; dividing this by the time step 
yields the vehicle total power reduction over that time interval. The aver- 
age velocity is used to confute the drag, rolling resistance and grade 
forces, and these are then added and multiplied by the average velocity to 
obtain the power lost to the road in the time interval. The difference be- 
tween these two powers is the power available for regeneration. The calcula- 
tion then proceeds in a manner similar to that used in the drive train power 
calculation during acceleration and cruise (Section E.3): the power is moved 
through each canponent (differential, gear, CVT, gear, motor, and chopper /- 
inverter), but with the difference that the component efficiencies multiply 
the output power to obtain the input power, rather than divide it as was the 
case earlier. When the battery is finally reached, the battery input power 
is known, but the current and voltage are not known. An iterative scheme is 
therefore again needed which guesses a battery voltage, uses that to calcu- 
late a battery current, and then uses this current with the battery charging 
model to calculate a battery voltage and state of charge. If the assumed and 
calculated voltages do not agree, the calculated value is used as a new 
starting point and the process repeated. This schene was also found to con- 
verge rapidly (less than three iterations). Care must be taken to check that 
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battery input power and current do not exceed maxinuins (which change with 
changing battery mass) . 
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APPENDIX F - BATTERY MODEL 


The nxxiel used for the battery is based on the model described in [3]. 
The batteries are assumed to be lead-acid batteries with a system voltage of 
168 volts, a cell voltage of 2.0 volts and 84 cells connected in series. 


F.l SCALING THE BATTERY MODEL 

In the optimization studies the battery mass is continually changed. 
This change is interpreted as a change in the mass of each cell; the number 
of cells in the battery is assumed to be constant at 84. To handle this sit- 
uation, all of the battery maps (discussed in the next subsection) were 
scaled (where necessary) so that they represented a battery comprised of one- 
kilogram cells. Calculations were therefore carried out on a per kilogram 
basis, and the results then scaled back up to the actual mass of each cell. 
The actual mass of each cell is equal to the total battery mass divided by 
84. The way in which the maps were scaled should be clear frcxn the discus- 
sions in the following subsections. 

F.2 CALCULATION OF BATTERY VOLTAGE AND STATE OF CHARGE DURING DISCHARGE IN 
HIGHWAY DRIVING AND IN THE ACCELERATION AND CRUISE "PHASES OF AN STC 

Using the battery model given in [3], the battery voltage and state of 

charge at any time instant in the simulation are calculated as follows. 

First, the average current since the start of the trip (assuming a fully 

charged battery at the start of the trip) is calculated as 


Ia 


(1/T) /i(t)dt 
0 


(F-1) 
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where T is the length of time into the trip, and i(t) is the current at any 
time t e [0, T]. In the simulation, was actually computed by dividing the 
ampere-hours (amp-hrs) out of the battery by T (ignoring regeneration). The 
current out of each cell (and out of the battery since the cells are in 
series) on a per kilogram basis is where M^, is the mass of a cell, 
i.e. = Mb/Nq where Mg is the battery mass and N^, the number of cells. 
This average battery output current per cell kilogram is then used in con- 
junction with the scaled Peukert curve. Figure F— 1, to determine the capacity 
per kilogram (amp-hrs/kilogram) , AHj^, available from each cell in the battery 
(and hence the capacity per kilogram available from the battery) if the bat- 
tery is discharged at a current equal to the average current per kilogram 
computed above. The total capacity from a battery with cells of mass is 
then AUg = Next, the net amp-hours out of the battery since the start 
of the trip is computed as 


AH 3 = TIa - AHg (F_2) 

where T and are as defined above, and AHg is the total amp-hrs regenerated 
into the battery (since the start of the trip) during the braking phases of 
the STC's. The battery state of charge, SOC, at the current time is, then, 

SOC = 1 - (AH 3 /AH 2 ) (F-3) 

The battery voltage is obtained from the battery discharge characteris- 
tic shown in Figure F— 2 . This discharge characteristic gives cell voltage as 
a function of both state of charge and cell current for a cell of a given 
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mass. Since the batteries are not used if the SOC drops below 0.2, the cell 
voltage is assumed constant below this SOC. If the mass of the cell changes, 
the internal resistance will change (increasing cell mass is interpreted to 
mean more plate area, hence a smaller interval resistance). This is handled 
by multiplying the battery current by (M3 'pnd/^^c^ where is the cell mass 
and Mstnd ^^e standard cell (taken to be 9.83 kg) used in gen- 
erating the discharge curve in Figure F-2. This scaled current is then used 
with the state of charge to enter the discharge characteristic and determine 
the cell voltage. Finally, the battery voltage is obtained by multiplying 
the cell voltage by the number of cells in the battery. 

F.3 BATTERY CHARGING AT EITO OF TRIP 

This section presents the details of the calculations involved in deter- 
mining the wall plug electrical energy required to charge the battery to a 
fully charged condition after a trip is completed. Calculations for regener- 
ation back into the battery during the breaking phase of an STC is covered in 
Section F.4. At the end of a trip the net amp-hours out of the battery is 
known; net amp-hours is the difference between total amp-hours out and amp- 
hours put back from regeneration during the braking phases of STC's driven in 
the trip. If AHfj is this net amp-hoiu’s out, then the amp-hours delivered per 
kilogram for a cell is AHn/J-Ic, where is the mass of a cell. Assuming a 
charging current of 25 amps, the charging characteristic for a cell (on a per 
kilc^ram basis) shown in Figure F-3 can be entered at A%A^c and an average 
charging voltage for a cell can be determined; multiplying this by the number 
of cells gives an average charging voltage for the battery. In the actual 
simulation, an additional map (AACV in ^Appendix A) was generated from the 
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charge characteristic which gave the average charging voltage (for a charging 
current of 25 amps) as a function of AHj^/Mq. Multiplying this average battery 
charging voltage by the net an^hoiu*s out determines the electrical energy 
needed frcm the charger to recharge the battery. Finally, dividing this en- 
ergy figure by the charter efficiency yields the total wall plv^ energy re- 
quired to recharge the batteries. 

F.4 CALCULATION OF BATTERY VOLT AGE AND STATE OF CHARGE DURING 

REGENERATION IN THE BRAKING PHASE OF AN STC 

Appendix E covers the details of the calculation procedure used to de- 
termine how much of the vehicle power at the wheels available for regenera- 
tion during the braking phase of an STC actually arrives at the battery ter- 
minals at each time increment in the simulation (taking into account effi- 
ciencies and maximum power levels of the drive train components). This sec- 
tion explains how that power alters the battery voltage and state of charge. 
As explained in Appendix E, at each time step in the simulation the regener- 
ated power calculation determines a regenerated power at the battery termin- 
als (PBOUT) and a battery current (lEFF). The current and power to be regen- 
erated must be limited by the maximum current and power that the battery can 
handle, and these limits must be scaled each time the battery mass is changed 
by the ratio of the present battery mass to a standard battery mass. At this 
point in the regeneration calculation, the calculation of the battery SOC is 
carried out; this calculation is exactly the same as described in Section 
F.2, with the exception that the most recent current increment in the calcu- 
lation of 1^ is negative (current is now into the battery) or, if amp-hours 
are used to calculate 1 /^, the regenerated amp-hrs are increased for the most 
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recent simulation increment, the net amp-hoiirs out of the battery is computed 
(total out minus regenerated) and calciilated from these net amp-hrs out. 

The battery voltage is calculated from the charge characteristic given 
in Figure F-3. Specifically, the net amp-hours used in calculating above 
is divided by the cell mass to determine the net amp-hrs per kilogram for a 
battery cell. Then, using both the effective battery current (lEFF) men- 
tioned above and this cell net amp-hrs per kilogram in conjunction with the 
battery charge characteristic in Figure F-3, a cell charging voltage is de- 
termined. Multiplying this by the number of cells in the battery yields a 
battery voltage for the present simulation time increment. In the actual 
battery calculation an iteration scheme is required because the current 
(lEFF) is not really known a priori to allow determination of the battery 
voltage from the charging characteristic (see Appendix E for a description of 
this iterative procedure). 
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APPENDIX G - CVT RATIO SELECTION PROCEDURE 


In this study the CVT ratio was selected (at each time step in the simu- 
lation) such that the heat engine was operated at its most efficient speed. 
The reason for employing this approach rather than the approach used in [1] 
in which the ratio was selected to operate the traction motor at its most ef- 
ficient speed was that the traction motor efficiency was always in the range 
of 90% or above (Appendix E) whereas the heat engine efficiency could be sub- 
stantially improved by proper choice of the ratio. 

The most efficient heat engine speed for a given engine output power was 
determined fron the Brake Specific Fuel Consumption (BSFC) map for a 75 kw 
engine shown in Figure G-1. Each characteristic on the map is engine BSFC 
(in grams/ sec of petroleum flow rate per megawatt of output power) as a func- 
tion of engine output power (in percent of rated power) with engine speed 
held constant. For example, if the engine speed is held constant at 104.7 
rad/sec , the BSFC decreases as the engine output power increases f rcxn zero to 
25% of rated output power (more power output frcm the engine requires a high- 
er engine speed; the characteristic essentially stops at 25% of rated power). 
Also, for a given output power, the speed corresponding to the characteristic 
with the lowest BSFC represents the most efficient operating speed since for 
that output power that particular speed produces the lowest fuel flow rate. 
For example, for output powers below 30%, the most fuel efficient engine 
speed is 104.7 rad/sec; at 50% it is 209.4 rad/sec; at 75% it is 314.12 rad/- 
sec; at 100% it is 418.9 rad/sec. For intermediate output powers linear in- 
terpolation was used to determine the most fuel efficient speed for each heat 
engine output power. The result is shown in Figure G-2. 

At each time increment in the simulation the CVT output power is calcu- 
lated (Appendix E). To continue the drive train power calculation and 
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Heat Engine Speed (Rad/Sec) 
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Figure G-2 ^fc>st Fuel-Efficient Engine Speed vs. Engine Output Power 


116 




determine the CVT input power, the CVT ratio must be known In order to deter- 
mine the CVT efficiency. However, selecting the ratio to run the engine at 
its most fuel efficient speed requires knowledge of the engine output power. 
Since this is not known, the calculation procedure estimates it with the val- 
ue which existed at the previous simulation time increment. Since the drive 
train power and power split are not changing rapidly from one time increment 
to the next, the heat engine output power should not change much either. 
Consequently, the calculation procedure is as follows: after the CVT output 
power is calculated, the previous heat engine output power and Figure G-2 are 
used to determine the best engine speed; this speed and the CVT output speed 
determine a CVT ratio (which is checked and limited to upper and lower bounds 
as appropriate). Knowing the CVT ratio, the CVT input power can be calculat- 
ed and split between the heat engine and traction motor, and the calculation 
proceeds as described in Appendix E. 

Modifications to the above approach occur under several conditions. 
First, if in the previous iteration the heat engine power was zero, the ratio 
is selected to run the traction motor at its most efficient speed. The map 
used for this purpose is shown in Figure G-3. This handles the case where 
the power split calls for all of the CVT input power to be met by the trac- 
tion motor. Second, each time the ratio is selected to run the engine most 
efficiently, a check must be made to make sure this ratio doesn't call for a 
traction motor speed which, at the present traction motor output power level, 
would cause the motor to stall. The stall out curve is shown in Figure G-4: 
for a given percent rated output power the curve yields the minimum motor 
speed to avoid stalling the motor. If the selected ratio would result in 
motor stall, the ratio is decremented by small steps until the stall-out con- 
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dition no longer exists. Finally, if the ratio selected to run the engine 
most efficiently calls for a motor speed which is excessive, then the ratio 
is changed to run the motor at its mazimum speed. 

A copy of the ratio selection program is given in Appendix I. 
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APPEITDIX H - AIRESEARCH P0V7ER SPLIT STRATEGY 


Since ttie A.i Research power split strategy [1] was used in parts of this 
study it is wortli’.-liile including it here for reference. A listing of the 
program is given in Appendix 1. The decision of how much of the CVT input 
power to assign i;o the float engine and how much to the traction motor depends 
on the battery state of cliarge and whether or not the vehicle is accelerat- 
ing. If the acceleration is zero, the state of charge above 20%, and the ve- 
hicle is not in a hill clitif) (CVT input power is < 20 kw) , the power split is 
100% electric. If the vehicle is in a hill climb, the split is 30% electric, 
70% heat engine;. If the acceleration is zero but the battery state of charge 
is below 20% (independent of liill climb), the split is 100% heat engine. 

If the vetiicle is accelerating (independent of whether the battery is 
above or below 20% state of charge) , the split is 30% electric and 70% heat 
engine. However, this split is modified if the CVT input power exceeds 80% 
of the maximum rating of the heat engine: in this latter case the heat en- 
ginp newpr 1.*=: nvT inpii't or Its retting, whichover Is 
smaller, and the remaining required power is delivered by the traction 
motor. 
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APPENDIX I - COMPUTER PROGRAMS 


riM.L 

C ADDING THE D OPTION HILL PRINTOUT DATA (LCC0ST,P0VRLD,S0CL) TO 
C FILE A8 FOR RUH2i=3 (VALUES OVER A SPECIFIED RANGE OF HEKH AND 
C MGBATT) . THE TWO WRITE STATEMENTS FOLLOWING STATEMENT 4109 MUST 
C ALSU 8E COMMENTED ( C ). 

PROGRAM HVPOP(3,100) 

C HYBRID VEHICLE PERFORMANCE OPTIMIZATION PROGRAM 

cmmm%%*mmummmuummmmmmmmmum%nmmm 

c 

C THIS MAIN PROGRAM READS IN INITIAL VALUES, SETS UP THE DRIVING CYCLE 
C ARRAYS (VELOCITY AND GRADE ANGLE), SCALES THE GRG2 VARIABLES, AND 
C CALLS THE OPTIMIZATION PROGRAM. 

C 

cuummumnmmunummmmmmmmmmmmunmum 

LOGICAL SOVRLD,SNOrC, SMOOTH 
REAL KMPRYR 
INTEGER RUN21 
COMMON RUN2I 

COMMON/ALP HA/ALP HAi , ALP HA2 , ALPHA3 

COMMON/CC/CDA,RDIFF,RGEAfi,RGRM,RriRE 

COHhON/MAP/DISCHG(iOi,6),EFFHTR(81,7),HEMAP(6,6), 

iGHEFF(ifll,U),VCHG(13,12),AT0RQ(6),AESP£D(6),ARATI0(i8) 

COMMON/OMISC/CHOOSE, SMOOTH 

COMHON/SCALE/3FX(17) ,3FG(20) ,NN 

COMMON/STC/V(03D,BET(83i),TRIPLN(8),TRIPNO(8) 

COHMON/T IMEC/DLT3TC, DELTHW , TACCEL , TCOAST , TBR AKE , T3T0P , TCYCLE 
COHMON/VARY/ UTL,INOCHG 

CGHHON/YEARC/PETYR ( i 5 ) , ELECYR ( IS ) , ELEC (IS ) , DISCH ( 15 ) 

COMMON GF 

DIMENSION G(20),XX(17),AVELA(iS) 

DIMENSION IBUF(?0) 

DIMENSION HE(3) ,LINE0(16I ) ,LINES(i6i > 

DATA HEKW/6S./, WGBATT/386./, XX(3)/.7/, XX(4)/0./, XX(S)/fl./ 

DATA AVELA/0., 3. 2777,5. 8611, 7. 9166,9. 9722, U. 333, i2. 583, 13. 722, 
614.81)5,15.777,16.75,17.639,18.5,19.333,20.111./ 

CALL LGBUF(IBUF,98) 

C 

C READING INITIAL VALUES OF OPTIMIZATION PARAMETERS 

1 REA))(5,«) HN,(XX(I),I=3,MN),HEKW,WGBATr,DHEKW,DWG8AT,NHEKW,NWGBAT 

2 READ(5,*) DLTSTC, DELTHW, WTL,CH003E,IN0CHG,RUN21,GF 

3 READ (5, 7) SMOOTH 
7 FORMAT(Ll) 

WRITE(6,2000) XX(3),XX(4>,>:X(5) 

201)0 FORMAT! 

6' CONSTANT TERM IN POWER SPLIT FUNCTION=' ,F8 .3,/ 
t,' TERM USED TO WEIGHT ACCELERATION=' ,F3.3,/ 

4' TERM USED TO HEIGHT STATE OF CHARGE--' ,F8. 4) 

WKITE(6,200S)HEKU,«GBATf 
2005 FORMAT!/, 

L' HEKW=',F8.4,/ 
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WGBATT=',F8.4,/) 

IF<SI100TH) «R(TE(6,2fl06) 

2006 FORMAK' SHOOTH=TRUE') 

IF (. NO r. SMOOTH) WRITE(6,2007) 

2007 FORMAT! ' SHOOTH=FALSE' ) 

WRITE(6,2008) GF 

2008 FORMAT!//' GAS FACTOR !GF)=',F7.3,//) 

«R I re ! 6 , 2 1) i 0 ) DLTSTC , DEL TH« , CHOOSE 

2010 FORMAT! 

6' SIC TIME IHCREMEMT !DLTSTC)=',F8.2,/ 

6' HIGHWAY TIME INCREMENT !DELTHH)=',F8.i,/ 

6'!CHOOSE=0 FOR ELECST COMPUTATIONS BASED ON 
6'ELEC!I) 6 ENRGY) CHOOSE=' ,F2.0) 

XX!i)=HEKU 

XX!2)=WGBATT 

C SCALING OF OPTIMIZATION PARAMEIERS 
DO 10 1=1, NN 
10 XX!I)=XX!I)/SFX!I) 

C SET UP 3AE DRIVING CYCLE. 

IMAX=IFIX!TACCEL/DLTSTC)+1 

r0EL=i 

IF!DLTSTC.GT.l.) IDEL=DLTSTC 
DO 20 I=1,IMAX 

C FIND DRIVING CYCLE VEL. AND STORE IN V!I) 
r=FLOAf!I-l)*OLTSrC 
TM!N=INT!T) 

TMAX=miN+FLOAr!IOEL) 

IDOT=INT!TMAX)+l-IDEL 

[TOt'=IBOr+IDEL 

)F!ITOP.LE.lS) GOTO 20 

IfUP=iS 

1B0T=1S-IDEL 

TMAX=14. 

YMlN=TMAX-FLOAT!IDEL) 

20 V!r)=!!T-THrN>*AVELA!ITOP)+!TMAX-T)«AVELA!IBOT)>/!TMAX-TMIN) 
x=i-i 

IM AX=K + IF I X! ! TCOAST-TACCEL ) /DLTSTC ) +1 
DO 30 I=X,IMAX 
30 V!I)=20. 

C DO 40 I=1,IMAX 

C T=DLTSTC#FLOAT!I-l) 

C40 WRITE!i,203fl) C,V!I) 

2030 F0RMAT!'T=',F5.2,3K,'V!I)=',F9.5) 

IMAX=rFIX! rSfOP/DLISTC) H 
DO SO I=i,IMAX 
50 bi:t(I)=o. 

C CALL GRG!Z,,Nf;ORE) 

IfiUN21=RIJN2Hl 

GOTO !101,iSl,202,303,404,50S),IRUN21 
101 CALL GCOMP!G,XX) 
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URITE(8,204C) 

2040 FORhAH' HEKW HGBATT LCCOST POVRLD SOCL ') 
«KtfE(b,20S0) XXa),XX(2),G(i),G<S),G<i3) 

2050 F0RMAT(F9.4,lX,F9.3,i(iX,F9.2),iX,GU.4,lX,F9.7,iX,F9.i) 
iFCoHOOrH) «rtITE(3,i06) 

IF (.NOT.SHOOTH) WRITE(8,iC7) 

106 FORMAfl' ShOOTH^TRUE' ) 

107 FORHAK' SHOOT H=FALSE') 

«Rire(3,083) G(2) 

883 FORMAT!// TtHPET=',F9.3) 

srop 101 

c 

c — 

151 CALL GCOMP(G,XX) 

STOP ISl 

C — 

C=:r===r== 

C UDH = 1 WHEN OVERLOAD REGION LIES BELOW CURRENT HEKW OVERLOAD THRESHOLD 
I'UKVe (SEARCH FOR FIRST ADMISSIBLE RUN, THEN GEF CURVE 
C =1 WHEN OVERLOAD REGION LIES ABOVE CURRENT HEKW OVERLOAD THRESHOLD 
C CURVE (SEARCH FOR FIRST OVERLOAD, THEN GET CURVE) 

C UDW = 1 WHEN WGBATT INCREASES IN THE SEARCH 
C --1 WHEN WGBATT DECREASES IN THE SEARCH 

C SQVRLD = TRUE TO SEARCH FOR AN OVERLOAD AFTER WGBATT HAS BEEN CHANGED 
C SNUTC = TRUE TO SEARCH FOR A NEW OVERLOAD THRESHOLD CURVE 
202 TWGBAT=WGBATT+DWGBAT*FLOAT(NWGBAT-l ) 

BWG8AT=WG8ATT 

WRITE(8,2060) 

2060 FORMAT!' HEKW WGBATT LCCOST POVRLD SOCL') 

UOW=i 

UDH=-1 

SNOTC=. FALSE. 

SOVRLD=. FALSE. 

DO 190 K=l,5 
1=1 

60 CONTINUE 

WRirE(i,2060) 

WRITE(9,2060) 

C 

C SEARCH FOR NEXT OVERLOAD AT WGBATT=30fl 
IF(K.EtJ.l) GOTO 90 
J=K“2 

1F(J/2*2.EQ.J) UDH=-UDH 
IF(K.GE.3> UDW--UDW 
HGBATT=3fl0. 

XR=IIE( 3 ) 

SNOTC=.TfiU£. 

SOVRL0=. FALSE. 

IF(K/2«2.EQ.K) GOTO 70 
SNUrC=. FALSE. 
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GOTO 9(i 


C SEARCH FOR FIRST HtM OVERLOAD CURVE WITH CONSTANT STEP SIZE 
70 XR^XR-DaXWOH 
GOTO ISO 

80 IE(G(S).LE.O.) GOTO 70 
3N0n>. FALSE. 

C 

C INCRENENT HATTERY HEIGHT i TEST IF ITS OUT OF RANGE 
90 m?.) =HGOATT 1 0 «G8AT#FLOA r< i-1 ) 3:UDW 

IF(XX<2).LE.IWGBAT) GOTO iflO 
I^MHGBAT 
GOTO 180 

100 IF(KX(2).GL.B«GBAT) GOTO 110 
1>-NHGL!A r 
GOTO 180 

lifl IF(SNOTC.OR.K.NE.l) GOTO ISO 
Ti-a.NE.i) GOTO ISO 
SOVRLD^. FALSE. 

DlILX=2. 

UDH=1 
XR--1S . 

GOTO 140 
C 

C SEARCH FOR OVERLOADS t, DOUBLE STEF SIZE WITH FAILURE 
120 IF(GTS).LE.0..> GOTO 130 

SOVRLD=. FALSE. 

GOTO 160 
C 

C DOUBLE STEP SIZE 
130 Dllx:=DELX*2. 

XR=XR-DELX*UDM 

3UVRL0=.TRUf:. 

GOTO ISO 
140 CONTINUE 
ISO XX(1)=XR 

CALL GCGt1F(G,XX) 

HRITE(1,2070) XX(1),XX(2),G(1),G(S) 

HRIfE(9,2O70) XX(1) ,XX(2) ,G(1 ) ,G(S) 

2070 F0RMAT(F9.3,1X,F9.3,1(1X,F9.2),1X,G11.4,1X,F9.7,1X,F9.1) 
IF(SOVRI.D) GOTO 120 
Jf (SNOTC) GOTO aO 

IF(G(S).LE.O..AND.DELX.LE.DHEKW) GOTO 170 
i:F(U(S).Gr.O.) GOTO 160 
DELX=^AINT(DELX*5.)/10. 

IFIDELX.LE.O.) DELX=OHEK«/2 . 

XR=XL+DELX«UDH 
GOTO 140 
160 XL=XR 

XR=XR+DELX*UDH 
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GOTO ISO 
170 CONTINUE 

«R1TE(8,2C70) XX(i ) ,XX(2> ,G(i ) ,G(S) ,G(i2> 

IK<XX(2).Eq.iOO.) HE(i)=XX<i> 
lf(XX(2).EQ.200.) HE(2)=XX(1) 

£K(XX(2).EU.300.) HE(3)=XX(i) 
iUO l=I+i 

IFa.LE.NWGBAT) GO TO 60 
190 CONTINUE 
STOP 202 

C-:====t== 

C-i"K++l+4 

303 CONTINUi: 

D MKI1E(8,2090) 

D2090 FORMAT!' HEKW UGBATT LCCOST POVRLD SOCL ') 

00 lOO I=l,KHi;KU 

XXd )=HEKM-DI€KW#FLOAT(I-l ) 

00 109 JdjNMGOAT 

XX ( 2 ) - WGB Al T-f DMGBAT SFLOAT ( J-- 1) 

CALL GCOIiP<G,XX) 

D WRITE(8,3000) XXd ) ,XX(2) ,G(i) ,G(S) ,G(13) 

03000 F0RMAr(F'y.3,lX,F7.3,t(lX,F9.2),iX,G11.4,lX,Fy.7,lX,F9.1) 
WRITE(8,3003) XX(l) ,XX(2) ,G(2> ,G(1) 

3003 F0RMAT(F9.3,3(5X,F9.3)) 

C LINE0(J)=1H. 

C II (G(S) .GT.O.) LII€U(J)=lHtf 

C IF(G(2).GT.1072.) LIME0(J)=1H9 

C LINES! 

C IF!G!13).L1..2) LINES! J)=1.H+ 

109 CONl'INUt-: 

C WRITES POVRLD GRID. COMMENT IF USING D OPTION (DATA PRINTOUT). 

C HKIfL!H,3005) !lIN£0!J) ,J=1,HWGBAT) 

C WRITES SOCL GRID. COMMENT IF USING D OPTION !DATA PRINTOUT), 
i: WfiITE(9,300S) (LINES! J),J=l,NWGiIAT) 

C WRITE!?, 300S) !LINEO(J) , J=1,NWGBAT) 

300S FORMA r!i;X,161Al) 

108 CONTINUE 
0 WRirE!6,3;.'.0) 

320 FORMAT!' FILE A8 CONTAINS THE POVRLD GRID:',/ 

6,' FUR NO OVERLOAD (POVRLD. LE.O )' ,/ 

i,' 9 WHEN AN OVERLOAD OCCURS (POVRLD. GT.fl)',/ 

6, 'FILE B9 CONTAINS THE SOCL GRID:',/ 

6,' ' WHIN ELECTRIC POWER IS AVAILABLE DURING ACCELERATION',/ 

6,' 5- WHEN ELB:. power IS NUr AVAILABLE DURING ACCELERAflON' ) 

STOP 303 
C4++i til 

404 DO 409 1=1,21 

XX!3)=.I)S»FL0AT!I-1> 

409 CALL GCOMP(G,XX) 
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STOP 404 
ctmit 
mmi9 
soy CONTINUII 

DO SIO I=i,48i 
XX(2)=WGBA rrti . 0*FLOAT( I-i ) 

IF(XX<2).LE.273.) XXa)=37.0-0.02S64<XX(2) 

IF < XX( 2 ) . G r . 273 . ) XX ( i )=77 . 2334-0 . 22677S6*XX ( 2 ) 
Ui.9819E-04m(2)*XX(2) 

CALL GCOHP(G,XX) 

HRITE(8,5H0) XXa),XX(2),G(l),G(7),G(8),G(9) 

S20 F0RHAT(F9.4,S(3X,F9.4)) 

SiO CONTINUE 
STOP SOS 

cesgggoo 

END 

END 
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nH4,L 

SUBROUTINf GCOHP(G.X'X) 

cmmmunmtmmmnnmmmmmmmmmtutmmmtmm 

L‘ 

C THIS EVALUATES THE CONSTRAINT AND OBJECTIVE FUNCTIONS FOR GRG2 GIVEN THE 
C IHE GRG2 FARAHETERS, WHICH ARE ROWER RATING OF THE HEAT ENGINE, BATTERY 
C WEIGHT, AND THE POWER SPLIT PARAMETERS, 

C 

cmmmmmunnmmmmtmmmmmmmmmtmtmtmtm 

REAL Mr,MHFCSr,KMPfiyR,Mri(MPK,LCC03T,l1AINTC,LINC0« 

DIMENSION G(20),XX(i7) 

ItffEGER RUN2I 
COMMON RUN2I 
COMMON GF 

COMHON/ALPHA/ALPHAI , ALPHA? > ALPHAS 
COMMON/MOTOR/HTKWPX 

COMMON/RSPLT/ACCEL,SOC,PCVT,CXaS),VAVE,PHEMAX, PHE,PELEC 

COMMON/SCALE/SFXI f/) ,SFG(?C) ,NN 

COMMON/YEARC/PETYR ( 15 ) ,ELECYR ( 15 ) ,ELEC ( 15 ) , DISCH (15) 

COMMON BTC0ST,6ATRYR,BATLF£ 

DIMENSION £LECST(10),PETCST(10) 

DATA KMPRYR/16D00./, WP/20./ 

C 

C THIS PROGRAM TAKES VALUES FOR THE IC ENGINE SIZE, BATTERY WEIGHT 
C AND POWER SPLIT FUNCTION PARAMETERS, AND COMPUTES THE VALUES OF 
C THE OBJECTIVE AND CONSTRAINT FUNCTIONS 
C 

HEKW=XX(1)*SFX(1) 

WGDATI>-XX(2)#SFX<?) 

TMAX=NN-2 
DU 115 I=-1,IHAX 
115 C:X(I)=XX(I+2)»SFX(I+2) 

C CX IS VECTOR OF UNSCALED POWER FLOW PARAMETERS 
CALL WGHT(HEKW,WGBATT, WT,M1,WGCHI,WP) 

CALL MFCST(WGBATT,WGCHI. HNFCST) 

CALL AQCST(MNFCST,WGBATT, BATCST,ACQCST) 

CALL MAINr(llEKW,HfKWPK,8ArCSr,KMPRYR,UP,WGCHI, BTLISf ,MAINTC) 

CALL REPAR(HE((W,MTKWPK,WP,WGCHI,KMPRYR, REPAIR) 

C TRAVL USES POWER FLOW PARAMETERS 

CALL 1RAVL(HEKH,WGBATT,M1, PETYR,ELECYR, ELEC, DISCH, POVRLD,PETV10. 
«iSOCL,BGPCVT) 

IF(RUN21.LE.l) WRITE(6,425) MT,WGBATT,WP, HNFCST, ACQCST,MAINTC, 
4fi£PAIR,HEKW,P0VRLD 

425 FORMAT! SUMMARY DATA ummumu'/ 

TOTAL MASS (MT) =',F9.3,/ 

BATTERY MASS (WGBATT) =',F9.3,/ 

POWER TRAIN WEIGHT (WP) =',F9.3,/ 
fx' MANUFACTURING COST (HNFCST) =',F9.3,/ 

ACQUISITION COST (ACQCST) =',F9.3,/ 

YEARLY MAINTENANCE COST (MAINTC) =',‘F9.3,/ 
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A ' YEARLY REPAIR 0)3 f (REPAIR ) ' ,F9 . 3 ,/ 

i' Ht POWER RATING (HEK«)=' ,F8.3,/ 

{.'SUN OF POWER OVERLOAO SQUARED (POVRLD)=' ,F3.3) 

CALL LIFE(PETYR,ELECYR, ELEC, DISCH,BTLIST,ACQCST,HAINTC, REPAIR, 
iBAiCST, LCOJSr,PErCSr,£L£CSr, fOTPET) 

C 

C NY CONVENTION IS TO SEPARATE GIVEN AND RETURNED VALUES BY A BLANK 
C AND TO COMHUNICATE CALC VALUES BY AKGUilEHr LIST AND NAPS BY CONHON 
C STATEMENTS 
C 

LINCOM=ALPHAitLCCOST+ALPHA2fJOTPET+ALPHA3*ACQCST 

G(l)=l.CCOST!lSFUa) 

G(2)=T0TPET)|:SFG(2) 

G(3)-=ACaC3T*SFG(..() 

G(^)--=LINC0M«SFG(4) 

G(S)=P0VRL0*3FG(S) 

G(7)=BTC0ST 

G(8)=DAfRYR 

G(9)=BATLFE 

G(i3)=^30Cl. 

G(i4)=BGPCVT 

RETURN 

RtrURN 
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■3UEIR0UTIHE TK AVLUIEKU , WGBATT .. HT, PETVR , E.I.LC YR , ELEC , DISCH , POVRLD , 
?<PErVlD,30CL,BGPCVT) 

cuumtmmtnummnnmnuumummmmmmmmumum 

c 

(: THIS SIMULATES BOTH STC AND HIGHWAY DRIVING FOR THE VEHICLE 
C IT CALCULATES THE YEARLY FUEL CONSUMPTION, ELECTRICITY USAGE, AMD 
C THE NUMBER OF DEEP DISCHARGES OF THE BATTERY. 

C 

cuutmmmmum%muunnmmm%ummummmnmmmm 

LOGICAL SMOOTH 
INTEGER RUN2I 
REAL INT£0,INTNO 

REAL MT , lAVE , IBATT , IVESTC ,KMPRYfi , IBMAX , lEFFMX , IBTOLD 
REAL K WCS r , LENG ri I , M f , MTGRKU , M IXWP K , MC , NC , I EFF , MCSTND 
DIMENSION PETYRaS),ELECYR(i5),ELEC(iS),DISCHaS) 

0IMEN3 CON EAKWHR < i’j) , AEKHAdS ) , AEKHP T IS) 

DIMENSION AHEKW<4),FFHEKW<4),SAT0RQ(6) 

COMMON RUN2I 

(:OiIMON/CC/CDA,RDIFF,RGEAR,RGRH,RTIRE 
C0MM0N/CHARG/AACV(3) , ACR (Q) 

COMMON/EFCVT/APCVnS) ,ARCVT(7) ,EFF(8,7) 

COMMON/LIMI r/ NTRIPX, IfERMX 

COHMON/M AP /D ISCHG( 1 01,6), EKFMTR (81 , 7 ) , HEMAP (6,6), 

60REFF (101,11), VCHG( 13 , 12 ) , A TORO ( 6 ) , AESPED ( 6 ) , AR ATIO( 1 0 ) 

C0hMQN/0HI3C/CH00SE,SM00TH 

COMHON/PHY/GRAV,RKO 

COMHON/PR INT/APETT( 15) , AELET A( IS) , AELET P ( 15 ) 
COHMON/RRArU/NriME,WG£AH,HCVr,PMOTOR,PMTRTl),PGEAR 
( :OMMON/R SP LT/ACCEL , SQC , PCVT , CX ( IS ) , V AVE , P HEM AX , P HE , PELEC 
ClJMH0N/STC/V(831 ) , BET(831 > , TR I?LN(8 ) , TR IPNO(8 ) 

COMHON/T IMEC/DLTSTC , DELTHW , TACCEL , TCOAST , TBR AKE ,TSTOP , TCYCLE 
COMHON/rRIt'/CHA8GE(36>,E3TOTL(.16),HOUT(36),IDOD(36> , 

(aVESTC(36) ,PETSTC(36) ,S0CSTC(36) ,REGEN(36) 

CUMHON/VAKY/ UrL,iNaCHO 
DATA AHEKW/37.,S6.,75.,112./ 

DA f A FFHEKW/1 . 0V718,1 . 05, 1 . , 0 . 954545/ 

DATA ALLSTC/SOOOl./, HUYVEL/2S. 0/, STCDIS/14S0 . 0/, CHGEFF/O . 05/ 

OA f A NC/84 . /,WDOH rE/116 . 35/,UG0RTE/l 16 . 35/,UM0R TE/419. / 

DATA N3TCFR/28/, PCVTMX/90000 ./, NCYLMX/35/ 

DATA PDRArE/9000C./, PINVMX/40COO./ 

DATA PGRATE/90000./, PMGRTE/40000 ./, REGMIN/lOOfl ./ 

DATA PMTMAX/40000./, S0CMAX/.3/, VCELL/2.057 

DATA VDELT/1./, MCSTND/9.83/, PBMAX/33600 . /, IBHAX/200./ 

DA TA EFFBC/ . 75/ , IEFFMX/7 0 0 . / 

C CALCULATE THE MASS OF A BATTERY CELL 
MC=WGBAT1 /NC 
C INITIALIZE VARIABLES 
I30DOD=0 
AHOUT=fl . 

AREGEN=D. 
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SOC=i.O 
SUCL=i . 0 
liGPCVT=(l. 

POV«L0=0 . 

)NOELE=0 

IH«y=0 

rSWTCH=0 

VBCHG=200. 

Vi]LTOr=0 . 

AHCHG=0 . 

IOV«LD=0 

]FLLBK=fl 

E010T=U. 

VBATr=NC*VC:ELL 

C SET TIME, TIME STEP AMD CYCLE COUNTER FOR DRIVING A SEGUENCE OF STC'S 
TIME=0 . 

Ht:YCLE=l 

DELT^DLTGTC 

PHEMAX=I000.K:HEKU 

FR SCLE=HEK« /7S . »SV INT < FFHEKW , HEKM , AHEK W , A ) 

01) S 1=1,6 

5 SATORO(I)=ATORG|(I)*(HEKW/75.) 

C 

C ACCELERATION AND CRUISE PHASE 

C AND HIGHWAY DRIVING 

C CALC. AVERAGE VELOCITY AND ACCELERATION OVER CURRENT TIME STEP 
10 Vi=V(NTIM£+l) 

VP.=V(NTIME) 

VAVE=<Vi+V2)/r!. 

ACCEL=(V1-V2)/DELT 

BETA=DET(NfIME) 

C CALC DRAG FORCE OF VEHICLE 

21) FD=CDA<rtHU»VAV£*VAVE/2. 

C CALC ROLLING RESISTANCE FORCE 

C«=(iO .+1) .01#VAVE+0.000flU*VAVE»VAVE)/1090 . 
fR=CR«HT*GRAV 

C INiriALIZe POWER SWITCH FLAC 
ISMTCH=0 

C CALC. GRADE FORCE 

FGR=MY»GRAV»SIN(BETA) 

C CALCULATE VEHICLE ACCELERATION FORCE 
)A=i.i«MT!|t ACCEL 
C CALC. POWER REG'O AT AXLE 
P AXLE= ( FD+( R+ FGR +FA ) »VAVL 
C CALC, AXLE SPEED AND TORQUE 
WAXLE=VAVE/RTIRE 

l'6V;i PsPAVI t-VUAyi F 

C CALC. DIFFERENTIAL INPUT SPEED, TORQUE, POWER 
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WDRIVL'^UAXLE/ROIFF 

(•»PCT=(PAXLE/PDRATE)«00. 

UDPCT=(«AXI.E/WOORrE)»iOO. 
|•DRIVE=PAXLE/lMTEQ(GREFF,PDPCT,«DPCT,^flO.,^00 , 

TDRIVE=PDRIVE/H»RIVE 

C CALC. TORQUE, SPEED, POWER AT INPUT TO GEAR BETWEEN CVT AND DIFFERENTIAL 
WGEAR=WDRIVE/RGEAR 
PGPCT=(PDRiyE/PGRATE)*iOfl . 

WGP C r= ( UDR I VE/WGQR TE > «i0 0 . 
l’GEAR=PDRIVE/INTEQ(GREFF,f'GPCT,WGPCT,iOO.,iOO., 

TGEAR=PGEAR/WG£AR 
C DETERHINE CUT RATIO (fiCUV) 

CALL RATIO 

C CALC. TORQUE, SPEED, POWER AT CVT INPUT 

PCyr=<PG'r;Afi)/(INTNQ(EFF,(PGEA»/1000.),RCVT,i00.,3.4, 

i.APCVT,ARCVT,8,7)) 

IF ( PCV r . GT . BGPCVf ) BGPCVr=PCV !' 

UCVT=WGEAR/RCVT 

TCVT=PCVr/«CVT 

C DETERHINE POWER SPLIT (RESULT IS PHE AND PELEC) 

CALL SPLIT 

C IS USE OF BATTERIES PROHIBITED 
IF( INOELE.NE.l) GOTO 33 
PELEC=0. 

PHE=PCVT 

C IS ELECTRIC POWER NEGLIGIBLE 
33 IFCPELEC.LT.Sfl.) GOTO 3i 
C IS BATTERY 8flX DISCHARGED 
IFISOC.GE.fl.iJ) GO TO 32 

31 PHE=PHE+PELEC 

If (PELEC. GT.O.) IFLLBK=IFLLBK+i 
PELEC=0 . 

PNOTOR=0. 

PINV=0. 

PBt)UT=0. 

IBATT=«. 

IELEC=0 
GU fO 80 

32 ir(SQC.GE.SOCHAX) GOTO 30 
IT(IHECHG.NE.l) GOTO 30 
GOTO 31, 

C CALC. INPUT QUANTITIES TO GEAR BETWEEN MOTOR AND CVT 
30 IEL0:=I 

WMOTOR=WCVT/RGRM 
P MGP C r- ( P ELEC/P MGR CE ) *1 0 0 . 

«MGPCT-(WCVT/WMORTE)*100 . 

PHOmH=PELEC/INT£Q(GReFF,?MGPCT,«MGPCT,10fl. ,t00 
(.i.,-I.,i0.,-10.,i0I,U) 
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lFGi=PELEC/PMOTOR 
TMUr08=P MOTOR /WHO rOR 
C DOES MOTOR OUTPUT POWER EXCEED PEAK 
li-(PMOTOR.LE.PMrMAX) GOTO SO 
C: SET MOTOR POWER AT MAX 
|■HOTOR=PHrMAX 

C CALC. NEW MOIOR TORQUE 
TMOrO«=P MOTOR /WHO rUR 

C CALC. NEW ELECTRIC POWER AT CUT INPUT 
40 PMGPCT=(PMOIOR/PMGRTE)tiOO . 

PELEC=PMmOR*INT£Q(GREFF,PMGPCT,WMGPCT,iOO.,iOO., 

LTGi!=PaEC/PMOrOR 

)rLLBK=IFLLBKH 

C SET ENGINE POWER TO MEET TOTAL REQ'D POWER 
PHh>-PCVr--PELEC 

C CALC MOTOR INPUT (INVERTER OUTPUT) POWER 
SO PMrPCT=(PMOTOR/PHTRTD)!fiOQ . 

P INV=PMOTOR/lNTEO(EFFMTR ,PMlPCT,WMOTOR ,200 . , 1466 . , 

^<2 . S , -2 . S , i 46 . 61 , 4 . 82 , 31 , •/ ) 

Ef-M1=PMOTOR/PINV 

C DOES INVERTER OUTPUT POWER EXCEED MAX 
IF(PINV.LE.PINVMX) GOTO S5 

C SET INVERTER POWER AT MAX AND COMPUTE NEW MOTOR OUTPUT POWER 
PINV=PINVMX 

PMIPCT=(PlNV/PMTRTD)*lflO. 

PMOTOR=P INVIINTEO(EFFHTR,PMTPCT,WMOTOR, 200., 1466. , 

62. S, -2. 5,146.61,439. 82, 81, 7) 

EFMl=PMOTGR/PINV 

C COMPUTE NEW ELECTRIC POWER AT CUT INPUT 
P HGP C T= ( P MO TOR /P MGR TE ) *1 0 0 . 
|■Ll.EC=PMOTOfi»INTEQ(GREFF,PMGPCl,WMGPCT,100.,100., 

M.,-1. ,10. ,-10. ,101,11) 

C COMPUTE NEW HEAT ENGINE POWER 
PIIE=PCVr-PELEC 
;ifLLBK=IFLLBK+l 

C CALC INVERTER LOSSES AND INPUT POWER USING OLD VALUE OF BATTERY VOLTAGE, 
C AND RESET ITERATION COUNTER 
55 P ILOSS=3 . 0* (P lNV/VBATT)+0 . 03SIP INV 

P80UT=PINV+PILOSS 
irrER=o 


t DISCHARGING BATTERY VOLTAGE LOOP 


C CALCULATE BATTERY CURRENT USING BATTERY VOLTAGE FROM LAST ITERATION 
41 1BATT=PB0UT/VBATT 
C SAVE OLD VALUE OF BATTERY VULTAGf 
VBTOLD=VBATT 
PV8ATV=VBArT 
)B10LD=IBATT 


133 



C Oi^LC 8ATrt'«Y AMH-'HRy OUT 

AHOUT=--AHOUT+(lBftTTIDELr/3600 . > 

0 I!i THIS MIGHHAY DRIVING 
IRIHWY.NE.I. ) GOTO 57 

C CALi: AVERAGL BAff CUfiRENT KOR SIC AMO HIGHWAY DRIVING 
1 AVE= ( 36 0 0 . < AHOUT ) / ( TIME+DELT ) 

GOrO 50 

C CALC AVERAGE BATT CURRENT SINCE START OF STC DRIVING 

57 IAV'c=( 3600 .SAHUUD/UriHE+OELn+FLOATINCYCLE-DiCTSTOR) 

C FIND CELL CAPACITY PER KILOGRAM 

58 CAPKG=PUKRfaAVE/MC) 

C CALC AHP-HRS PER CELL 

CAPCTY=CAPKG*Mi: 

C CALC ACTUAL CAPACITY WITHDRAWN AND STATE OF CHARGE 
AAC rL=Al (()U r-AREGEN-AHCHG 
SUC-I .-(AACTL/CAPCTY) 

IF ( sue . L r . SOCL . AND . ACCEL . NE . 0 . ) S0CI.-30C 
C CHANGE CELL CURRENT TO ACCOUNT FOR CHANGE IN CELL RESISTANCE 
IEFF= [0ATTi!:(f1CSfND/Mt;) 

C CALC BATTERY VOLTAGE FROM DISCHARGE CURVE 

VBA r I =NC)FIN TEIK OISCHGjSOC , lEFF, 1 . , lEFFMX , . 8 1 , - . Oi , 1 0 U , -100 , 

6101,6) ) 

C DOF.S CALCULATED BATTERY VOLTAGE EQUAL BATT VOLTAGE ASSUMED AT START OF CALC 
lF<A3S(VBAf(-V8rOLl)).GT.VDELT> GOTO S? 

WOO 71 

C HAVL WE ITERATED THE MAX NO. OF TIMES 
5V IF(IirER.GE.rrERMX) GOTO 68 

C RESET AMP-HRS OUT OF BATT. AND INCREMENT HERAT ION COUNTER 
AITOU r=AHOIJT-( IBA F ftDELT/3600 , ) 
iITER=IITER+l 
GOTO 41 

C SET BATTERY VOLTAGES AND CURRENTS TO OLD VALUES 
60 VBA IT=VBT0IJ) 

IBATT=IBTOLD 

lEFF-IBAmiMCSfND/HC) 

C WAS THE EFFECTIVE BATTERY CURRENT ABOVE ITS MAX 
71 IFdEFF.LtMEFFMX) GOTO 70 

C COMPUTE BATTERY CURRENT CORRESPONDING TO MAX EFFECTIVE CURRENT 
I3ATl = IEfFMX?<MC/flCSTNO) 

C COMPUTE NEW BATTERY OUTPUT POWER 
PBOUT-IBArTSVBATr 

C ESTIMATE NEW INVERTER OUTPUT POWER 

P INV=P80U r-3 . OCIBATT-O . 035«(PBOUT-P ILOSS ) 

C ESTIMATE MOIOR OUTPUT POWER AND ELECTRIC POWER AT CVT INPUT 
PMrPCT=((PINV*0.V)/PMrRTO)»10O. 

PM0T0R==P1NV!|!INTEQ(EFFHTR,PMTPCT ,WMOTOR,200 . ,1466 . 

W. 5, -2. 5, 146. 61, 437. 82, 81,7) 

PMGPCT=(PMOTOR/PMGRTE)*iOO . 

PaEC=PMOriJR»IHTeQ<GrteFF,PMGPCT,WHGPCr, 100. , 100 , 

61. ,-l. ,10., -10. ,101,11) 
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PrLLeK=IFLLBK+l 
C CALC NEU HE POWER 
PHE=PCVr-PeLEC 

C IF BATT IS DISCHARGED SET 80% DOD FLAG 
70 IF(30C.LT.0.2) I30DOD=i 

C CALC. BATT. ENERGY USAGE FOR THIS INCREMENT 
B£NRGY=(PBOUr*DELT)/l.E+06 

C CALC BATT INTERNAL VOLTAGE AND ENERGY WEIGHTING FACTOR 

VINf=NC5|tINTEO(DI3CHG,SOC,0.,l.,IEFFHX,.01,-.01,100.,-i00.,i0i,6) 

I»SEFF=VINT/VBATT 

C CALC. TOTAL BATT ENERGY USAGE SINCE START OF TRIP 
£B TO r=E BTOT+BENRGY5D ISEFF 

C 

t he charging COHPUrATIONS 

C 

C IS HE CHARGING INHIBITED 
80 IFIINOCHG.EQ.I) GOTO lOI 
C IS BATT DISCHARGED 

1FISOC.lt. 0.2) GO TO 90 
C IS BATT PRESENTLY BEING CHARGED 
IF(IHECHG.NE.i) GOTO 101 
C IS BAff BELOW MAX CHARGING LEVEL 
IF(SOC.LT.SOCHAX) GOTO 90 
C STOP CHARGING 
I1IECHG=0 
GOTO 101 

C IS THIS HIGHWAY DRIVING 
90 IF(IHWY-l) 102,103,103 
C INHIBIT ALL FURTHER USE OF BATTERIES 
103 INOELE=l 
GOTO IBS 

C IS THIS THE CRUISE PHASE OF THE STC 
J02 IFITIME.LT.TCOAST. AND. TIME. GE.TACCEL) GOTO 105 
101 PCHG=0. 

GOTO 140 

C SET- BATT CHARGING FLAG 
105 IHECHG=1 

C IS POWER SPLIT BEING CHANGED 

IFIPELEC .GT. 0.) IFLLBK=IFLLBK H 
C SET HE POWER TO TOTAL 
PM£=P1IE+PELEC 
PBOUT=0. 

I‘LLEC=0 . 

PMUTOR=0 . 

PINV=0. 

C HAS ELECTRIC POWER BEEN CHARGED TO BATT THIS ITERATION 
IF (lELEC.NE.l) GOTO 106 
C RESET BATTERY QUANTITIES 

AHOIJI>-AHiJUf-(IBATT*O£LT/3600 . ) 
i BTOT=EBTOT-BENRGY»DISEFF 
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IBATT-0 , 

(: CALC CHARGING POWER 
106 PCI1G=PHEHAX“PHE 

IFCPCHG .GL. 0.) GO TO 100 
Pi:HG=« . 

GOTO 140 

C CALC POWER SUPPLIED TO ELECTRIC HOTOR 
100 «MGPCT= ( WCVl /HMOR TE ) $1 0 0 . 

PHGPCf=(PCHG/P MGR 104:100. 

PMOTOR=PCHG*1NTEQ(GREFF,PMGPCT,WMGPCT,100 . ,100 . , 

61. , -l ,10. ,-10. ,101,11) 

C IS THE CHARGING POWER GREATER THAN IHE MAX MOTOR POWER 
IF(PMOrOR.LE.PMrMAX) GOTO 110 
C SET POWER AT MAX AND CORRECT POWER CHARGED TO HE 
PHOrOR=PMfMAX 

l'HGPCT=(PMOTOR/PMGRTE>4:lOfl . 

IT.HG=PMTMAX/INTEQ(GREFF,PMGPCT,WMGPCT,iOO. ,100. , 

61. , -l. ,10. ,-10. ,101,11) 

C CALC POWER FROM GENERATOR TO INVERTER 
110 WMOrOR=WCVT/RGRM 

l'HTPCT=(PMOTOR/PMTRTD)*100 . 

PINV=PMOTOR«INfF.Q(EFFMrR,PMrPCT,WMOTOR, 200. ,1466. , 

62. 5, -2.5,146.61,439. 82, 81, 7) 

C IS INVERTER POWER TOO GREAT 

)i(PINV.LE.PINVMX) GOTO 120 
C SET ROWER AT MAX AND CORRECT POWER CHARGED TO HE 
PINV=PINVMX 

PMTPCT=<PINV/FMTRTD)*10fl. 

P CHG=P I NV/ ( IN 1 EO ( EFFHTR , PM TPC T , WMOTOR , 2 0 0 . , 1 466 . , 

62. 5, -2.5,146.61,439.82,81,7) 

6 ?IHT£Q(GRefF,PMGPCT,WMGPCT,100.,100.,l.,-l.,10.,-10.,101,ll)) 

C INITIALIZE ITERATION FLAG 
120 IICHG=0 

ETHi=PINV/P HOTOR 

C 

C HE CHARGING BATTERY LOOP 

C 

C CALCUi.ATE 8 ATT INPUT POWER AND CURRENT 
123 PBOUl =-P lNV+3 . 0»(P INV/VBCHG) 4 0 . 035*PINV 

IBArT=PBOUr/VeCHG 
PVBATT=VBCHG 

C IS CHARGING POWER TOO GREAT 

IF(PBHAXi|t(MC./MCSTND)+PPOUT.GE.G.) GOTO 121 
C SET CHARGING POWER AT MAX AND CAl.C NEW BATI CURRENT 
PB0Ur=^-PBHAX»(MC/MC3TND) 

IBATT=PBOUT/VBCHG 
C IS CHARGING CURRENT TOO LARGE 
521 IRIBMAXUBAn.GE.O.) GOTO 122 
C SET BATT CURRENT AT MAX AND RECALCULATE CHARGING POWER 
IBATT= -IDMAX 
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PSOiJT=i:BATr#.VBCHG 

C SAVt OLD VALUE OF BATTERY VOLTAGE 
12E V.On]LO=VliCMG 

C GALC. CHARGING AMP-HRS AND ACTUAL AHP-HR3 
AHCHG=AHCHG-aBATT*DELT/3A0Q . ) 

AAC I L“ AHOU f-AHCHG-AREG£N 

C CALC AVER BATl CURRENT FOR SIC AND HIGHWAY DRIVING 
IFdHUY.NE.l) GO 10 125 
IAVE=3600 .#AHOUT/(TIME+DELT) 

GOTO 130 

C CALC AVER BATT CURRENT SINCE START OF STC DRIVING 
I2S IAVE"< 3600 , *AHOU f > /( ( T IHE H>£L f ) +FLOA KNCYCLE-i ) * fSTUP > 

C CALCULATE STATE OF CHARGE AND BATTERY VOLTAGE DURING CHARGING 

130 CAPKG=PUKRr<IAVe/rtG) 

CAPCTY-CAPKGiitMC 

SOC= 1 . 0 -(AACrL/CAPCrY) 

lEFF=--IBAUit:(hCSTND/MC) 

VBCHG=NC*INTE0<VCHG,(AACTL/HC),IEFF,15. ,300 . ,1 . ,0 . ,25. ,0 . ,13.12) 
KIGEFF=NCY.INTEG|(DISCHG,SOC,0 . ,1 . ,IEFF«X , . 01 . 01 ,100 . ,-lOii . , 

6101,6) 

l/VBCHG 

C DUES CALC VOLTAGE EQUAL VOLTAGE ASSUMED 

)F(AB3(VBCHG-VI!T0LD).LE.VDELT) GO TO 133 
G HAVE G£ ITERATED MA.XIHUM HUMBER OF TIMES 
IFdlCHG.GE.lTERMX) GO TO 13A 

C INCREMENT rfERATlOH COUNfER AND RESET CHARGING AMP-HRS INTO BATT 
IICHG=IICHG+1 

AIICHG=AHCHG! (IBATKDELT/3600 . ) 

Gl! TO l2o 

C SET BATTERY VOLTAGE TO OLD VALUE 
1,54 VBCHG=VBTOL0 
C LIMIT SOC TO MAX 
533 IFISQC.GT.1.0) SOC=l. 

C STOP CHARGING IF BATTERY HAS REACHED MAX CHARGING LEVEL 

131 IldSOC.GE.SOCMAX) IH£CHG=0 

(: CALC TOTAL BATT ENERGY USED SO FAR ON THIS TRIP OR FROM START OF STC'S 
£B fiyr=EB fOTK (PBOU rUDtlT ) /I . E ^06) *REGEi-T 
C INCREASE REQ'D HE POWER TO INCLUDE CHARGING POWER 
PIE^PIIE+PCHG 

C IS HE POWER TOO LARGE 

140 ITIPHE.LE.PHEMAX) GOTO 160 
C IS BATT DISCHARGED 

IRSOC.GE.0.2) GOTO 141 

IOVRLD=IOVRLlHl 

PSIlUlT hPCVr-PHEMAX-PELEG 
IF(PSHORT) 136,136,137 
136 WRITE(7,i072) PSHORT 

1072 FORMAT!" AN ERROR HAS OCCURED IN COMPUTING THE", 

6' SUM OF THE SQUARED POWER OVERLOAD (P0VRLD=*,F9.2,*)*) 

GUfO 160 
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13? POVRLD=POVRLD+PSl(Oin*PSHOR1 
lillfO 160 

t HAVE UE ATTEHPTfD TO SUITCH MORE POWER TO El.E'CTRIC MOTOR 

141 IfdSWTCH.NL.i) GOTO 142 
lOViiLD-^lOVRLOH 
PSHORT-PCVT-PHEMAX-PELEC 
IF(PSHORT) 136,136,137 

0 3E! SWITCH FLAG 

142 ],SWTCH=1 

t; DID WE USE ELECTRIC POWER IH THIS INCREMENT 
IF(PELEC.LE.l).) GOTO ISO 
C RESET AMP HOURS OUT Of BATT 

AI IOUf=AHGUT- (IBArr<0aT/36fl() . ) 

C BTOT--EBTOT-BLNRGY<DISEFF 
C SHirCH EXCESS POWER TO ELECTRIC MOTOR 
ISO PHE^T'HEMAX 
PELEC=;'CVT-PilE 
11 LLBK=1ELLBK+1 
GOTO 30 

C CALC HE TOROUE, SPEED, FUEL FLOW RATE 
160 WllE=WCVl 
T((E=PHE/WHE 

F UELR T^INTNO ( HEMAP , THE , WHE , PHEMAX/20 0 . , 840 . 0 , SATORQ , AESPED ,6 , 6 ) 

•i ti'ilSCLE 

BSFC1=FUELR T/PHEX:4178722 . 1 
C CALC HE FUEL VOLUME USED SINCE START OF TRIP 
V0LT0T=V0LT0T4-(FUELRT)|tDELT) 

C [HCREMEtir TIM!,. 

nME=-TIME+DEL.T 

NT[ME=NriME-tt. 

Ct 

C HIGHWAY PRINTOUT 

I F ( 1 A . EQ . I A /S?#S7 . AND . R UN2 1 . EQ . 1 ) WR IT E ( B , 4 0 1 ) 
ia=j:a+i 

401 FORHATI'l'/NCYCLE TIME VAVE UDRIVE WGEAR WCVT 

■</WIIE THE WHO FOR lAVE lOUDUO IFLLBK lOVRLD 

6' ISWTCH IITER IHECHG IICHG IIREGN') 

CFIRUHEl ,E0.1) WfilfECdjSOUO) HCYCLE, TIME, VAVE, WDRIVE, WGEAR, WCVT, 
6WHE,THE,WM0T0R, 

6 1 AVF; , 180 DOD , If LLEiK , lOVRLD , ISWTCH , I ITER , IHECHG , IICHG , I IREGN 
•5000 FiJRMAra4,2X,F'8.0,lX,FiJ.2,lX,F8.3,lX,F'c3.2,lX,FB.3,lX,F8.3,lX, 
6FB.2,iX,F8.2,lX,F8.3,lX,0a8,iX)) 

TF ( I B . EO . r 8/S7ri7 . AND . RUM21 . EM . 1 ) WRITh*(V,402) 

402 F(JRMAT('1','NCYCLE TIME PAXlE PDR.IVE PGEAR PCVT ^ 
■V RCVl PHE PCHG PELEC PMOFUR PINV ', 

6'’ PBOUT IBATT VBATT SOC PROAD PREGEN 
6 •• 8SF C ' , 4X , ' FUELR T ' , 'SX / Ei-Ml ' ) 
ni=-IB+l 

II- ( R UN21 . EGi . 1 ) WR ITE (9 ,50 1 0 ) NCYCLE , T IME ,P AXLE , P DR I VE ,PGEAR , 
M'CVTjRCVr, PHE. PCHG, 
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iPELEC,PHOTOfi,PIMV,PBOUr,IBAVT,PVBAn,SOC,PROAD, PREGEN, 
J,6SFCi,FUELRT,EPHl 

SOifl FORMATa4,2X,F8.fl,lX,F8.i,lX,F8.i,iX,F8.1,iX,F8.i,iX,F0.3,lX, 
iF8.i,lX,F8.2,iX,F8.i,iX,F8.i,iX,F8.2,iX,F8.2,iX,F8.2,iX, 
iF8.2,iX,F8.4,iX,F8.2,lX,F8.?.,lX,F9.4,iX,Fifl.8,iX,F7.6) 

C$ 

C 18 THIS HIGHWAr DRIVING 
iFaHwr.EQ.i) Goro r/o 

C HAS THE 'COAST' PHASE OF THE STC BEEN REACHED 
tF< riHf-TCOAST) 10,130,180 
C HAS HIGHWAY DRIVING BEEN COMPLETED FOR THIS TRIP 
170 CONTINUE 
C 

If (TIME-TRIPTM) 20,300,300 

C- 

C COASf PHASE 



C SET HEAT ENiilNE CHARGING POWER TO ZERO 
180 PCHG=0. 
nrER=« 

) ICHG=0 

pyBATT=NC«INTEQ(DISCHG,SOC,0.,l.,IEFFMX,.fll,-.01,100.,-lfl0., 

f.101,6) 

C CALC CURRENT VELOCITY 
VELO=V(NnME) 

C CALC VEHICLE DRAG FORCE 
190 F0=CDA»RHl]»VEL0*VtL0/2 . 

C CALC ROLLING RESISTANCE FORCE. 

f R = M nicl ;R AV Y. < 1 0 . + 0 . 0 1 *VELlH-0 . 0 0 0 08» VE LO» VELU ) / 1 0 0 0 . 

C CALC GRADE FORCE 
8eTA=B£TafriME) 
fGR=HT*GRAV*SlN(BETA) 

C CALC VEHICLE VELOCITY AT END OF TIME INCREMENT 
VN£XT=VELO- ( DELTjKFGR+FD+FR )/HT ) 

C CALC AVERAGE CURRENT AND VELOCI TY 

TAVE- <3600. MAHOUT ) /( T IME+FLOAT(NCYCLE-l ) JdTSTOP ) 
VAVE=.S*(VELOfVNEXT) 

C ZERO ALL DRIVE TRAIN POWERS, ETC 
PAXLE=U. 

PDR1VE=0. 

PGEAR=0, 

PCVl-0. 

f’ELEi>l) 

PHE=^0. 

pMoros=o, 

fINV=0. 

PILOSS-U . 

JBATT=0. 

PBOUT-l) . 

IME=0. 
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C CALC DRIVE SHAFT SPEEDS 
UHE=0 . 

«DRIVE=VAVE/(RTIRE)(tRDIFF) 

WGEAR^UDRCVE/RGEAR 
WCVT-WGEAR/RCVT 
«Hl]rOR=HCVT/RG«H 
C INCREMENT TIME 

TIME=TIME+DELr 

HTIHE=NTIME4i 

FuaRr=8. 

C4 

I F a A . EQ . I A/S7*5? . AMD . R UN2 i . EQ . i ) HR I TE ( 8 > 4 0 i ) 

IA=IAa 

IF(RUN2i.EQ.i) HRITE(8,5000) NCYCLE,TIME,VAVE,UDRIVE,«GEAR, 
iyi:vr,HHE,rHE,uHaTOR, 

4 .1 AVE , 180 DOD , IFLLBK , lOVRLD , ISHTCH , 1 1 TER , I HECHG , IICHG , I IREGN 
IF< 1 8 . EQ . IB/S7*S-; . AMO . RUN21 .Eq . 1 ) WRITE< V , 402 ) 
ni=-IB+i 

I RR UN2 i . EtI . i ) HR ITE (9 , 50 i 0 ) NCYCLE , T IME , P AXLE , P DR I VE , PGEAR , 
ii'cvr,Rcvr,pi!E,PcnG, 

&PELEC,PMOTOR,PIMV,PBOUT,IHATT,PVBATT,SOC,PROAD,PREGEN, 

iBGFCijFIJeLKT,EF«l 

c-t 

C SAVE VELOCITY 

V‘;NrrME)=VNExr 

C HAS THE BRAKING PART OF THE STC BEEN REACHED 
IFd'IME.GE.TBRAKE) GOrO 195 
C INITIALIZE VELOCTITY FOR NEXT TIME INCREMENT 
VELO^VNEXT 
GOTO 190 

^ 

C 8RAKING PHASE - 

C 

C CALC VEL DECREASE PER FIME INCREMENT DURING 3REAKING 
195 DELV= ( VNEXT*DELT ) /( TSTOP-TBR AKE ) 

C CALC THE AVERAGE VELOCITY OVER THE HEX! TIME INCREHEMT’ 

200 V1=V(NTIME) 

V2=V1-DELV 

VAVE=(Vl+V2)/2. 

C CALC POWER LOST DUE TO ROAD LOSSES FOR THIS INCREMENT 
FO=CDA!f:RHO*VAVE*VAV£/2 . 

I R=MT*GRAV«( 10 . +0 . 01*VAVE+C . 00008IVAV£1iVAVE)/l 000 . 
BErA=B£r<NriME> 

FGR=HT1KGRAV*SIN(BETA) 

PROA0=(FGR+FR^FD)IVAVE 

C CALC THE POWER DUE TO VEHICLE ENERGY CHANGE 
PWRSTP=(Mr»l.l/2.)»((Vl*Vl)“(V2IV2))/DELr 
C CALC POWER AVAILABLE FOR REGENERATION 
PREGeN=PWR3TP-PROAO 
C IS REGEN POWER ABOVE MIN 
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lF(PKEG’c«-8£i;mN) 231), 2,50, i!;5i 
C StT REGENERATED POWERS 10 ZERO 

230 PAXU>1) 

PDRIVE=0 

PGEAR=0. 

PCVT=(I . 

PHlJTOK=0 

PINV=fl. 

P [1.033=0. 

PB(JI!T=0. 

PHE=0 . 

UiATT=0. 

PEl£C=(). 

GOTO 235 

I.: 13 BAT TER Y FULLY CHARGE.)) 

231 IF(SOC.GE.l.O) GOTO 23(1 
C CAU; AXLE SPEED AND TORQUE 

WAXLE=-yAyE/fiTIRF 

rAXLE=PREGE:N/WAXL.E 

PAX1E=PREGEN 

C OAl.C DR WE SHAR TORQUE, SPEED, POWER 
WDRi:Vr>UAXl..LVKDi'r>‘ 

Pl)PCT=(PREGLN/PDaATE)liOO . 

WDPCT= ( WAXLE/WOORTE )#10 0 , 

fI)KiyE=PREGEN*lNTEQ(GRErF-.PDPCT,WDPCT ,100 . ,100 . , 
?>l.,-l.,lfl. ,-10.,i0i,l.l) 
rDRIVE=PDRWE/WDRW£ 

C CALC I1UANTITIE3 AT OIlTf'L SIDE Oi' COT 
WGEAR=WDRIVt/RGEAR 
PGPi;T=a'DRi;OE/PGRAT'£)«iOO 
UGf'CT=(WDR WC/WGORTE)*100 . 

P GEAR --P OR 1 VE* IN TEtUGREFF , PGPC f , WGP C f , 1 0 0 . , i 0 0 . , 
i'j.,-1., 10. ,-10. ,101,11) 

rGEAR=PGEA«/WGEAP 
('; CALC CVT RATIO 
PHt.=0 . 

GALL RATIO 

C CALC QUANTITIES ON ENGINE SIDE OF CVT 
ui:vr=wcEAK/RCvr 

P G VT= ( P G t AR ) X (1 NTNO ( E FF , ( P GEAR /1 0 0 0 . ) , R CV T , 1 0 0 . , 3 . 4 , 
f.APCVr,ARCVT,e,7)) 

TCVT=PCVT7WGVT 
C LIMIT CVr POWER fU MAX 

1F(PCVT.G( .PCVTMX) PCVr=PCVlHX 
PELEC=PCVf 

C CALC POWER AT MOTOR SIDE OF MOTOR GEAR REDUCER 
WMOflJR=UCVT/RGRH 
Pt!(;PCT=(PCVT/PMGRTE)*100. 

WMGP CT= ( WCV r/WMOR TE ) *10 0 . 
|■|■^(lTOR=PCVr*lNTE£!(GREFF,PMG?GT,WMGPCT,100.,100., 
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ii., 

1H0r0fi=PH0T(JR/Wh'0T0R 

C If- RLGLN£RftT10H P(i«£;R IS GREATEIR THAN MAX MOTOR POMtR SET POWER AT MAX 
IK(PMOrUR .Of .PMTMAX) PMOTiJR=PHTMAX 
C CALC POWER FROM GENERATOR TO INVERTER 
PM rPCi -• ( PMO rOR/PH TRTD ) »iC 1) . 

PlNV=PMOTORt 

ilNTEQ(EFEMIR, PhTPtT>WMOTOR, 200., 1466. , a, S, -2. 146. 61,439. 02,81,7) 

C INiriAi.IZE ITERATION COUNTER 
IIREGN=0 

C IF INVERTER POWER IS TOO GREAT SET POWER AT MAX 
IF(PINV .CT . PINVIiX) PINV=PINVhX 
liMl=PINV/PMOTOR 

c CAu: PArfERr input power and curreht 


c REGENERAriVE BATTERr VOLTAGE LOOP 

C 

24S PBOOf-- PINV•^3 .«*(PINV/V8CHG)-!-0 . «35ifP IHV 
niATT=Pf)OUT/VBCHG 
PV3ATT=VBCHG 

C IS CHARGING POWER TOO LARGE 

IF(PBMAX^(MC/MCSrND)-«-PBOUr.GE.U.) GOTO 23V 
C SET POWER AT MAX AMD RECALCULATE BATT CURRENT 
? GOU r=-P BMAXf. ( HC/MC3 CND ) 

IHATT^PBOUT/VDCHG 
C IS CHARGING CURRENT TOO LARGE 
239 irdBMAX+IBATT.GE.O.) GOTO £41 
C SET CURRENT AT MAX AND RECALCULATE CHARGING POWER 
1BATT--IBMAX 
PBOIjr=ieATI»VflCliG 
C SAVt OLD VALUE Of BATT VOLTAGE 
241 VBTGLD=VBCHG 

C CALC BATT AMP-HRS OUT AND AVER BATT CURRENT SINCE START OF STC'S 
AREGEN ^AREGEN- ( IBAT f*OEL T/360 « ) 

AACTL=AHOUT - AREGEN-AHCHG 
IAVE=(360fl .?AHOUT)/< riMF>FLOAT<NCrCI.E-i)TTSTOP ) 

C CALC STATE OF CHARGE AND BATTERY VOLTAGE DURING REGENERATION 
CAPC rY=MC»PUKRTa AVE/MC ) 

SOC=I.fl-(AACTL/CAPCTY) 

I£Ffr-IBArri!(MCSTND/MC) 

VBCHG=NC«INrEO(VCHG,(AACTL/MC) ,IEFF,IS. ,300 . ,1 . ,0 . ,25. ,0 . ,13,12) 

R EGEFF=IICY IN TGfl ( 0 ISCHG .. SCC , 0 . , 1 . , lEFF MX , . 0 1 . - . 0 1 . 10 )) . , -1 0 0 . , 

6101,6) 

6/VBCHG 

C DOES CALC BATIERY VOLTAGE EQUAL ASSUMED VOLTAGE 
IFIABSIVBCHG-VBTOLD) .LE. VDELT) GOTO 243 
C HAVE WE ITERATED THE MAX NO OF TIMES 
IFCIIREGN.GE.ITERMX) GOTO 244 
C INCREMENT IIERATION COUNTER AND RESET A-H INTO BATT 
IIREGN=IIR£GN+i 
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AREGt:N=AREGEN+aBATT*DELT/3600 . ) 

GOTO 2« 

C GET BATT VOLTAGE TO OLD VALUE 
2A4 VBCtlG=VBTOLD 
C LIMIT SOC TO MAX 
243 IFTSOG.LE.i.O) GOTO 242 
SOC=i.O 

C CALC BATTERY CURREMT REQUIRED TO FULLY CHARGE BATTERY 
AR£GEN=AREGEM+(IBATT*DELT/360l) . ) 

If ILL=360 0 . /DELT*( AHOUT-AREGEM-AHCHG ) 

IBArT=-IFILL 

AREGEH=AREGEN-aBATT*DELT/36flO . ) 

A ACTL= AHOU l-AREGEN- AHCHG 
I LI F=-IBATT*(MCSTND/MC) 

VBCHG=HC*IKTEfJ(VCHG, (AACTL/HC) , lEFF , IS . ,300 . , i . , 0 . ,25 . , 0 . , 13,12) 
REGEFF=HC*INTEU(DISCHG,SOC,fl.,i.,IEFFMX,.01,-.01,10U.,-\00., 
Gi01,6) 
i/VBCHG 

PBOUT=VBCHG«IBATT 

PVBATT=VBCHG 

C SET TOTAL ENERGY OUT OF BATTERY TO ZERO 
£BfOr=0. 

GO TO 235 
242 CONTINUE 

C CALC TOTAL BATT ENERGY USED ON THIS TRIP OR FROM START OF STC'S 
EB rOT=EBTO r+ < ( P 80UT*DELT ) /I . E+06 ) IREGEFF 
C INCREMENT TIME 

235 TIM'r>TiME+D£LT 
NTIME=NTIME+1 

Cl 

IF ( I A . EQ . I A/S7457 . AND . RUN21 . EQ . 1 ) WR ITE(8 , 40 1 ) 

IA=IA+1 

IF(RUN21.£Q.l) WRITE(8,5fl00) NCYCLE,TIME,VAV£,HDRIVE,WGEAR, 
iHCVT,yHl.:,THE,WMOTOR, 

i I AVE , 180 DOD , IFLLBK , lOVRLD , ISUTCH , I ITER , IHECHG , I ICHG , I IREGN 
lEdB . E(J. 18/57*57 . AND . RUN21 .ElU ) UR ITE(9 , 402 ) 

3B=iB+l 

IF(RUN21.EQ.l) «RITE(9,50i0) NCYCLE,TIMF:,PAXLE,PDRIVE,PGEAR, 
<.PCVf,RCVr,PHE,PCHG, 

iPELEC,PMQTOR,PINV,PBOUr,IBATT,PVBATT, SOC, PROAD, PREGEN, 
f.B3FCl,FUELRT,EFMl 
Cl 

C HAVE HE REACHED THE STOP PHASE OF THE STC 
rFTTIME.GE.TSrnP) GOTO 236 
V(NTIME)=V2 
GO TO 200 

C CALC AVER CURRENT FOR ALL STC'S 

236 IAVE-(AHOUT*36flO . )/(FLOAT(NCYCLE)*fSTOP ) 

C AT END OF THIS STC, SAVE: BATT SOC, AVERAGE CURRENT, PETROLEUM VOLUME, 
C ELECTRICAL ENERGY, AMP-HR3 OUT, AMP-HRS REGEN AND BATT. STATE- 
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C ALL SINCL START OF STC'S 
IDOD(NCya£)=I30DOD 
IVFSTC(NCYCLF.)==IAVE 
HOUl(NCYCLD=AHOUT 
RE(iEN(NCYCLE)=AREGEN 
£BTDTL(N(:YCLE)=EBTOr 
FET3TC(MCYCLE)=V0LT0T 
G0 CSTC(HCYCLl')=30C 
CHAKGE(NCYCLE)=AHCHG 
C# 

I6=-I6+l 

ira6/NCYLHX*NCYLHX . EQ . 16 . ANB . RUN21 . LE . i ) WR ITE(6 , 404 ) 

4(1 ,t6 rOTLa ) ,HOUT( I ) , I DODd ) , IVF.SrC( I ) ,PErsrC( C > ,REGEH(r) , 
4SUCSTC (I ) ,CHARGE( I ) , 1=1 ,NCYCLE) 

404 FlJRHAr<' ',/lX/ I*' ,SX, 'EBTOrL(I) ' ,4X, 'HOUm')' ,4X, 'IDOO(I)' , 
44>:/IV£STa)',3X,'PETSTCa)',3X,'REGEN(I)',2X,'S0CSTC(I)',3X, 
4'CHARGE(I)',/36a2,3X,2X,F9.S,2X,F9.5,3X,18,3X,Fy.S, 
64(2X,F?.5)/)) 

C* 

C IS THIS THE END OF THE STC WHICH PRECEEDS ALL HIGHWAY DRIVING? 

IF (NCYCLE-NSTCFR) 238,237,238 
C SAVE PRESENT VALUE OF BATTERY VOLTAGE 
237 VHt)BTT=VBATT 

C INITIALIZE POWER REGENERATION VARIABLES 
238 PREGEN=0 . 

PROAD=0. 

I[REGN=« 

C HAVE WE GONE THE HAX NO OF STC'S 
IF ( NC YCLE-NC YLHX ) 240 , 2S0 , 2S0 
C INCREMENT CYCLE COUNTER AND INITIALIZE TIME 
240 NCYCLE=NCYCLE+1 
TIME=fl. 

NfIHE=i 
GOTO 10 

C INITIALIZE TRIP CONSUMPTION VARIABLES 
2S0 PETROl-0. 

1HWY=0 
£N!;LEC=0 . 

NTR1P=1 

C DETERMINE LENGTH OF CURRENT TRIP 
260 LENGTH=rSIPLH(Nf«IP) 

C IS TRIP ALL STC'S 

IF(LENGTH.GT.ALLSTC) GOTO 290 
C CALC NUMBER OF STC'S IN TRIP 
STCNQ=L.ENGTH/STCD1S 
NOSrC=STCNO 

I)EI.STC=STCNi)-FLOAT(NOSTC) 

C CALC AMP -HRS FOR fHIS TRIP 

AACTLl=HOyT(NOSTC>-REGEN(NOSTC)-CHARGE(NOSTC) 

AACTL2=H0Ur(N0STC^•^)“REGEN(N0STC+l)-CHARGE(N0STC^•1.) 
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fiACTL=AACTLi-KAACTL?.-AACTU)!«:DeLSTC 
C CALC PETROL CONSUMPTION FOR THIS TRIP 

PETSOL=l’LTSTC(NOSTC)+DCLSTC*<PE1STC(NOS'lC+i)-PETSTC(NOSTC>) 

APETT(NTRIP)=PErROL 

t CALC ELEC ENERGY CONSUMPTION FOR THIS TRIP USING POWER AND USING AMP-HRS 
FMLEC-EB rOTUNOSTC>+D£LSVCa(E0TOTL(NOSI'C+i )-EBTOTL(NOSTC) ) 
AI:lETP(NTfiIP)=EN£LEC 

ENRGy=NCIAACTL*3.6E-03*SUINV(AACV,AACTL/MC,ACR >8) 

AELETA(HIRIP)=F.NRGY 

C-* 

li 

V, DETERMINE NO OF SUCH TRIPS/YEAR 
T«IPYR=TRIPNO(NfRIP) 

C DETERMINE YEARLY PETROL AND ELECT CONSUMPTION FOR THIS TRIP 
PLiYR(NIRCP>=rRtPYfi*?ETROL 
1 1 1 CYR (NIRIP >=fRIPYR»ENEL£C 
£L£C(NTRIP)=TR.(PYRi!:fc.NRGY 

C WAS THE BATT DISCHARGED SOMETIME DURING THIS TRIP 
IFdOUDCNOSi'C) NE.l) GOTO £6S 
DlSFAC=-i 
GOTO 267 

21, S DISFAC=(i .-SOCSTC(NOSTC) I/O .8 

U CALC NO OF DISCHARGES PER YEAR FOR THIS IRIP TYPh 
267 D ISCH (N I R I P ) =DISFAC* VR IP YR 

C IS THIS mi: LAST TRIP TYPE 

270 IRNTRIP.GE.NIRIPX) GOTO 275 
HTHIP=NIRTPi-1. 

CO TO 260 

C DLTERMINE WALL PLUG ENERGY REQ'D TO MEET ELECT . ENERGY CONSUMPTION 

271 00 201) .NTRIPX 
ELEC(I)=ELEC(I)/CHGEFF 
AELE r A T I ) = AELET A ( I > /CHGEP F 
AtLElPm^ALLETPID/EFFBC 

28(1 i;i.ECYR 1 1 )=ELECYfi (I) /EFFOC 
IF(R0N2i.LL.i) WRITE(6,I26) 

126 FORMA I <1X , ' T ,7X , ' OISCHI I ) ' ,SX , 'PETYR ( I) ' , OX / ELECYR (I ) ' , 

6SX/ELEC(lV.SX,'TRIFN0a)',3X/TRIPLNa)',lflX, 

63X, 'PETROL' ,6X, 'ENELEC' ,6X, 'ENRGY' ) 

DU 6S79 I"i,8 

6S77 IF(ftUN21.LE.I) WlUTE(6,408) I,DlSCH<I),PETYRa),ELECYRa), 

^i;lec< I ) , tr ipnoi i > , ir iplnt i > , apett ( n , aei.e n* ( i ) , aelet at i > 

AOO FGRHAT(iX,I2,SX,F9.A,2X,F9.A,3X,F9.4,6X,P9.A,3X,F6.0,AX,P9.0,lCX, 
6,lX,F10.4,2X,FiD.4,2X,,Fi().A) 

C-i +++!++++++++++++++++++++••+++++■•+++++<++++++ ++■'■ 

SOMA -0 . 

EAKHY=0 , 

SUMP=0. 

P(..TVY=0. 

PETV=0. 

DO 83341 1=1,8 
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t;('iKwiiR(C)-.:Lfi:ca)/3.6 

Al'KHA(l)=Af.LtTAa)/3.6 

A£Klifn.[)=AELi:rPa)/3.6 

Pt1V=l'ETV+AP£TT(J.) 

pf..ivY=?Em+PEVYaa) 

SUdA-SUrtA+AEKHAd) 

'3UilP=SUilP+AEKIlPa) 

83341 EAKHY=EAKHy+fAKWHR(I) 

£F(RUK£!1.N£.4) GO 10 i2Z 
)F(CX(1).NL,0.) GOTO 33244 
aRirE(7,3346) 

)((RUN21.tO,4) WRITE(8,3346) 

I F ( H UN2 1 , to . 4 ) aa I T£ ( V , 3346 ) 

WRITE (1(1, 3346) 

3346 FORMA! (2X,'PS%lflX/l',liX/2MlX/3‘',llX/4SliX/S',ilX, 
t'6',ltX, '7', llX/O',10X/3Urt', /'«#«') 

33244 CONTINUE 

WU1T£(7,3432) CX(l) ,(AEKHA(1) ,1=1,0) ,SUMA 
M',!I TE ( 0 , .3 432 ) CX ( 1 > , ( AP E i f ( I ) , £ =1 , 8 ) , P F. r V 
WRITE(9,3432) CX(l) jlAEKHPd) ,1=1,8) , SUMP 
HR n E( 1 R , 3432) CX( 1 ) , ( EAKWURd ) , [=1 ,8 ) ,CAXHY 
3432 FORMAT ( 1 X , F4 . 2 , 1X,3<F11 . 5, IX, > ,F11 . S ) 

323 CONTINUf 

C»TH+4++4+++++++++++l+++++t++++++++H ++++++++++++++++ 

lF(RUN2i . LE . 1 ) HR CfE(6,3431 ) £OVRLO,POVRLO,Pr; iVY,EAKHY. INllCHG 
3431 r0RMAT(d0VRLD=',IS,2X,'P00Sl.li=',G13.7,2X,d0IAL PETKOLFJUM 
6'VllLUME PER YEAR (PE(yY)='' ,F?.2,' LITERS/ YEAR ',5X, 

6'TQTAl. ELECiRlC ENERGY CONSUMPTION PER YEAR (EAKHY)=' ,F9.2, 
h' KIUIWATI iiOUR-3/YEA«','jX,dNOCFiG=',Il) 
i^ETV10=PEf 01*10. 

RE TURN 

C CALC TIME IT TAKES TO MAKE THIS TRIP 

290 IR£PfM=FLiJAr(NSrCFR)*TSfi)PK(LENGrH-FLOAr(NSrCFR)*SrCDr>)/ilHYyEL) 

C IS THIS THE FIRS! TRIP REQUIRING HIGHWAY DRIVING 
IFdHUY.NE.O) GiJfO 20 

C INITIALIZE BATT MODEL AND FUEL CONSUMPTIONS 

C veAlT-OHUBn 

AH0UT=H0(mNSTCFR) 

l!JODaD=Il)UlHNSrCFR) 

)AVE=IV£STC(NSTCFR) 

V0Lf0f=PETSTC(N3rCFR) 

LBT01=EBT0TL(HSTCFR) 

AREGEN=R£GEN(NSrCFR) 

AHCHG=CHAKGE(NSTCFR) 

C SET VELOCITY, ACCELERATION AND GRADE FOR HIGHWAY DRIVING 
VAVE=HHYVEL 
ACCEL =0. 

BEfA=0. 

C INIIIALIZE TIME TO START OF HIGHWAY DRIVING 
r £ ME---FLU AT ( HS IFF R ) » TS TUP 
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c initialize; battery model quantities for start of highway driving 

CAPCTY=HC*PUKRT(IAVli;/HC) 

AACTL=AHOUT-AREGEN-AHCHG 

S(JC=i.O-<AACfL/CAPCTY> 

C SET HIGHWAY FLAG TO HIGHWAY DRIVING CONDITION 
(HWY=i 

C SET HWY DRIVING TIME INCREMENT 
l)EIJ=Dl-:LfHW 
GOTO 20 

C CALC YEARLY PETROL CONSUMPTION FOR THIS TRIP TYPE 
300 fRIPTR=IRIPNO(NT«IP) 

PETYR(NTRIP)^TRIPYR*V0LT01 

pi:r:iOL=voLfUf 

APETHNTRIP)=PETROL 

C CALC YEARLY ELEC ENERGY CONSUMPTION FOR THIS TRIP TYPE 
i:i.h‘OYR(NrRip)=iiUPYR*Eerof 
ENtLEC=tBT01 
AaEllMM[RlP>^--tHtLhC 

Lt.EC(NrRlP) = IRlPYR*NC*AACTL«3.6E-03TcSVINT(AACV,AACTL/MC,ACR/d) 

LNRGY=hL£C(NIRLP)/iSIPYR 

A(:LOA(NTiUP)“LNfiGY 

C SET DISCHG FACTOR TO TULL DISCHG IF DATT WAS DISCHARGED DURING (HIS TRIP 
IFasaDUO.N'E.T) GOTO 3iC 
D3.SFAC--i . 

GOTO 320 

C CALC DISCHG FACTOR 
.510 DISFAC=(1 . -SllC)/0 .0 

C CALC GO OF YEARLY DISCHARGES FOR THIS TRIP TYPE 
320 DlSCSUH IKTP )---Oi'SFAC*fR I'PYR 
GOTO 270 
LND 
I ND 
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•jUBROUirNE RATIO 

ctM%untnutunuurmm%tununtnuuunintuunnuuunnnu* 

i; 

C THIS DETtRHINES RATIO OF THE CONTINUOUSLY VARIABLE TRANSNISSIOM BA<3ED ON 
C THE LAST VALUE OF HIE GEAR POWER. AT VEHICLE START UP TIME THE VALUE 
C IS SET TO A MINIMUM OF .3. THE RATIO VALUE IS LIMITED TO PREVENT 
C EXCESSIVE MOTOR SPEED. THE MAXIMUM VALUE IS 3.4. 

C 

cumumumuuuttuuutunnnmtnunmmnnunnnuunntu 

REAL HVRSPOTS.' ,APRP<S) >AWMT<6) ..STALL(6) 

COMMQN/CC/CDA,RDIFF,RGEAR,RGRM,KTIRE 

COMMOH/RR A fiJ/N T [ME,UGEAR .RCV T , PMOTOR . PH TRTD , PGEAR 

COMMON/RSP LT/ ACCEL , SQC , PCVT , CX (IS ) , VAVE , P HEMAX , P HE , PELEC 

DATA RCVTMX/3.4//RCVTMN/U.3/,HMTRMX/i67S.S/, DELR/tl.i/ 

DATA BSFC/4.2/, WHEMN/i04.72/ 

DATA MTRSP D/S06 . 43i , 077 . 646 ,1172 . 86 , 1 31 7 . 47 , i 466 . 03.F 
DATA AWMT/fl . ,586 . 43i ,879 .646 ,ii?2 .86,1319.47,1466.08/ 

DATA 3 TALL/fl . , 50 . , 75 . 84 , 109 . 57 ,129.77,200./ 

DATA APRP/19.4,36.7,64.3,75.5,101./ 

C 13 VEHICLE JUST STARTING UP 
IF(NTIME.NE.l) GO TO 10 
C SET RATIO AT MINIMUM 
RCVT=RCVTMN 
RETURN 

C HAS HE POWER USED IN THE LAST ITERATION 
10 IFTPHE.GT.O.) GO TO 20 
C PICK RATIO TO MAXIMIZE MOTOR EFF 
PMT=PGEAR/(0.73*0.98) 

C COMPUTE MOTOR POWER IN PERCENT OF RATED POWER 
PliTPCT=(PHT/PMTRTD)*100. 

RCVT=WGEAR/(SVINT(MTRSPD,PHTPCT,APRP,5)*RGRM> 
i; Lltlir RATIO SO MOTOR SPEED IS NOT EXCESSIVE 
WMT=WGEAR/(RCVT*RGRH) 

IF(WMr.LE.UMTRMX) GO TO 40 
WHT=WhTRHX 

RCVT=WGEAR/(WHT!j:RGRM) 

GO TO 40 

C ESTIMATE A HE OUTPUT POWER FOR THIS ITERATION 
20 PH£ST-(PHE/(PHEH'ELEC))»(PG£AR/0.93) 

C CALC MOST EFFICIENT HE SPEED 

WHEFF- 8 SFC»(PHESr/PHEHAX)# 100 . 

J F(WHEFF . LT . WHEMN) WH£FF=WHEMN 
C CAL RATIO 

RCVT=WGEAR/WHEFF 
C CAL MOTOR SPEED 

WMT=WGEAR/(RCVT*RGRM) 

C CALC REQUIRED MOTOR POWER AT THIS SPEED 
PMT=(PGEAR/0.93)--PHEST 
PMTPCT=(PMT/PMriITO)#10». 

C CAN MOTOR DELIVER THIS POWER AT THIS SPEED 
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JO i:K<^’«n’c^.L^.svI^^r(;^lAa,WKr,A«m^ 6 )) gu to 40 

(,: DLCREASE RATIO TO INCREASE MOTOR SPEED 
RCVr'-RWr-DEi.R 
yHT=(JGEAR/(RCyf*RGRH) 

ii'(Rcof.Gr,ac‘*riiN) go ro jo 

C IS RATIO AliOVE. UPPER LIMIT 
40 ci-'<Ri:vr.Gf.RCviMX) rtcvi'=RCvr(ix 
IFTRCVT.LT.ROVTMN) RCVT=RCVVMN 
RETURN 
t.Hl) 

'END 
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sumoutine: spi. n 


FTK4,L 


zt 

i; 

c 

c 

c 


THIS DETERMINES THE POHER SPLIT BETWEEN THE HEAT ENSINE AND THE 
hLECIRIU MOTOR. IT IS A EUICIION OF THE ACCELERATION AND THE STATE 
or CHARGE OF THE BATTERY. 


zmmmmuuutmnutmttmmmummmunmmnmnmnm 

COMMON/R SP L T/ACCEL , SOC , PCV T , CX (IS ) , VAVE , PI lEM AX , P HE , PELEC 
C PS ^0 WHEN ONLY HEAT ENGINE USED 
C PS ■■i WHEN ONLY BATTERIES USED 
DATA VMINHE/1.34/ 

C IS VEHICLE VELOCITY BELOW KIN 
IFTVAVE-VMINHD 10,10,20 
10 l’S'5. .0 

GOTO 30 

20 PS=CX(1)+CX(2)SACCEL+CX(3.>*G0C 
TFd'S.GT.i. ) ps=l. 

)F(FS.LT.0.) PS=0. 

30 P£LEC^=PCVT)f.PS 
PHF.=PCVT<(l.fl-PS) 

RFTORN 
! ND 
END 
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FTN4;L 

SUBROUTINE SPLIT 

cttttummtmmtmmmutmmtnmmtmmmummtmnmm 

i; 

C THIS I'ETEKnINES THE POWER SPLIT BETWEEN THE HEAT ENGINE AND THE 
C ELECTRIC MOTOR. If IS A FUNCTION OF THE ACCELERATION AND THE STATE 
C or CHARGE OF THE BATTERY. (AIRESEARCH POWER SPLIT) 

C 

mumnmmmtmnumntmmmmmnmmmtmmmmttmt 

C0HMi]N/RSPLr/ACCEL,30C,PCVr,CX(lS),VAVE,PHEMAX, PHE,PELEC 
C PS =0 WHEN ONLY HEAT ENGINE USED 
C PS -I WHEN ONLY BATTERIES USED 
DATA VMINHE/i.34/ 

C IS VEIIICI.E VELOCITY BELOW MIN 
IRVAVE-VHINHE) 10,10,2(1 
10 PS^^l.O 
GOTO 80 

20 IFIACCEL.NE.O.) GOTO 40 
IRSOC.LE. .2) GOTO 70 
C IS THIS A HILL CLIMB 

IF(PCVr.GE.20flfl0.) GOTO 30 
PS=1 . 

Goro ilO 
30 PS=.3 

GOTO 80 

C IS REQUIRED POWER TOO LARGE 
40 IFIPCVT.GE. .8$PH£MAX) GOTO SO 
PS=.3 
GOTO 80 
SO PHE=PCVT 

IFIPHE.LE.PHEMAX) GOTO 60 

PHE=PHEMAX 

PELEC=PCV1-PHEMAX 

PS=PELEC/PCVT 

RETURN 

60 PHE=PCVr 
PELEC=0. 

p:t=o . 

RETURN 
70 PS-=i] . 

BO PELEC=PCVr»PS 

PHE=PCVT*(1.0-PS) 

RETURN 

END 

(.'ND 
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(^L'AL FUNCTION PUKRT(IAVE) 

C 

C CALCULATES CAfACITY PER KILOGRAfT OF BATTERY BASED ON THE AVERAGE 
C CURRENT PER KILOGRAH OUT OF fHE BATTERY. IT IS LIMITED TO A 
C MINIMUM OF O.S. 

C 

£,uitu%tnuuiuuunutt%u%u*%nnnnu*uu*uui*u%***t%uutuutti 

REAL lAVE 

)r(IAVE.LT,i0.173) GOTO iO 
plJKRT=2a.;588il - 3.7%?48*AL0G<IAVE) 

.UTPUKRT.LT. .5) PUKRT=.S 
RETURN 

1 0 PUKRT=i6 . 98831- . 3806714*1 AVE 

RETURN 
END 


FUNCTION SVINT(ARRAY,X,XARRAY,IDIM) 

C***************************************************************************** 

c 

C SINGLE VARIABLE INTERPOLATION SUBROUTINE. SPACING BETWEEN DATA 
C POINTS OF THE INDEPENDENT VARIABLE NEED NOT BE EQUAL. 

C 

C***************************************************************************** 
DIMENSION ARRAYaOIM),XARRAYO;DIM) 

IDMi=IDIM-i 
DO S IMAX--2,IDMi 
)r(X.LE.XARRAY(IMAX)> GO TO 7 
S CONTINUE 
7 IMIN=IMAX-i 

SVINT=((X-XARRAY(IMIN))*ARRAY<IMAX)+ 

4(XARRAY(IHAX)-X)*ARRAY(IMIN))/(XARRAY(IMAX)-XARRAYaMIN)) 

RETURN 

IND 

END 
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l;rAL FUNCTION INTNO(ARRAY,>;K.YY,XTOf>YTQP,)(AHRAY,YARRAY,IDIM,JniH) 

atm%n%m%%n%u%%*-mmmtmnm%nt%'tmnn%i%tnmntnmunnt 

c 

C UOUBLt VARIABLF INTERPOLATION SUBROUTINE. BOTH VARIABLES NEED KOI 
C HAVE EIJUAL SI'ACIHO BEIWEEN DATA POINTS. 

C 

znnut*t*n%utntx%un%ununntntnntn*uu%t*%ttu%utt%tnuntui 
C THCS LINEAR INTERPOLATION MUST TAKE INTO ACCOUf VALUES OF X AND Y LYING 
C ABOVE AND BELOU THE HAXIMUH AND MINIHUH VALUES IN THE ARRAY 
DIMENSION XARRAYI IDIH) , YARKAYI JDIM) , ARRAY! lOIH; J'DIM) 

X--XX 

Y=YY 

)f(X.GT.XTOP) X=XTOP 
IF(Y.GT.YrOP) Y=YTOP 
C RIND INDICES ON X ARRAY 
IDMi^IDIM-i 
DO 5 IMAX=2,IDMi 
If(X.LE.XARRAYaMAX)) GO TO 7 
S CONTINUE 

7 IMIN=IMAX-i 

C FIND THE INDICES OF THE VALUES IN THE YARRAY WHICH BRACKET THE 
C ACTUAL Y VALUE 

JDHl=-JDIH-i 
DU iS JliAX=2,JDMi 
inY. LE. YARRAY! JMAX)) GO TO 17 
IS CONriNUt 

17 JMIN=JMAX-i 

C INTERPOLATE TO FIND THE ARRAY VALUES CORRESPONDING TO THE 

C BRAi:KEriNG X VALUES AND THE ACTUAL Y VALUES 

ALOWY=! !X--XARRAY! IMIN) )*ARRAY!IHAX, JMIN)+ 

?,!XARRAY!IMAX)-X)*ARRAY!IHIN,JMIN))/!XARRAY!IMAX)-XARRAY!IMIN)) 

AHIGHY=! ! X-XARRAY! IMIM) )*AKRAY! IMAX , JMAX )+ 

.\!XARRAY!IMAX)-X)1iARRAY!IMIN,JMAX))/!XAfiRAY!IMAX)-XARRAY!IHIN)) 

C INTERPOLATE BETWEEN THE ABOVE. ARRAY VALUES TO FIND THE ARRAY VALUE 
C CUR8ESPONOING TO THE ACTUAL Y VALUE 

INTNG|^!!Y-YARRAY!JMIN))*AH1GHY+!YARKAY!JMAX)-Y)*AL0WY)/ 

! YARRAY ! JMAX ) -YARRAY! JI1IN ) > 

RETURN 

END 

LND 
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l.'l,AL FUNCTION INTEO(ARRAY,XX,YY..XTOP,YTQP, 

4XC0N). . XCON^ , YCONl . YCUH2 . IDIM . XDIH ) 

mmmm%mm$mtutnmnu%mmmmmmm*mnnmmu%m 

c 

C DOUBLE VARIABLE INTERPOLATION SUBROUTINE, BOTH VARIABLES MUST HAVE 
C EOUAL SHACINU BETyCEN THE DATA POIMIS. 

C 

cntmmm%mmmuumm%%%m%nmmum%%u%mumnmtmm 

DIMENSION ARRAY(IDIM,JDIM) 

C 

X=XX 

Y=YY 

;if (X.GT.XTOP> X=XTOP 
IT-'(Y.Cr.YrOP) Y-YTUP 
C X>FLOAT(I):i:X(:ONi+XCONB 
C Y-F LOA r (I ) aYCONl >YC0N2 
C 

C FIND THE INDICE OF THE X VALUE WHICH BRACKETS THE ACTUAL X VALUE 
It1IN=<X-XC0N2)/XC0Nl 
IFdMIN.LT.i.i IHlH=i 
lEdMIN.ijE. IDIH) IMI«=II)IH“1 
lMAX=IMINd 

X(1IN=FL0AT(IMIN)SXC0Ni4XC0N2 

XMAX^XMINvXCONi 

C FIND THE INDICE OF THE Y VALUE WHICH BRACKETS THE ACTUAL Y VALUE 
JM[N--(Y-YC0N2)/YC0N1 
)F(JMIN.LT.i) JMlN-1 
IF ( JM IN . GE . JDI M ) JHIN-- J DIM-i 
JHAX=JhINd 

YNlN=FL0AT(.JMIN)*YC0Ni+YC0N2 
YhAX=YMIIH YCONl 

C INTERPOLATE TO FIND THE ARRAY VALUE CORRESPONDING TO THE BRACKETING 
C X VALUES AND THE ACTUAL X VALUES 

ALOW Y= ( ( X-XIil N ) ?ARRA Y ( IMAX . JMIN ) -i 
A(XMAX-X)»ARRAY( IMIN,JMIN) )>(XMAX XMIN) 

AHIGHY=((X-XMIN)T;ARRAY(IMAX,JMAX)4 
i ( XMAX-X ) $ARKAY( IMIN , JHaX) ) /(XMAX- XI1IN ) 

C INTERPOLATE TO FIND THE ARRAY VALUE CORRESPONDING TO THE ACTUAL 
C Y VALUE 

INTEt!= ( ( Y- YMIN) *AHIGHY+ ( YMAX-Y)»ALOWY ) /( YMAX--YMIN ) 

RETURN 

END 

END 
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O O o 


'ifJBROUTIHC WGHT({I£;KW,WGBATT, WT,tir,WGCHI,WP) 

{:%n%n%%nn%uuu%uuutunni%*nu%nun*tunu%tntuutttutt%%nu 

THIS CALCULATES THE TOTAL MASS AND CURB MASS OF THE VEHICLE. 

cunuun%*ututuu*unnuuuuuu%uuuun%uuuunttttnut%ut* 

HEAL Mf,MTRCl,MrRKH,MTGRK« 

INTEGER RUN2I 
COMMON RUN2I 

COHMON/K ILOG/UGENG ,«GCVT , WGCHGR , WGGRM , WGGRDF 
4 , WG D I FF , UF , WP L , MPROP , HS , UGMO T8 , WGCLT , UCUR B 
COHMON/VARV/ WTL,INOCHG 

DATA CVTCi/i.07/, CVTKW/90.0/, HECVO.O/, HEC2/2.84/ 

DATA HTRCt/0.6/, MTRKU/20.0/, CHICi/0.9i/, CHIKW/40.0/ 

DATA MGEXTR/O.fl/ 

DATA DIFFCi/0.66/, DIFFKW/90 . 0/, GRCi/0.3/, MTGRKH/40.0/ 

DATA DFGRKU/90.0/, CHGRCi/0.0/, CHG8KU/14.S/ 

C COMPUTE CONSTANT TERM IN VEHICLE WEIGHT FORMULA 
C (FROM MT1',APPENDIX8) 

UCONST=( (0 .23«WPL+WF)/0 .77>+MTL 
C CALC. PROPULSION SYST. COMPONENT WEIGHTS 
C CVT 

HGCVT=CVTCi*CVTKW 

C HE 

WGENG=(HECH!:HEKW+HEC2)*HEKW 
C MOTOR 

HGM0TR=--MTRC1)|:MTRKU 
C CIIOPPER/INV. 

WGCHI=CHICiTCHIKW 
C DIFFT'L 

WGDIFF=DIFFCi«PIFFKH 
C GEAR BETWEEN MOTOR AND CVT 
WGGRH=GRCi*MTGRKW 
C GEAR BETWEEN DIFF t, CVT 
WGGRDF=GRC1*DFGRKW 
C CHARGER 

WGCHGR=CHGRCi«CHGRKW 
C CALCULATE PROPULSION SYSTEM WEIGHT 

WPROP=UGENG+WGCVT+WGMOTR+WGCHI+WGDIFF+WGGRH4WGGRDF+WGCHGR 
.i+WGEXTR+WGBATT!WGCLT 
C CALCULATE VEHICLE TOTAL WEIGHT 
Wr=HCONSr + WPROP/O.77 
C CALCUt.ATE VEHICLE TOTAL MASS 
MT=WT 

C CALCULATE CURB MASS OF VEHICLE 
WCURB=WT-WTL 

C CALC. POWER TRAIN WEIGHT 
WP-TWPROP-WG8ATT) 

A=WGCVT+WGGGRM+WGGRDF+WGDIFF 

IF(RUN21.LE.i) WRITE(6,200i) WCURB,WF>WGBATT, WGCHGR, WGCHI,WGCLT, 
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h I4GCVT , yCOIFK ,«GeNG ,«GEXTR ,UGGRDF , WGGRH , WGHO IR , UP ,UPL .UPltOP , 
4Wr,WTL,A 

2001 FORtiAK/'mmmm vehicle masses unnmn',/ 

CURB MASS OF VEHICLE <WCURB)=' ,F9.3,/ 
FIXED VEHICLE MASS (WF)=',F9.3,/ 
BATTERY HEIGHT (WGBATn=SF9.3,/ 
CHARGER MASS (WGCHGR)=',F9.3,/ 
CHOPPER INVERTER HASS (HGCHI)=',F9.3,/ 
CLUTCH MASS (WGCLT)=',F9.3,/ 
CVr HASS (HGCVT)=',F9.3,/ 
DIFFERENTIAL HASS (WGDIFF)=',F9.3,/ 
ENGINE HASS (UG£NG)=' ,F9.,5,/ 
i'EXTRA VEHICLE PROPULSION COMPONENT HASS (WGEXTR)=' F9.3,/ 
i'HASS OF GEAR BETHEEN CVT i DIFFERENTIAL <UGGRDF)=' ,F9.3,/ 

U' MASS OF GEAR BETWEEN MOTOR L CVT (WGGRM)=' ,F9.3,/ 

f'' ELECTRIC HOrOR HASS (WGHOTR)=',F9.3,/ 

POWER TRAIN HASS (HP)=^,F9.3,/ 
MAXIMUM PAYLOAD MASS (WPL)=' ,F9.3,/ 
PROPULSION SYSTEM MASS (WPR0P)=',F9.3,/ 
TOTAL HASS (WT)=^F9.3,/ 
lEST PAYLOAD MASS (HTL)=',F9.3,/ 
h ' UGCVT+WGGGRM fWGGRDF+WGO I FF= ' , / 

TRANSHISSION A DRIVE TRAIN WEIGHT=',F9.3) 

RETURN 

END 

END 


156 



SUPfiOUllME MrCST(«GBATT,«GCHI, MNFCST) 

c 

C llilS DETERMINES THE MANUFACTURING COST OF THE VEHICLE AS A FUNCTION 
C OF THE VARIOUS MASSES, 

C 

cuuunut%uunnu%tnt%nnntnuu%tuttuu%%tuutuu%ntu%%n%t%t 

REAL MNFCS f ,MTRCER ,«TRCST 
INTEGER RUNS! 

COMMON RUNS! 

COMHON/K IL OG/WGENG , WGCVT ,UGCHGR , UGGRM , UGGR DF 
, UGDIFF , Ur , URL , WPKOP , US ,UGMOTR « UGCL f ,UCUR B 
DATA HECER/I.4S/, CVfCER/i .42/, CHGCER/S.09/ 

DATA C'HICEH/lA.li/, GRCER/2.43/, DIFCER/2.42/, UFCER/3.13/ 

DATA WSCER/?./, MTfiCER/iS.9S/, CLTCER/0.0/ 

DATA ashi:er/().v;o4/ 

C CALCULATES MANUFACTURING COSTS FOR PROPULSION SYSTEM COMPONENTS 
C 

HECOST=HECER*WGENG 

CVTC3T=CVTC£R*WGCVT 

Ch'GCST^CHGCERifUGCHGR 

CHH:ST=CHICER!f.UGCllI 

GRMCST=GRCER!»WGGRM 

GRDCSr=GRCERTUGGRDF 

))ICCST=DIFCER)ltUGDIFr 

UFC03T=WFCER1;«F 

WS--=0.3*(UF'L+WPROP+MF) 

USC03T^H3CER»:UG 
mR(;3T=r.HTRCER$WGMOTR 
CLTC3I-CLTCERSUGCLT 
AGMCST=ASMCER»UCURE 
C CALCULA FE TOTAL MANUFACfURING COST 

A=CVrCST+GRMCST+GRDCST+DIFCST+CLTCST 
MNFCS r^^HfXOS T<CVTCST+CHGCS r^-CH^CSr+GfiHi:ST+GRDC3 T 
4iDlFCST+WFCQST+«SC03T+MTRCST+CLTCST+A3MCSI 
IF(RIJN21.LE.l) WRITF(6,i900.) US 
iVOO FORMAT! 

STRUCTURE 4 CHASSIS MASS <WS)=^' ,F9 .3) 

IF(RIJM 2 i .LE.i) URI rL(6,20C0) ASHCST,CHGCSr,CHICST,CLrCSr,CVTCSr, 
4DlFCST,GRDCST,GRMCST,HEC0ST,MNFC3T,MTfiCST,UFC0ST,USCQ3T,A 

2000 mmu/"mnnntn manufacturing costs follou nnmxnuv ,/ 

4' ASEMBLY COST (ASMCST)=SF9.3,/ 

y ClFARGER COST (CHGCSf)-' ,F9.3,/ 

i' CHOPPER INVERTER COST (CHICST)=' ,F9.3,/ 

A' CLUrCII COST (CLTCGf) = ',F?.3,/ 

4' CV1 COST (CVTCST)=',F?.3,/ 

4' DIFFERENTIAL COST (DIFCSD=' ,F9.3,/ 

4' GEAR BETWEEN DIFF 4 CVT COST <GRDCST)=',F9.3,/ 

4' GEAR BEfWEEN MOTOR 4 CVT COST (GRMCSn=' ,F9.3,/ 

4' HEAT ENGINE COST (HEC0ST)=',F9.3,/ 
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// fOrftL HAN!JFACfU«ING COST (HNFCST)=',F9.3,/ 

FLEriRIC hOrtlR COST (MTRCST)='',F9.3,/ 
fi' COST OF VEHICLF FIXED WEIGHT (UFCOST)=' ,F9.3,/ 
COST OF STRUCTURE 4 CHASSIS (WSCOST)=^' ,F9.3,/ 
4^ CUTCST+GR«CSr+GROCST>DIFCS T'rCL rC3T-',F9 .3) 

RETURN 
END 
END 
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^:UBROUTINE. HAlNT(HEK«,MTKWPi(,BATCST,KMPRYR ,UP ,WGCHI,BTLIST,HAINTC) 

cu%n%ntnu%nnnutnuutn%tn%%%nnt%nu%nnnn%%n*%tt'n%tttttt*t 

c 

C THIS CALCULATES THE YEARLY MAINTENANCE COST OF THE VEHICLE AS A FUNCTION 
C OF HEAT ENGINE, eATTERY, AND VEHICLE MAINTENANCE, AND ANNUAL DRIVING 
C DISTANCE. 

C 

CMtuutuu%utuuuuuuutuuutnn%iunuuuni%utuut*iiUtu%tt 

REAL K«CST,MTKWPK,KMPRYR,MAINTC,HTMAIN 

INTEGER RUH2I 

COMMON RUNEi 

COMHON/BATT/BTHKUP 

COMHON/COS r/G ASCST (i Q ) , KUC3T (iO > 

C CALCULATE COMPONENT MAINTENANCE COSTS <$/KM) 

HEM A IN= ( D . i O>0 . 0fl3»(H£KM/fl . 746 ) ) /I 60 . 73 
HTHAIN=(Q . 06+0 . OOESIMTKMPK/O .746) )/I6fl .93 
BrLISr=8ATC3T?(t . +BTMKUP ) 

BTMAIN=(0 . 0004*BTLIST)/160 .93 
CVrHAN=fl. 063/160. 93 

C (AR-HP. 78) i/2 HAINTANENCE i/2 REPAIR 

CHIMAN=((36.83/2)43qRT(UGCHI))/(ifl.*KMPRYR) 

C CALCULATE MAINTENANCE COST FOR ONE YEAR 

MAINTC=(H£HAIN+MTHAIN+8THAIH+CVTHAN+CHIMAN)^KMPRYR 
JF(RUN2i.LE.i) WRITE(6,2020) BTLIST 
2020 FORMAK' LIST COST OF BATTERIES (BTLISf)=' ,F9.3) 

IF(RUN2i.LE.i) MfiITE(6,2030) BTHAIN,CVTMAN,CHIMAN,HEHA1N, 
iMAINTCjMTHAIN 

2030 FORMATI/'ltSCWCWm MAINTENANCE COSTS FOLLOW UtUtUtUU*' ,/ 

6' BATTERY MAINTENANCE COST PER KILOMETER <BTHAIN)=' ,F9.8,/ 

tx' cvr MAINTENANCE COST PER KILOMETER <CVTMAN)=' ,F9.8,/ 

6'CHOPPER INVERTER MAINTENANCE COST PER KILOMETER (CHIMAN)=' ,F9.8,/ 

6' HEAT ENGINE MAINTENANCE COST PER KILOMETER (HEMAIN)=' ,F9 .8,/ 

■V TOTAL MAINTENANCE COST (MAINTC)-=',F9.3,/ 

6' ELECTRIC MOTOR MAINTENANCE COST PER KILOMETER <MTMAIN)=' ,F9.8) 

RETURN 

END 

END 
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%'BRQUTINE Rf.PAR(HEKU,MrKUPK,«P,WGCHI,KHPfiYR, REPAIR) 

cum%mimtmm%%m*%%nmmmttmuuunmntmtmnnmnm 

c 

C THIS CALCULATES THE YEARLY VfHICLE REPAIR COST AS A FUNCTION OF THE 
C ANTICIPArtO REPAIR COSTS FOR THE HEAT ENGINE, MOTOR, CVT, POWER TRAIN, 

C AND CHOPPER INVERTER. 

C 


cmmnmmmmumnnnmmmmmmmmmmmmmmuu 

REAL KHPRYRjMTKUPK 
INTEGER RUN21 
COMMON RUN2I 

REPRHE=(0 .28+0 . fl08*(HEKW/0 .746) )/I60 .93 
REPMfR--( 0 . 09+0 . fl02*(HTK«PK/0 . 746) ) /i 60 . 93 
REPCVT=(0 . 05+0 . fl0i3i(t(MTKWPK/fl .746) )/I6fl .93 
REPC1II= ( ( 36 . 33/2 . ) *30R f ( WGCHI ) )/ (i 0 . *KMPR YR ) 

C CALCULATE YEARLY REPAIR COST 

REPAIfi=(fiEPRHEtREPHTR+RePCVT+R£PCHI))):KMPRYR 


) r(RUN2i . LE. . I ) HRITE(6,2020 ) REPAIR ,REPCHI ,REPCVT,REPMTR ,REPRKE 

2021 ) FORMAr(/'4!m**m«*m« repair costs follow tutnttutnt' ,/ 

TOTAL REPAIR COST (REPAIR)=',F9.3,/ 
6'CHOPPER INVERTER REPAIR COST PER KILOMETER (REPCHI )=' ,F9.8,/ 

CVT REPAIR COST PER KILOMETER (REPCVT)=',F9.8,/ 

6' ELECTRIC MOTOR REPAIR COST PER KILOMETER (REPMTR)=" ,F?.0,/ 

HE REPAIR COST PER KILOMETER (REPRHE)=',F9.8) 

RETURN 


END 

END 
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'3UIiR0UTIN[-; A6CST<HNFC3T,WGBATT, BATCST,ACQCST) 

c 

C THIS DETFRhlNSS THE ACQUISIUON COST OF THE VEHICLE AS A FUNCTION OF THE 
C MANUFACrUHING COST, BATTERY MANUFACTURING COST, AND THE MARKUP FACTOR. 

C 

REAL MARKUP, MNFCST 
INTEGER RUH21 
COMMON RUN2I 
COMMON/BATT/BTMKUP 
DATA BATCER/1.87/ 

C CALCIJI.ACE THE MARKUP FACTOR 

MARKUP=(29 . i76E.-5)*MNFCST+i .40833 
C CAl.CULiUE BATTEtTY COST 
.()ATCST=BATCER*MGBATT 

C CALCULATE ACOUISITIOH COST (COST TO CONSUMER) 
ACQCST=MARKUP?.MNFCST+BATCST$(i . +BTMKUP ) 
lF(RUt«l.LE.i) WRlTETAjiUl) BATCST, MARKUP 
iOv lORMATC 

Tr-HANUFACTIIRING COST OF BATTERIES (BATCST)=' ,F?.3,/ 

A40X/Ue 00 W. DEALER MARKUP FACTOR (MARKUP)^' ,Gii .b) 

RETURN 

END 

(,:nd 
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F1H4,L 

SIJDROUTIME UFE(PETYR,ELFCYR,ELEC,DISCH,BTLIST,ACQCST,«A]HTC, 
/.RD’AER,DAfCST,LCC03T,Hf:TCST,EU-;CST,TlJrPEn 


c 

C ims CALCULATES THE LIFE CYCLE COST OF THE VEHICLE BASED ON TEN YEARS 
C USAGi-:. IT IS A FUNCTION OF THE ACQUISITIUN, MAINTENANCE, REPAIR, 

C BATTERY, FUEL, AND ELECTRICITY COSTS OVER THAT PERIOD. 

C 

mtmmmnmtt%nm*mnnmmm%*tmmunmutmmmnmn 

LOGICAL SMOOTH 
REAL KWCST,MAINTC,LCCOST 
INTEGER RUNPi 
COMMON RUN2T 
COMMON GT 

CUHMON/BATT/DTMKUP 
CGMH0N/C03T/GA3C3T<ifl ) jXWCSmO ) 

CaHMON/LIHIT/NTfiIPX,nERMX 
CO MMON /OM 1 3C/CH00SE , 3MOO TH 
COMMON BTC03T,BATRYR,BATLFE 

DIHEN3I0N PE f YR US) ,£LECY« ( 15) , ELEC US ) , DISCHU5 ) , 

J.PETCSTUfl),ELEC3TUfl) 

DIMENSION VHC03TUO) 

DATA DISCNT/fl . 02/,DISLFE/84() ./ 

C CALC PETROL AND EcEC COSTS-. INITIALIZE YEAR INDEX AND TOTAL COSTS 
DO 5 lYEAR^ljlO 
PETC3raYEAR)=D. 

S ELECST(IYEAR)=--fl. 

DO 10 IYEAR=i,lQ 
DO 10 JTR1P=1,NTRIPX 

C ADD COSTS FOR mi3 TRIP TO TOTAL COSTS 

PLTCST ( I YEAR )=PETCST(IYEAR ) +PE TYR ( JTR IP > *GASCST< I YEAR )*GF 
10 ELECST ( lYEAR )=aEC3TI lYEAK )+(CHOOSE«ELECVR ( J'TRIP ) + 

4(1. -CHOOSE )i!(£LEC( JTR IP ) ) *KHCST( lYEAR ) 

C CALC NO. OF DEEP DISCHARGES OF BATT. DURING ONE YEAR 
DEEP^O , 

DO 20 JTRIP=1,NTRIPX 
20 0EEP=DEEP+0I3CH(JTRIP) 

C CALC NO. OF YEARS IN DATT LIFE 
3ATLF£^DISLFE/D££P 
C CALC BATT REPLACEMENT COST 

B rCOST-BATCSTS U . -p'2 . SB TMXUP ) 

C CALC TOTAL VEHICLE COST FOR ONE YEAR 
COlJHV-1. 

TOTBAT^O . 

BATRYR=--( ( UO ./BATLFE)-! . )*BTCflST)/10 . 

DU 40 IYEAR-1,10 

VHC03T (I YEAR )=PETCST( lYEAR )-»MAINTC+ELECST ( I YE AR )+REPAIR 

IF(. NOT. SHOO ni) GOTO 3S 

VHCOST( 1 YEAR )=VHCOST (lYEAR )-HiATRYR 
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TO r8AT=ra fBAr+B(!iT«Y«/<DiscN rn . )»*aYi-:AR-i ) 

GO TO 40 

C IS THIS A YEAR FOR BATT REPLACEHEMT 
3S IF(FLOAraYEAR).LE.COUNTMArLFE) GOTO 40 
C ADD IN BATT REPLACEMENT COST 

VHCOS r ( I YEAR > =VHCOS I' ( I YEAR ) +0 . 9»BTC0S T 
TOTBAT=TOTBAT+BTCOST/( (DISCNT+1 . )**(IYEAR-1 ) ) 

COUNT=COUNT+i. 

40 CONTINUE 

C CALC VEHICLE AND POWER TRAIN AND BATT SALVAGE VALUE 
3ALVGE=0 . LiKACtlCSr-BfLIST) 
l!ATSAL=fl .5*BTLIST»(COUNT-(iO ./BATLFE) ) 

IFISMOOFH) BAfSAL^^O. 

C CALC LIFE CYCLE COST 
LCCOST=0 , 

TOTPET=Q. 

DO SO IYEAR=i,iO 

T0TPET=T0TPET+PETC3Ta YEAR )/( (DISCNT+i . )#*(IYEAR-i ) ) 

SO LCC0Sr=LCC0Sr+VHC03TaYEAR)/((DI3CNT+i . ) W(IYEAR-i ) ) 

C SUBTRACT SALVAGE VALUES 

LCCOST=LCCOST-(SALVGE+BATSAL)/((DISCNT+i , )«*i0) 

C ADO IN AOUISrriON COST 
LCC03T=LCC0ST+AC0CST 

IF(RUN2i.LE.i) WRITE(6,13i) BATLFE,BATSAL,LCCOST,BTCOST,SALVGE 
131 FORMAT! 



BATTERY LIFE (BATLFE) 

=',F9.3,/ 

4' 

BAT TER Y SALVAGE VALUE (8ATSAL) 

-',F9.3,/ 

i' 

LIFE CYCLE COST (LCCOST) 

=',F9.3,/ 

^ / 
•'I 

BAfTERY REPLACEMENT COST (BTCOST) 

=',F9.3,/ 

i/ 

SALVAGE VALUE (SALVGE) 

=^',F9.3) 


lF(RUN2i.LE.l) 

«.MRrTE(6,404) (I,FETCST(l>,ELECST(I),VHCOSTa),I=l,lfl) 
404 FtlRMAl !' ' ,/lX,'Y£Afi'' ,SX/PETROLEUM C0ST'',3X, 
i' ELECTRICITY COST ',3X/ VEHICLE COST',/ 

•V I ' , lOX /PFTC3 I'd) ' ,3X, 'ELECSTd ) ' ,3X, 'VHCOS 1(1)', 
&/iO(lX,I2,2X,3(8X,F9.5)/)) 

RE i URN 

END 

END 
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BLOCK DATA 

znmummmmmmnmmmmmmmmnnm 

c 

C BLOCK DATA 

c 

cmmmuuuuutmmmummnmmnmmum 

UJOCCAL 3H0UTH 

kl:al kwcsi,htkhpk,ivestc 

COHHON/ALPHA/ALPHAi , ALPHA2, ALPHA3 
CO«hON/BATT/BTKI(U? 

CUM(1UH/CC/C0A,R0TFF,KGfc;A« ,RUf!l1,R riRf 
COHrtON/COSl /GASCST (i 0 ) , KWCS I (1 0 ) 

COrlliON/tFC V f /Ai' CV T (3 ) , AROV F ( 7 ) , £FF < 8 , 7 > 
t:0«M0N/K I1.0G/MGENG , HGCVT , WGCHGR , WGG8 H , WGGR DF 
i , WGDlTF , UF , WPL , WPROP , W3 , WGMOTR , WGCL F , UCUR B 
COhHON/LrMlT/HTRIPXjITERHX 
DJi'iKON/rtAP/DISCIiG(ifli ,6> ,EFFHTR(8l,^ ,»'r:iiAP(6,&) , 
iCREFP (ifli ,li > ,VCHG(i3,i2) ,AT0R0(6) ,AESPED(6) jARATIOdO) 
COMMON/OfFIGC/CI IGiiSE , 3«00 ( 1 1 
COHAON/hOTOR/firK'WPK 
cohNUH/PHr/Gi?Ay,ani] 

COh'HON/P R I HT/APETT( 1 S ) , AELETA (iS ) , AEL.ETP ( i S ) 

Ci]H«03/RRATa/f'i.riil£,yGhAR,aCVr,PllLlTaR,PtlTfirD,PGtAR 

COHMON/RSPt.T/ACCfl,SOC,PCVT,CX(iS),VAVE,PHb.MAX, PH£,PELEC 

COHAON/SCAL£/SFX( 17) ,SFG(28 ) ,NN 

COHMON/S TC/y ( 031 ) , BET ( 831 ) , TR IPLN ( 8 ) , TR IP NO ( 8 ) 

COilMtJN/CHAHG/AAOVFOljACRlB) 

CGKMON/TIMEC/PLTSTC , DELTHW , TACCEL , TCOAST , TBR AK£ , TSTOP , TCYCLE 
CUHNOH/ F fi IP /CHA«G£(36 ) , EB FO F L ( 36 ) , HOU F ( 36 > , I DOD ( 36 ) , 
6)VE8TC(36),PETSTC(36),S0CSTC(36),REGEN(36) 

CGArtON/VARY/ OFL^INOCHG 

C;OHfiON/Y£ARC/PErYR(15),ELECYR(i5),ELEC(lS),DISCH(lS) 

C 

0 ALPHA 

DATA ALPHAl/I./, ALPHA2/0./, ALPHA3/0V 
C BAFT 

DATA BTMKUP/.3/ 

i: CC 

DATA CDA/.6/, RDIFF/1./, RGEAR/. 08333333/, RGRH/0.2857/ 

DA FA R FIRE/. 3/ 

C COSF 

DATA CASCST/C. 32,0. 35, 0.38,0. 41, 0.44, 0.47, 0.5,0. S3, 0.56,0.59/ 
DATA KWCSF/.fllll, .0125, .0139, .0153, .0167, .0181, .0194, 
6.0208, .0222, .0236/ 

C KILOU 

DATA WrVSlO./, MGCLT/fl./, WPL/41S./ 

C LIMIT 

DATA ITERHX/5/, NTR IPX/8/ 

C HO FOR 

DATA HTKWPK740./ 
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c OiiiGi: 

DATA CHOOSt/C . /, SMOOTH/ . TRUE . / 

c pur 

DATA GRAV/9.8U7/, RHO/i.?.fS/ 

0 RiTATU 

DATA PHTRTD/tIfiCOO./ 

C OCALf 

DATA SFX/20<i./, NN/S/ 

DA (A 0rG/2t)'*l./ 

c ore 

DATA PET/830X0./ 

DATA TRIPLN/10C00.,?000C. ,50006. ,80000. ,130000. ,160000 
6500000. ,01)0000./ 

DATA TRrPK0/t3C.,85.,57.,54.,l2.,?.,3.,i./ 

C (E3T 

DATA AACV/2. 265, 2. 2225,2.1993,?. 1827, 2. 1702, 2, 1607, 2. 1524, 2. 1455/ 
DA TA ACR /O . , 1 . 0 1 7 , 2 . 0-34 , 3 . 052 , 4 . 069 , 5 . 0 86 . 6 . 1 04 , 7 . 1 21 / 
c 

DATA DLTirit/l./, DElTHW/3'0 . / 

DATA 1ACCEL/14.C/, TCOAST/64.0/, TDRAKE/74.0/, TSTOP/8.3./ 

DATA rCYCLE/lOU./ 

C DA'/r 

DATA WTL/2G7./, INOCHG/1/ 

1 ) 

DAT A DISCHG/1 .30,1.30,1.30,1.30,1.30,1.30,1.30,1.30, 

61. 3000, 1.3000, 1.3000, 1.3000,1. .3000, 1.3000, 1.3000, 1,30(10, 

<’-j .3000,1.3000,1.3000,1.3000,1.3000,1.6800,1.7400,1.7570, 
61.76/0,1.7800,1.7930,1.3015,1.8100,1.8198,1. 8295 , 1 . 8393 , 

6 5 . 8490 , 1 . 8565 , 1 . 8640 ,1.8715,1. 8790 , i . 8345 ,1.8900,1. 8955, 
61.9010,1.9048,1. 9085 ,1.9123,1. 9160 ,1.9197,1. 9235,1 . 9272, 

6 1 .9310,1. 9349 , 1 . 9388 , 1 . 9426 , 1 . 9465 ,1.9504,1. 9542 , 1 . 9581 , 

61 . 9620 , 1 . 9645 , 1 . 9670 , 1 . 9695,1 . 9720 , 1 . 9745 ,1.9770,1. 9795, 

61 . 9820 , 1 . 9845,1 . 9870 , 1 . 9895 , 1 . 9920 , 1 . 9942 ,1.9964,1. 9986 , 
62.0008,2.0030,2. 0052,2 . 0074 ,2.0096,2.0118,2.0140,2. 0159 , 
62.0178,2.0197,2.0217,2.0236,2.0255,2.0274,2.0293,2.0313, 

62 . 03.32 , 2 . 0351 ,2 . 037J ,2 . 0.502 ,2 . 0394 , 2 . 0406 , 2 . 0418 , 2 . 0430 , 

62 . 0442,2. 0454,2 . 0466,2 . 0478,2 . 0490 , 

61 .3000,1 .3000,1.3000,1.3000,1.3000,1.3000,1 .3000,1.3003, 

61 .3000,1.3000,1.3000,1.3000,1.3000,1.3000,1.3000,1.3000, 

61 . 30 00 , 1 . 30 00 , 1 . 3000 , 1 . 3000 ,1 . 3000 ,1 . 5200 , 1 . 6320 , 1 . 63.70 , 
61.7190,1.7420,1.7570,1. 7685 ,1.7800,1. 7883 , 1 . 7965 >1.8047, 

61.8130,1. 8185 , 1 . 824 9,1. 8295 , 1 . 8350 , 1 . 8395 , 1 . 8440 , 1 . 8485 , 

61 . 8530 , 1 . 8565 ,1.8600,1. 8635 , 1 . 3670 ,1.8705,1.8740,1. 8775 , 

61.8810,1. 0834 , 1 . 0058 , 1 . 0081 ,1.8905,1. 0929 ,1.0953,1. 0;76 , 
61.9000,1.9020,1 .9040,1.9060,1.9080,1.9100,1.9120,1.9140, 
61.9160,1.9180,1 .9200, 1.9220, 1.9240, 1.9260, 1.9280, 1.9300, 

61 . 9320 , 1 . 9340 , 1 . 9360 , 1 . 9380 ,1.9400,1. 9420 , 1 . 9440 , 1 . 9453 , 

61 .9465,1 .9478,1 .9490,1 .9503,1 .9515,1 .9528,1 .9540,1 .9553, 
61.9565,1.9578,1.9590,1. 9598 ,1.9606,1.9614,1.9622,1.9630, 

61 . 9638 ,1.9646,1. 9654 , 1 . 9662 , 1 . 9670 , 
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41.3000,1 .3000, 1.3000, 1.30004-3000, 1.3000, 1,3000,1.3000, 

41.3000.1 3000,1.3000,1.3300,1 ,3000,1.3000,1.3009,1.3000, 
41 ,30CG,1.3000,i.3GOO,1.3oGO,1..300C,l .4iJ00, 1.1.530,1. 60U0, 
41 . 6-5'fOO , 1 . 661 0 , 1 . 6V50 , 1 . 6855 , 1 . 69.50 ,1.7025,1. 7090 , 1 . 7155 
4 1 . 7220 , 1 . 7273 , 1 . 7325 ,1 . 7378 , 1 . 7430 ,1.7473,1.7515,1. 7557 , 
41 . 750 0 , i . 7625 , 1 . 7650 , 1 . 76?5 ,1.7700,1 7725 , 1 . 7753 , 1 . 7/75 , 
4 1 . 780 0,1. 7825 , 1 . 7850 , 1 . 7875 4 . 790 0 , 1 . 7925 , 1 . 7956 , 1 . 7975 , 

41. 1. ! 000 , 1 . 8017 , 1 . 8035 , 1 . 8052,1. 8070, 1 .UO 884 . 8105, 1.8123, 
41 .8140,1.8150,1.8175,1.8193,1.8210,1.8229,1.8248,1.8267, 
■U . 8285 , 1 8-505 , 1 . 8-524 , 1 . 8343 , 1 . 8352 , 1 . 8,581 ,1.8400,1. 8408 , 
41.8417,1.8425,1.8433,1.04424.8450,1.8458,1,8467,1.8475, 
41 , 848-5 , 1 . 8472 , \ . 850 0,1. 051 0,1. 8520 ,1.8530,1. 8540 , 1 . 8S5.i , 
4 5 , 8560 , 1 . 8S7G , 1 . 8580 ,1.8590,1. 060 0 , 

41.3000,1.3000,1.3000,1.3000,1.3000,1 .3000,5. .-5000,1.3000, 
41.3000,1.3000,1.3000,1 .30 00 >1.30 06, 1.3000, 1.3000,1.3000, 
41 . 30 0 0 , 1 . 30 0 0 , 1 . 30 0 0 , 1 . 30 0 0 , 1 . 30 0 0 , 1 . 41 3 0 , 1 . 460 0 , 1 . 433 0 , 
45 . 5030 , 1 . 521 0 , 1 . 5370 , 1 . 5405 ,1.5600,1.5698,1. 5795 , i . 5893 , 
41 . 577 0,1.6050,1.540,1.6170,1. 6230 , 1 . 6277 , 1 . 6325 , 1 . 6372 , 
4 5 .6420 ,1 . 6463,1 . 6505,1 .6547,1 .6590 , 1 . 6633,1 .6675,1 .6718, 
41 . 5/60 , 1 . £'770 , 1 . 6820 , 1 . 6850 , 1 . 6880 , 1 . 6710 ,1 .6940,1. 67 7 0 , 
45 .7600,1.7026,5 .7052,1 .7677 >1 .7103,1 .7129, 1 .7155,1 .7131 , 
41 . 7207 , 1 . 7233 , 1 . 7258 , 1 . 728-1 , 1 . 731 0,1. 7326 , 1 7342 , 1 7358 , 
4 5 . 7374 , 1 . 7390 ,1.7406,1. 7422 , 5 . 7438 , 1 . 7454 ,1.7470,1. 7486 , 
41 .7490,1. 7500 , 1 . 7510 , 1 . 7520 , 1 , 7530 ,1.7540,1. 7550 , i . 7560 , 
4 5 . 757 0,1. 758G ,1 . 7590 , 1 . 7598 ,1.7606,1.7614,1. 7622 , 1 . 7630 , 
4 1 . 7538 ,1.7645,1. 7554 , i . 7652 , 1 , 767 G , 

45 .3060.1.3066.1.3000.1.3066.1.3000.1.3006.1.3000.1.3000, 

41.. 10JG, 1.3000, 1.-3000, 1.3000, 1.3000, 1.3000, 1.3000, 1.3000, 
45 .3006,1.3006,1.3000,1.3000,1.3006,1.3670,1.4000,1.4186, 
41 . 433G , 1 . 4530 ,1 . 4660 , 1 . 4755,1 , 4870 , 1 . 4935,1 . 5000 , 1 .5055, 
45 .5130,1 .51 86,1. 523C- , 1 . 5286 , 1 . 5330 , 1 . 5380 ,1.5430,1. 5480 , 
41 . 5530 , 1 . 5571 , 1 . 5612 , 1 , 5654 , 1 . 5695 , 1 . 57.16 , 1 . 5/78, 1 . 5319 , 
4 i . 5860 , 1 . 5894 , 1 . 5927 , 1 . 5961 , 1 . 5995 ,1.6029,1.6062,1. 6096 , 
41.6130,1.6154,1. 6173 , 1 . 6202 , 1 . 6227 , 1 . 6251 , 1 . 6275 , 1 . 6299 , 
41 . 6323 , 1 . 6347,1 . 6372,1 . 6396 ,1.6420,1. 6445 ,1.6470,1. 6495 , 
41 . 652 J ,1 . 6545 , 1 . 6570 , 1 . 6595,1 . 6620 , 1 . 6645 , 1 .6670,1. 6686 , 
4 1 . 6702 ,1.6718,1. 6733 ,1.6749,1. 6765 ,1.6781,1. 6797 ,1.6813, 

41 . 6828 . 1 6844 , 1 . 6860 , 1 . 6878 , 1 . 6896 , 1 . 6914 , 1 . 69-32 , 1 . 6950 , 
41 .6968,1.6986,1.7004,1.7022,1.7040, 

41.3000.1.3000.1.3000.1.3000.1.3000.1.3000.1.3000.1.3000, 
41 3000.1.3606,1.3000,1 .3000,1.3060,1.3000,1.3060,1.3000, 

41.. 30 JO, 1.3000, 1.3000,1. 3000, 1.3000,1. 3160, 1.3310, 1.3430, 
41.3540,1. 3666 . 1 . 3720 , 1 .381 0 , 1 .3900,1. 3965 ,1.4030,1. 4095, 

4 1 . 4 1 60 . 1. . 42 05 , 1 . 4250 , 1 . 4295 , 1 434 0,1. 4388 , 1 . 4 435 , 1 4483 , 

41 . 4530 . 1 . 4565.1 .4660.1. 4635 . 1 . 4670 . 1. . 4705 . 1 . 4740 . 1 . 4775 , 

41 . 4810 . 1 . 48.39 . 1 . 4868 . 1 . -4896 ,1.4925,1. 4954 , 1 . 4983 , 1 . 501 1 , 
41.5040,1.5064,1.5088,1.5112,1.5137,1.5161,1.5185,1.5209, 
41 . 5233 , 1 . 5258 , 1 , 5282,1 . 5306 , 1 . 5330 , 1 . 5354 , 1 . 5373 , 1 . 5482 , 
4 1 . 5426 ,1.5450,1. 5474 ,1.5498,1. 5522 ,1.5546,1. 5570 ,1.5591, 
41 . 5612,1 . 5633,1 . 5653,1 . 5674 , i . 5695 , 1 . 57l 6 , 1 . 57.57 , 1 . 5758, 
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. 5778 , 1 . S7V9 > i . 5820 , i . 5845, i . 5870 , i . 5895 ,1.5920,1. 5945, 

. 5970 , i . 5995,1 . 6020 , 1 .6045, 1 . 6070/ 

D A1 A ATORQ/(i . , 19 . ,47 . , 94 . , 140 . , 178 . / 

OAT A AES7EO/0 .,209.419, 261 . 799 ,314.159,366.519,413. 879/ 

DATA EFFrtlR/O . , .8094, .91, .9175, .9250, .9280, .9310, .9315, .9320, 

6.9295. . 9270. .9240. .9210. .9175. .9140. .9035. .9030. .8780, 

6 . 8930 . . 8860 . . 8790 . . 8720 . . 8650 . . 8580 . . 851 0 . . 8442 . . 8375, 

6.3303. . 8240. .8170. .8100. .8030. .7960. .7890. .7820. .7750, 

6 . 7680 . . 761 0 . . 7540 . . 7470 ..7400.. 7332 . . 7265 .. 7198 .. 7130 , 

6.7060. . 6990. .6920. .6350. .6780. .6710 . .6640. .6570. .6500, 

6 . 6430 . . 6360 . . 6290 . . 6222. . 6155. . 6087 . . 6020 . . 5950 . . 5880 , 

6 . 531 0 . . 5740 . . 5670 . . 560, 0 , . 5530 , . 5460 , . 5390 , . 5320 , . 5250 , 

6.5180 .. 5112. .5045. .4977. .4910. .4840. .4770. .4700. .4630, 

it 

60.000. .7761. .8870. .8983. .9095. .9158. .9220. .9253. .9285, 

6 . 9293 . . 930, 0 , . 7295 , . 929 0 , . 9275 ,. 926 0 ,. 923 0 ,. 92 0 0 ,. 91 73 , 

6.9145. . 9108. .9070. .9026. .8983. .8939. .8895. .8846. .8798, 

6 . 874? ,.8700,. 3643 , . 3595 , . 3543 ,. 3490 ,. 8438 , . 3385 , . 8333 , 

6 . 8280 . . 8227 . . 81 75. . 8123. . 8070 .. 8020 .. 7970 .. 7920 . . 7870 , 

6.7317. . 7765. .7713. .7660. .7607. .7555. .7502. .7450. .7393, 

6 . 7345 . . 7292 .. 7240 .. 7189 . . 7138 .. 7086 .. 7035 . . 6982 . . 6930 , 

6.6377. . 6325. .6774. .6723. .6671. .6620. .6563. . 6515. .6463, 

6 .6410. . 6359. .6308. .6256. .6205. .6152. .6100. .6048. .5995, 

6 

60 .000. .7428. .8640. .8790. .8940. .9035. .9130. .9190. .9250, 

6 . 9290 . . 9330 . . 9350 . .9370 . .9375. . 9330 . . 9375. . 9370 . . 9365 , 

6.9360.. 9355 . . 9350 . . 9333 . . 931 5 . . 9298 . . 9280 . . 9250 . . 9220 , 

6. 9190. . 9160. .9125. .9090. .9055. .9020. .8935. .3950. .8915, 

6.8080. . 8845. .8810. .8775. .3,40, .8707, .8675, .8642, .8610, 

6,85V5, .8540, .3505, .3470, .3435, .8400, .8365, .3330, .8295, 

6. 8260.. 8225.. 8190.. 8155.. 8120.. 8085.. 8050.. 8015.. 7980, 

6.7945. . 7910. .7373. .7845. .7813. .7780.. 7745. .7710. .7675, 

6 . 7640 . . 7605. .7570 . .7535. .7500 . .7465. .7430 . . 7395. . 7360 , 

6 

6 0 . 0 0 0 , . 6862 , . 81 65 , . 8400 , . 8635 , . 8765 , . 8895 ,. 8980 ,. 9065 , 

6.9130. . 9195. .9230.. 9265. .9288. .9310. .9323. .9335. .9343, 

6.9350.. 9350 . . 9350 . . 9344 . . 9338 . . 9331 .. 9325 .. 9309 . . 9293 , 

6 . 9276 .. 9260 .. 9240 . . 9220 . . 9200 .. 9180 . . 9156 . . 9133 . . 9109 , 

6.9085. . 9056. .9028. .8999. .89,0 , .8943, .8915, .8888, .8860, 

6.8830. . 8800. .8770. .8740. .8709. .8678.. 8646. .8615. .8585, 

6 . 8555 . . 8525 . . 8495 . . 8464 . . 8433 ..8401.. 8370 . . 8340 . . 831 0 , 

6 . 8280 . . 8250 . . 8220 . . 8190 . . 8160 . . 8130 . . 81 0 0 , . 8070 , . 804 0 , 

6.8010. . 7979. .7948. .7916. .7885. .7855. .7825,'. 7795, .7765, 

it 

60 . 0 0 fl , . 6295 , . 7690 ,. 801 G ,. 8330 , . 8495 ,. 866(1 , . 8770 , . 8880 , 

6 . 8970 . , 9060 , .9110 , .9160 , . 9200 , . 9240 , , 9270 , . 9.300 , . 9320 , 

6 . 9340 . . 9345 . . 9350 . . 9355 . . 9360 . . 9365 .. 9370 .. 9368 . . 9365 , 

6.9362. . 9360. .9355. .9350. .9345. .9.340. .9328. .9315. .9302, 

6.9290. . 9268. .9245. .9222. .9200 . .9173. .9155. .9132. .9110, 

6.9085. . 9060. .9035, 9010, .8933, .8955, .8928, .8900, .88/5, 


167 



6 . 8550 . . 8828 . . 88,0 , . 8773 , . 8748, . 8718 , . 8690 , . 866S, . 8640 , 

6.8618. . 8870. .8863. .8838. .8808. .8480. .8458. .8430 . .8408, 

6 . 8380 . . 8353 . . 8325 . . 8298 . . 8270 . . 8245 .. 8220 .. 8195 .. 8170 , 

4 

6 0 . 0 0 0 , . 6295 , . 769 0 , . 80 1 0 , . 8330 , . 8495 ,. 8660 ,. 8755 ,. 8850 , 

6.8920. . 8790. .7035. .9080. .9115. .9150. .9180. .9210. .9225, 

6.9240.. 9255 . . 9270 . . 9280 . . 9290 . . 930, 0 , . 931 0 , . 931 0 , . 931 0 , 

6.9510. . 9310. .9310. .9310. .9310. .9310.. 9305. .9300. .9295, 

6 . 9290 . . 9282 . . 9275 .. 9267 .. 9260 .. 9247 . . 9235 . . 9222 . . 921 0 , 

6.9197. . 9185. .7173. .7160. .7148. .7135. .7122. .9110. .9095, 

6.9080. . 9065. .9050. .9035. .9020. .9005. .8990. .8975. .8960, 

6.8745. . 8730. .8717.. 8705. .8892. .8880. .8865. .8850. .8835, 

6 . 8820 .. 8805 .. 8790 . . 8775 . . 8760 . . 8745 . . 8730 . . 871 5 . . 870, 0 , 

6 0 . 0 0 0 , . 4946 , . 6520 ,. 7055 ,. 7590 , . 7845 , . 81 0 C , . 8245 , . 8390 , 

6 .8530. . 8670. .8745. .3820. .8870. .8960 . .9010. .9060. .9095, 

6 .9130, .9160 , .9190, .9210, .9230, .9250, .9270, .9280, .9290, 

6.9300. . 7310. .9315.. 9320. .9325. .9330. .9332. .7335. .7338, 

6 .9340. . 9343. .9345. .9347. .9350 . .9347. .9345. .9342. .9340, 

6.9332. . 9325. .9317. .9.310. .9302. .9295. .9287. .9280. .9267, 

6 .9255. . 9242. .9230. .9217. .9205. .9192. .9180 . .9162. .9145, 

6 . 9127 . . 911 0 . .9090 . .9070 . .9050 . . 9030 . .901,l , . 8990 , , 8970 , 

6 .8950. . 8927. .8905. .8382. .8860. .8835. .8310. .8785. .3760, 
i)f'i FA AH CV (70 . , 2 . 85,4 . 42 ,9 . 32 ,60 . ,80 . ,90 . , 10 0 . / 
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.97.59, .9742, .9746, .9750, ,9?53, 

. 9771 ,. 9773 ,. 9776 , . 9779 , . 9781 , 

. 9793 , . 9/95 , . 9797 , . 9798 ,.9800, 
.9809,. 9810,. 9312,. 9313,. 9815, 

. 9821 , . 9821 , . 9322, . 9823, . 9824 , 

. 9327 ,. 9827 ,. 9828 ,. 9828 ,. 9829 , 
.9833,. 9833,. 9834,. 9335,. 9835, 
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1. 9311. . 9812. .9813. .9814. .9815. .9816. .9817. .9818. .9819. .9820, 

4 . 9821 . . 9821 . . 9822 . . 9822 .. 9G23 ,. 9823 ,. 9823 ,. 9824 , . 9824 , . 9325 , 
6.9825/ 

DATA HtrtAP/0.,C.,C.,0.,0.,0., 

40.. 9.7141-;-(}4,i.371E-03,2.17ie-03,3.15£-03,3.943E-03, 
4(l.,1.2E-G3,1.714t-C3,2.686E-fl3,3.966E-03,5.038E-03, 
«.0.,1.436E-03,2.0S7E-8.5,3.267E-03,4.871£-03,6.16E-03, 

60.. 1.714E-03,2.366E-03,3.9ilE-fl3,S.829E-fl3,7.27S£-03, 

40.. 1.929e-33,2.686E-03,4.714E-03,6.857E-i)3,8.87;5E-03/ 

DATA VCHG/2. 100,2. 154,2. 134, 2. 121, 2. 114, 2. 104, 

62. 098, 2. 088, 2. 079,2. 069, 2. 057, 2. 045,2. 030, 

62.188,2.171,2.156,2.144,2.134,2.128, 

62.118,2.109,2.099,2.089,2.079,2.067,2.054, 

62.256,2.225,2.204,2.189,2.179,2.171, 

62.160,2.152,2.142,2.132,2.122,2.112,2.101, 

62.324,2.279,2.251,2.234,2.224,2.214, 

62.202,2.194,2.185,2.175,2.165,2.153,2.148, 

62.421,2.365,2.323,2.295,2.276,2.260, 

62.246,2.236,2.226,2.217,2.209,2.205,2.198, 

62.519,2.452,2.394,2.356,2.327,2.306, 

62.290,2.278,2.267,2. 259 , 2 . 254 ,2.251,2.247, 
62.617,2.539,2.466,2.416,2.379,2.352, 
62.335,2.320,2.308,2.301,2.293,2.297,2.297, 
62.715,2.626,2.538,2.477,2.431,2.397, 

62 . 379 , 2 . 362 , 2 . 349 , 2 . 343 , 2 . 343 , 2 . 344 , 2 . 346 , 
62.788,2.719,2.645,2.582,2.525,2.480, 
62.450,2.429,2.415,2.409,2.407,2.408,2.411, 
62.861,2.812,2.753,2.688,2.619,2.563, 

62 . 522 ,2.496,2.482,2.474,2. 4/1 , 2 . 472 , 2 . 4 75 , 

6 2 . 9 1 0 , 2 . 8 7 3 , 2 . 8 23 , 2 . 7 62 , 2 . 6 9 7 , 2 . 6 4 3 , 

62 . 599 , 2 . 567 , 2 . 544 , 2 . 526 , 2 . 51 7 , 2 . 512 , 2 . 5 1 S , 
62.959,2.934,2.892,2.837,2. 775 ,2.722, 

62 . 676 , 2 . 638 , 2 . 605 , 2 . 573 , 2 . 563 , 2 . 553 , 2 . 554 / 

END 

EDO 
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